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Measurements are reported of the proton and deuteron spin structure functions 
and gf at beam energies of 29.1, 16.2, and 9.7 GeV, and gff ^nd (7^ at a beam energy 
of 29.1 GeV. The integrals Tp = gl{x,Q'^)dx and Vd = Jq gf{x, Q'^)dx 
were evaluated at fixed — 3 (GeV/c)^ using the full data set to yield Tp = 
0.132 ± 0.003(stat.) ± 0.009(syst.) and Ta = 0.047 ± 0.003 ± 0.006. 
The dependence of the ratio gi / Fi was studied and found to be small for 
> 1 (GeV/c)'^. Within experimental precision the g2 data are well-described 
by the twist-2 contribution, g^^ . Twist-3 matrix elements were extracted and 
compared to theoretical predictions. The asymmetry A2 was measured and found to 
be significantly smaller than the positivity limit \fR for both proton and deuteron 
targets. is found to be positive and inconsistent with zero. Measurements of 
gi in the resonance region show strong variations with X and , consistent with 
resonant amplitudes extracted from unpolarized data. These data allow us to study 
the dependence of the integrals Tp and T^ below the scaling region. 

PACS Numbers: 13.60.Hb, 29.25.Ks, ll.50.Li, 13.88.+e 
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I. INTRODUCTION 



Inelastic lepton scattering from nucleons has been used over the past thirty years to 
obtain an ever-increasing knowledge of the distribution of the partons that make up the 
nucleon, namely gluons and up, down, strange, and perhaps charmed quarks. It is one 
of the great successes of QCD that the same parton densities can be used to describe the 
unpolarized inelastic structure functions Q^) and F2(x, Q^) as well as many other 

physical processes, such as the production of jets in pp collisions. The parton densities 
depend on the fractional momentum of the parton x = Q'^/2Mv, where —Q^ is the 
four-momentum transfer squared, M is the nucleon mass, and v is the lepton energy transfer. 
The measured dependence at fixed x of Q^) and ^2(0;, Q^) has been shown to be 
in very good agreement with the QCD-based evolution equations 

The Fi and F2 structure functions are sensitive to the helicity-averaged parton densities. 
Recent improvements in polarized lepton beams and targets have made it possible to make 
increasingly accurate measurements of two additional structure functions, gi^x.Q"^) and 
g2{x,Q'^), which depend on the difference in parton densities with helicity either aligned or 
anti-aligned with the spin of the nucleon. Measurements of Qi have been made using electron 
beams at SLAC 0J|,|,^ and muon beams at CERN ||^,|^ , while (7" has been measured both 
using polarized deuteron targets at SLAC ||,^ and CERN 0, and a ^He target at SLAC 
TO| , p!I[| and DESY [jl2[. Measurements have also been made of g2 for both the proton 



and deuteron |T0|JT3| , p!^JT5|] , although with limited statistical precision compared to the gi 
measurements. This paper reports final results for g^, gf, g2, and g2 from experiment E143 
at SLAC, and includes more details of the analysis procedure, as well as some auxiliary 
results not covered in the original short publications P,|5|,|8| JT3|JT6[| . 

The earliest experiments sparked considerable interest in the spin structure 

functions when it was reported that, contrary to the quark model expectation, the quarks 
contribute very little to the proton's spin (the so-called "spin crisis"). Subsequent precision 
measurements are consistent with the original experimental results (with improved QCD 
corrections applied), but the theoretical interpretation has become more complex. It is 
now believed that in addition to the quarks, the orbital angular momentum and gluons may 
contribute significantly to the proton's spin. There is still the unanswered question as to how 
much the gluons alone really contribute. The gi and g2 structure functions are interesting 
not only in opening a new degree of freedom with which to explore the detailed structure of 
the nucleon, but also for making a precise test of QCD via the Bjorken sum rule which is a 
strict QCD prediction [p!7[] . 

In this paper we describe the theory and phenomenology of spin structure physics, 
and detail the SLAC experiment E143, which measured both Ay and A± for proton and 
deuteron targets over a wide range of kinematics. The theory and experimental apparatus 
are described in Sections II and III. The analysis procedure is detailed in Section IV. Results, 
their interpretation, and a discussion of systematic errors are shown in Section V, and finally 
we present a summary and conclusions in Section VI. 
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II. INTERPRETATION AND THEORY 



A. Formalism 



The structure functions gi{x,Q'^) and g2{x,Q'^) are typically extracted from asymmetry 
measurements. Longitudinally polarized leptons are scattered from a target that is polarized 
either longitudinally or transversely. The longitudinal (Ay) and transverse {A±) asymmetries 
are formed from combining data taken with opposite beam helicity: 



The polarized structure functions can be determined from these asymmetries: 

1 



9i{x, Q') = [All + tan(e/2)Aj 



92{x,Q ) = 



2d' 



E + E'cos(e) , 



E'sm{e) 



(2) 



where E is the incident electron energy, E' is the scattered electron energy, 6 is the 
scattering angle, y = {E - E')/E, d' = [(1 - e)(2 - y)]/[y{l + eR{x,Q^))l e^^ = 
1 + 2[1 + 7-2] tan2(6'/2), 7 = 2Mx/y/Q^, M is the nucleon mass and R{x, Q^) = OL/ar is 
the ratio of longitudinal and transverse virtual photon-absorption cross sections. Q^) 
is related to the spin-averaged (or unpolarized) structure functions Q^) and F2(x, Q^) 
by R{x, Q2) + 1 = (1 + 72)F2(x, Q^)/[2xF,{x, Q% 

The virtual photon-absorption asymmetries Ai and A2 are related to the measured 
asymmetries by 

A\\=D{A^ + riA2), 

A^ = d{A2-CAi), (3) 
where the photon depolarization factor D = {1 — E'e/ E)/{1 + eR), rj = e\f(^ l{E — -E"e), 



d = D\l2e/{1 + e), and ^ = 77(1 + e)/2e. Ai and can be expressed as: 

v4|| iqA^ 9i{x,Q^) -i^g2{x,Q^) 



Ai 



T T 
^1/2 ~ ^3/2 


(JxT 


^1/2 + ^3/2 




2a LT 




^1/2 + ^3/2 





D(l + r/C) d{l + r^C) F^{x,Q 



2-] 



^ = r±L = C^ll , ^ j{gi{x,Q^)+g2{x,Q^)) 

' • - D{l + r^C) d{l + r,C) F,{x,Q') ' 

where ^^"^ '^3/2 virtual photoabsorption transverse cross sections for total helicity 

between photon and nucleon of 1/2 and 3/2 respectively, a^T is the interference term 
between the transverse and longitudinal photon-nucleon amplitudes, ctt = + '^3/2)/'^^ 
and cttt = (o"i/2 '^3/2) f^- from Eq. ^ that for low x or high (where 7 << 1), 



^1 ~ Qi/Fi. Positivity constrains |y4i| < 1 and \A2\ < y R{x,Q'^). For the case where 
only the longitudinal asymmetry is measured, and a model is used for g2, Ai and gi can be 
expressed as 
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^11 



1 + 



xMy 



9i 



d' [' ' E + E'cos{e) 
A||Fir E + E' 
d' [e + E'cos{6) 



92_ 



4xME'cos2(|; 



+ 92 



p{E + E'cos{e)) 
2Mx 



E + E'cos{e) 



(5) 



In the resonance region, gi and g2 are well defined but are more properly interpreted in 
terms of the helicity structure of the resonance transition amplitudes. The 'j*NN* vertex 
for electro-excitation of the resonance A^* is generally given in terms of three amplitudes, 
^1/2 (Q^), ^3/2 (<5^) and 5*1/2 (Q^) Here, A denotes transverse photon polarization 

and S indicates longitudinal photons. The index 1/2 or 3/2 refers again to the total 7*A^ 
helicity. The virtual photon-nucleon cross sections for an isolated resonance can then be 
written in terms of helicity amplitudes as 



a 



1/2 



^3/2 



C^l/2 



0^1/2 



(i^i + 91 



2Mx 



92) 



KM ■ - V 



M.. . 



C^LT 



K 



^(1 + -)-^ 



(^1 + 92) 



2n—b—\Si 
W q 2 



W q* 2 2 



(6) 



K Mv 

in which K is the incoming photon flux which is chosen using the Hand convention such that 
the invariant mass squared of the final state is W"^ = + 2MK, b is the resonance line 
shape (unit area), and g*^ = Q"^ + {W^ — — Q^)/AW'^ is the squared magnitude of the 
3-momentum transfer measured in the resonance rest frame. The electron scattering cross 
sections are then written 



where 





da 




dE'dVL 


da^^ 




dE'dn 


dE'dVL 


da^^ 


da^^ 


dE'dn 


dE'dVL 







Tv[(yT + ec^i], 

2VyD{l + eR)[aTT + riaLTl 
2Vvd{l + tR)[aLT-C<yTTl 

a K E' 2 



(7) 



An^Q^ E 1-e 



(8) 



B. The Deep-Inelastic Spin Structure Function gi{x,Q'^) 

As will be shown below, the first moment of the spin structure function gi{x, Q"^) is related 
to the net quark helicity AS which contributes to the proton spin. Angular momentum 
conservation requires that 

l = ^ + AG + L,, (9) 
where AG is the net gluon helicity, and is the orbital angular momentum. 
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1. The Quark-Parton Model 



In the naive quark-parton model (QPM) the nucleon is composed of quarks which have 
no orbital angular momentum, and there are no polarized gluons present. In this simple 
picture, the unpolarized structure function Q^) and the polarized structure function 
gi{x,Q^) can be simply expressed as the charge- weighted sum and difference between 
momentum distributions for quark helicities aligned parallel (g^) and antiparallel (g^) to 
the longitudinally polarized nucleon: 

^ i 

9i{x) = \Y.^l[ql{x) - q\{x)] = E^'^g^x) . (10) 

i i 

The charge of quark flavor u,d, and s is denoted by Cj, and (x) are the quark plus 
antiquark momentum distributions. The quantity Jq Aqi{x)dx = Aqi refers to the helicity 
of quark flavor i in the proton, and AS = An + Ad + As is the net helicity of quarks. Since 
AG = and = 0, it follows from Eq. ^, that AS is expected to be unity in this model. In 
a relativistic quark-parton model pO| , |2T|j2^ (with no polarized gluons), the orbital angular 
momentum contribution is no longer zero and the quark helicity contributions to the proton 
helicity are suppressed by a factor of about 0.75. 



2. Perturhative QCD and the Role of the Gluons 



The quark-parton model is useful for understanding some properties of the nucleon such 
as charge and isospin. However, it fails to adequately describe all properties, and it falls short 
in explaining the dynamics of particle interactions. For this we need a more comprehensive 
theory such as Quantum Chromodynamics (QCD) which can account for gluons and their 
interactions with the quarks. 

The operator product expansion (OPE) p3tp4| , p5| is a useful technique within QCD 
because it separates the physics into a perturbative part that is easily treatable and 
a non-perturbative part that is parameterized in terms of unknown matrix elements of 
Lorentz-covariant operators. At leading twist the first moment of gi{x, Q"^) can be expressed 
in terms of singlet (oq) and nonsinglet {a^ and ag) proton matrix elements of the axial current: 



/■I , /I 1 \ 1 

ri(Q^) = gl{x,Q'^)dx = ag + —asjCns + -aoCs, 



rim = j\t{x,Q')dx =[l- 



1 ^ 1 ^ 
6b 9 



(11) 



where ujd is the D-state probability in the deuteron, and the factors Cns and Cg are the 
Q^-dependent non-singlet and singlet QCD corrections, which are discussed in more detail 
below. 

Assuming that there are no polarized gluons contributing to the proton spin, the singlet 
and nonsinglet proton matrix elements given in Eq. |ll] can be related to the quark helicities: 
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ao = Au + Ad + As = AS, 

as = Au - Ad = F + D, 

as = Au + Ad- 2As = 3F -D. 



(12) 



Here, F and D are weak hyperon decay constants which can be extracted from data assuming 
SU{3) symmetry [g6],^ 



F + D = gA= 1.2601 ± 0.0025, 
= 0.588 ±0.033. 



(13) 



The error quoted above on 3F — D is the experimental error assuming SU(3) symmetry. 
It may be an underestimate because possible SU{3) symmetry breaking effects could 
be significant. There have been a number of attempts to estimate these effects 

symmetry breaking effects in the past have 



30|gT|,|32|. According to Ratcliffe 



always been found to be at most 10%. Assuming a generous 20% systematic error from 
symmetry breaking combined with the above error in quadrature yields an error of 0.12 
on 3F — D. This error is somewhat smaller than the range of possible values (0.40-0.84) 



presented under various assumptions p8| , |29| , |31| , |32| , some of which have come under criticism 



After combining Eq. |TT] and Eq. |T2| ||2^ it is straightforward to extract the singlet matrix 
element from the measured first moments of the proton: 



ao 



and the deuteron: 



ao 
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9 



l-|^D 



1 



(3F + D)C„ 



(14) 



36 



(3F - D)C„ 



(15) 



The nonsinglet QCD correction Cns ||33| calculated in the MS scheme to order three for three 
quark flavors is given by 



- 20.22 



TT 



TT 




(16) 



where as{Q'^) is the strong coupling constant. Fourth order QCD corrections have been 
estimated ||3^ to be small at the kinematics of this experiment. The singlet QCD correction 
exists in two forms I^S], one which yields a Q^-dependent ao(Q^) in Eqs. p^p!5| and one which 
yields Oq"'' which is the asymptotic high limit of ao(Q^). These singlet QCD corrections 



have been calculated in the MS scheme []33| : 



( asm 



TT 



0.3333 



7r 



0.5495 



TT 



TT 



(17) 
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The contribution of AG to (and thus to the first moment of gi) is a factorization 
scheme-dependent quantity. In a gauge-invariant scheme such as the MS scheme gluons do 
not contribute to the first moment of gi which means that Oq = AE. For chiral-invariant 
schemes such as the Adler-Bardeen scheme |^ the gluons do contribute to Fi. The physical 
quantity gi{x) is independent of the factorization scheme, however. In the Adler-Bardeen 
scheme 



35[| , the quantity ao 7^ AS, and is instead written as 



AE 



(18) 



This contribution of AG is called the gluon axial anomaly |36| or the Adler-Bell-Jackiw 
anomaly (as applied to QCD from QED). The product as{Q^)AG{Q^) is independent of 
in leading order which implies that AG(Q^) grows in like l/asiQ"^), -^2 compensates 
to satisfy Eq. ^ Physically, this means that as each quark radiates a gluon with some 
preferential helicity the orbital angular momentum of the quark-gluon system must increase 
to conserve the total angular momentum. Thus, as more gluons are emitted, both AG and 
Lz will grow, but with opposite signs. 

Other quantities of interest are the helicity contributions from the individual quarks. 
These quantities can be extracted from the measured oq, but may be subject to possible 
gluon contributions as in Eq. |18|. Allowing for the possibility of gluon contributions, these 



quark helicities are calculated using: 



Au = i(ao + 3F + D) + ^asiQ^)AG{Q^) 
6 Ztt 



Ad = -(ao - 2D) 



27T 



asiQ^)AG{Q' 



(19) 



As = -(ao -3F + D) + —as{Q^)AGiQ^) 
o Ztt 



If we include a contribution of AG(1 (GeV/c)^ 
0.118 ± 0.003 [H], we calculate a,AG/27r 



-- 1.6 ± 0.9 0, and use as{Mz) = 
0.13 ± 0.08 and find good agreement with 
existing data. This model along with quark-parton model expectations are summarized in 
Table | and can be compared with data from this experiment in Table pCXIX| . Note that the 
value used above for AG agrees well with a theoretical prediction based on QCD sum rules 
p8| which yields AG(1 (GeV/c)^) = 2.1 ± 1.0 and on an earlier parameterization |^ which 
yields AG(1 (GeV/cf) = 1.7. 

There are a number of other theoretical models which attempt to explain how the quark 
helicity is distributed within the nucleon. Non-perturbative effects enhancing the role of 



intrinsic sea quarks have been proposed by several authors. Halperin and Zhitnitsky EO 



argue that a large portion of the nucleon spin comes from charm quarks by adding a term 
2Ac to the ao term in Eq. |11|. Brodsky and Ma |^1[| contend that asymmetries in the 
light quark sea could generate the observed AS. The Skyrme model predicts that 
AS = Au + Arf + As=0 and AG = and = 1/2, and should be accurate to 0(1/A^c) 
where Nc = 3, the number of colors. Within its uncertainty this is consistent with the small 
observed value of AS. Other models include the chiral bag model the chiral quark 
model calculations based on QCD spectral sum rules or Pauli-exclusion principles 
45| , and also lattice QCD predictions [^| 



8^49 
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3. The Bjorken Sum Rule 



This sum rule was originated by Bjorken using current algebra and isospin symmetry. 
It has since been re-derived in QCD and is a strict prediction made by this theory. It relates 
the integral over all x at fixed of the difference between gi{x,Q'^) and gi{x,Q^) to the 
well-measured neutron beta decay coupling constant = 1.2601 ± 0.0025 ||26| , 

nm - r-,m = J {g'li^, q') - 9ii^, q')) dx = c^s. m 

An experimental test of this sum rule provides a test of fundamental QCD assumptions. In 
addition, it is possible to use the measurement to extract a relatively accurate determination 
of asiQ"^) at low (on the order of 2 to 10 (GeV/c)^) [^]. A significant difference from 
other as{Q^) determinations could indicate the presence of interesting new physics. 



4- The Ellis- Jaffe Sum Rule 



The other sum rules of interest for gi, although less rigorous than the Bjorken sum rule, 
are the Ellis- Jaffe sum rules ||51|] which were derived using SU(3) symmetry and assuming 
the strange sea in the nucleons is unpolarized. These sum rules, including the necessary 
QCD corrections, follow naturally from Eqs. |TT] and with As = AG = such that 
ao = as = 3F — D: 

r?(g')= l\l{x,Q^)dx = ^[Cns{:iF + D) + 2Cs{?>F-D)] , 
Jo lo 

r?(Q') = /' 9i{^. Q')dx = I l-DCns + CsiSF - D)] . (21) 
JO 9 



5. Dependence: Evolution and Higher Twist 



The quark-parton model does not inherently include gluons, and it is the interaction 
between the quarks and gluons which generates the observed dependence of both the 
polarized and unpolarized nucleon structure functions. The QCD theory which describes 
the quark-gluon dynamics gives predictions about how the parton distribution functions 
(and thus structure functions) evolve in in the perturbative limit of small a^. The 
evolution of the polarized parton densities is governed by the DGLAP equations which 
embody the emission of gluons by quarks. This gluon emission is responsible for the leading 
logarithmic dependence. In addition, there are higher-twist contributions to the 
dependence which are suppressed by powers of 1/ ^/Cp. Higher twist corrections to gi have 
been estimated to be small p2| , p3| , p4| , p5| while higher twist corrections for g2 have been 

Fits to Au{x,Q'^), Ad{x,Q'^), As{x,Q'^), and 



49.52.56.57 



estimated to be significant 
AG{x, Q^) have been made using next-to-leading-order (NLO) DGLAP equations 

||60| . The results indicate that NLO fits are more sensitive to the strength of the polarized 
gluon distribution function AG{x, Q"^) than leading order fits. 
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C. The Deep-Inelastic Spin Structure Function g2{x,Q'^) 



1. Physical Interpretation 



The interpretation of g2 in the naive parton model is less straightforward than that of gi. 
Feynman related the quantity gT{,x) = gi{x) + (72(2;) = ^2-^1/7 to the distribution of quark 
polarizations aligned parallel (fc^) and antiparallel {k^) to that of a transversely polarized 
proton p4| , |6l| by the expression 



9t{x) 



2xM 



x) 



(22) 



where rriq is the quark mass. Leader and Anselmino ||6^ subsequently derived the parton 
model expressions for gi and g2 for a nucleon polarized at an arbitrary angle 6 relative to 
the incident electron direction. Evaluated at ^ = their expression for gxix) is 



9t[,x) 



2xM 



[q]{x) - qlix)], 



(23) 



where qi are the same as in Eq. [1^. Jaffe and Ji [2^| pointed out that claims that gr is 

small were generated by setting = in Eq. For consistency, if the quark momenta are 
taken to be along the longitudinal direction then = xM [|^,^, which yields g2{x) = 0. 
Because of Fermi motion, however, the quarks are off-shell and rriq 7^ xM in general. It is 
this large off-shell nature of the quark which produces large twist-3 effects in the MIT bag 
model calculation |2i . 

The naive parton model does not include transverse momentum or quark-gluon 
interactions. For this we can turn to a more advanced light-cone parton model, |]65|j66[] 
or an OPE analysis, p4|,i3 "which indicates that there are three components (up to twist-3) 
contributing to g2- These components include the leading twist-2 part gY^{x,Q'^) p8| , 
coming from the same set of operators that contribute to gi, another twist-2 part coming 
from the quark transverse-polarization distribution hj^ix^Q'^), and a twist-3 part coming 
from quark-gluon interactions ^(x, Q^): 



^^|(>(.,g^)+e(.,g% ^ 



(24) 



The term containing hriy, Q"^) is usually neglected because it is suppressed by the quark 
mass rrig, and the g^'^ expression of Wandzura-Wilczek is given by 



-gi{x,Q^' 



9i{y,Q^ 



dy 



(25) 



2. OPE Sum Rules and the Twist-3 Matrix Element 



Keeping terms up to twist-3, the OPE analysis of gi and g2 yields an infinite number of 
sum rules: 
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in) 



p(n) 
J- 2 
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(26) 



2 n + l 

where are the twist-2 and dn are the twist-3 matrix elements of the renormahzed operators. 
The OPE only gives information on the odd moments of the spin structure functions. Note 
that contributions involving ntq/M (see Eq. have been left out of Eq. ^ as have target 
mass effects discussed below. The twist-3 matrix elements follow from Eq. |2B: 
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(27) 



where g^ = g2 — g2^^- We see from Eq. ^ that if all dn 



then g2 is completely determined 
by gi because there are an infinite number of sum rules. This is how the quantity g^^ was 
originally derived. 

There are a number of theoretical predictions for d2 for both proton and neutron targets 
49| , |52|j56| , |57|j58| . Some are based on bag models |l56| , |57| , others on QCD sum rules ||52| , |58 



and there is also a lattice QCD calculation Many of these models have predicted large 
values for d2 which means there could be significant twist-3 contributions to g2- This makes 
the study of g2 particularly interesting. 



3. The Burkhardt-Cottingham Sum Rule 
The Burkhardt-Cottingham sum rule for g2 at large Q^, namely 

J\2ix)dx = 0, (28) 

was derived from virtual Compton scattering dispersion relations. This sum rule does not 
follow from the OPE since the n = sum rule is not defined for g2 in Eq. ^ The validity of 
the Burkhardt-Cottingham sum rule relies on g2 obeying Regge theory at low x, which may 
not be a good assumption. A non- Regge divergence of g2 at low x would invalidate this sum 
rule although such a divergence could be very difficult to detect experimentally. 



4- Target Mass Effects 



The OPE sum rules as given in Eq. ^ were derived in the limit M^x^/Q^ —>■ 0. These 
target mass effects can become significant when M'^/Q'^ is of order unity which is certainly 
the case for a subset of the data presented in this paper. These target mass effects for 
polarized electroproduction have been determined [^,^]. The corrected Bjorken sum rule 
derived from these formulae is given by 



dx 



e 



5 + 4JI + 4M2xVQ2 



gl{x,Q')-g",{x,Q' 



1 M^e 





gP,ix,Q')-g^ix,Q' 



77 gA Cms, 



(29) 
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where the Nachtmann variable ^ = 2x/(l + ^1 + AM'^x'^/Q'^). This sum rule is now 
dependent on g2. The size of the target mass effects |]72| , |73| to the uncorrected Bjorken 
sum rule formula are estimated to be of the same magnitude as higher-twist effects which 
are typically small. The target mass effects for the calculation (missing in Eq. ^) have 
been investigated and are negligible for our kinematics ||7i| . 



D. Resonance Region Polarized Structure Functions 



The values of Ai in the resonance region are a combination of the asymmetries for 
individual resonances and for the nonresonant background. Resonance helicity amplitudes 
Ai/2 and are reasonably well measured at = for the prominent resonances p6[ 



Sparse data exist also for virtual photons |j75|. The excitation of the A (1232) resonance 
(spin-|) includes both | and | spin projections. At low Q'^ the A(1232) excitation is 
expected to be primarily a magnetic dipole transition for which ^3/2/^1/2 = and 
Ai = (IA1/2P - |A3/2n/(|Ai/2p + |A3/2n = -i. For real photons A3/2M1/2 = 1.064v^ 
p6| . Perturbative QCD predicts that the ratio ^3/2/^1/2 should go as l/Q^ and Ai should 
approach unity as ^ 00. However, a recent analysis of pion electroproduction data 
shows that the magnetic dipole transition still dominates at = 3.2 (GeV/c)^. On the 
other hand, the 5*11(1535) resonance has no spin-| projection, so Ai should be unity at all 
Q'^. Data |75[ from Bonn, Daresbury and DESY have been used to extract y4i/2 and ^3/2 up 
to = 3 (GeV/c)^ for the 6*11, D13, and F15 resonances. Because of the large uncertainties 
of these extractions, our knowledge of the dependence of the helicity amplitudes is still 
rudimentary. The asymmetries Ai for both D13 and F15 make a transition from Ai ^ — 1 to 
y4i ^ 1 somewhere in the range < < 3 (GeV/c)^ 



Less is known about a^T and A2. The positivity limit A2 < ^R{x,Q'^) constrains A2. 
The world average value for R in the resonance region is 0.06 ± 0.02 for 1 < < 8 



(GeV/c)^ and W"^ < 3 GeV^. Since this is smaller by half than the deep-inelastic fit to 
R{x, Q"^) 1^ extrapolated into the resonance region, one might argue that -R(x, Q"^) for the 
resonances themselves is small. However, little is known for < 1.3 (GeV/c)^. 

A complete mapping of gi{x,Q'^) at low where the resonances dominate is useful 
for two reasons. First, these data provide important input for radiative corrections 
of the deep-inelastic data. Second, the evolution of the integral (defined to exclude 
elastic scattering) Vi{Q'^) = Jq gi{x,Q'^)dx for should be determined by the 

Gerasimov-Drell-Hearn (GDH) sum rule ||7^ for real photons: 

in which u is the photon energy, is the threshold energy for pion production, k is the 
nucleon anomalous magnetic moment, and M is the nucleon mass. A simple change of 
variables from z/ to x in Eq. [3^ and a reformulation of aTT in terms of gi and g2 yields 

r ri(g^) 
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One important feature of Eq. |3T] is the sign. for the proton is positive for above 
3 (GeV/c)^ as measured in the deep-inelastic regime. However, the GDH sum rule predicts 
that should become negative at small Q^. This implies that somewhere in the region 
< < 3 (GeV/c)^, V\ must cross zero. Exactly where this occurs depends crucially on 
the evolution of the resonance helicity amplitudes which are presently not well known. 
Predictions about how ri(Q^) goes from the deep-inelastic values to the GDH limit have 
taken one of two paths: 1) theoretically motivated interpolation and 2) computations that 
include all available knowledge of the resonance behavior. Both need to be checked with 
direct measurements of ri((5^) at low Q^. 

Recent theoretical work |80| indicates that at low Q^, ri((5^) corrected to include the 
elastic contribution, 

rf (g^) = iFi(Q2)[Fi(Q2) + F^{(f)\ (32) 

provides the twist-4 (l/Q^) corrections to the Ellis- Jaffe sum rule. Fx and F2 in Eq. p2| are 
the Dirac and Pauli elastic form factors (not the deep- inelastic structure functions). 



III. THE EXPERIMENT 



The goal of the E143 experiment was to determine the longitudinal and transverse 
cross-section asymmetries via deep-inelastic scattering of longitudinally polarized electrons 
from polarized protons and deuterons. Over a period of three calendar months data were 
taken at beam energies of 29.13 GeV (122 million events), 16.18 GeV (56 million events) 
and 9.71 GeV (58 million events). The longitudinal asymmetry A\\ was obtained with the 
target polarization parallel to the beam momentum, whereas the transverse asymmetry 
(at F = 29.1 GeV only) was obtained with the target polarization transverse to the beam 
momentum (right or left of the beam). 

The experimental apparatus employed consisted of five components: the polarized source, 
the accelerator and beam transport, the M0ller polarimeter to measure the beam polarization 
in the End Station A (ESA), the polarized proton/deuteron target, and the two spectrometer 
arms to detect the scattered electrons. These components are discussed in the sections which 
follow. 



A. The SLAC Polarized Electron Source 

A polarized electron source for the SLAC Linear Accelerator was first developed in the 
early 1970's for experiments on the spin structure of the proton. Since 1978 the SLAC 
polarized electron source has been based on the principle of laser photoemission from a 
gallium arsenide (GaAs) photocathode. Strained GaAs photocathodes, which effectively 
doubled the polarization obtainable from an unstrained GaAs photocathode, were developed 
in 1991 |81|,8^ and first used in a SLAC experiment in 1993. The design and operational 



characteristics of the SLAC polarized electron source are fully described in Ref. fS^ . 



Figure. |I| is a schematic of the layout of the laser and gun structure at the SLAC injector. 
Electrons are photo-emitted from a GaAs photocathode by illuminating the surface with 
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a laser. The electrons are polarized with a helicity defined by the sign of the circular 
polarization of the incident laser light. Spin reversals are achieved at the source by reversing 
the circular polarization of the laser light with a Pockels cell. The pattern for the sign of 
the polarization is chosen to be a known pseudo-random sequence, permitting validation of 
the sign of each pulse in the offline data stream. 

The polarization of the electrons is a consequence of the band structure of GaAs and 
the angular momentum selection rules that apply to this system. The presence of strain 
changes the lattice constant of the GaAs, shifts the energy levels, and breaks a spin state 
degeneracy in the valence band. Excitation by circularly polarized light near the band gap 
edge (A 850 nm) will then result in only one set of spin states populating the conduction 
band provided the strain is sufficiently large. The strain in the active GaAs layer is achieved 
by growing a 100 nm thick epitaxial layer of GaAs on a substrate of GaAs(i_x)Px (x=0.28). 
The lattice spacing for the GaAsP is about 1% smaller than for GaAs, and the resulting 
lattice mismatch puts the GaAs epitaxial layer under a compressive strain sufficiently high 
to remove the spin state degeneracy. Under these conditions one expects that photo-emitted 
electrons will have close to 100% polarization. In practice, the electron polarization averaged 
85% for the E143 experiment. 

The laser system was designed and built at SLAG. It consisted of a flash- lamp pumped 
titanium sapphire rod, producing light pulses which were optically chopped to a 2.3 /xsec 
long pulse. The laser beam was transmitted through a lens system which allowed for steering 
and focusing on the cathode. For the E143 experiment, the amount of laser power available 
was larger than needed, so the power was attenuated to about 10 watts peak, yielding 
approximately 4 x 10^ electrons per pulse. At this low electron intensity, the accelerator 
control system was unable to sense the presence of beam. To allow the accelerator controls 
to operate, the beam was intentionally intensified to about 2 x 10^*^ in one of the 120 pulses 
generated per second. This "witness pulse" was then sent into a beam dump before reaching 
the target, and the experiment operated on the remaining 119 pulses per second. 

Possible systematic errors associated with reversal of the electron spin are important to 
this type of experiment. Correlations between beam current, beam energy, beam positions, 
and beam angles on the target were available to the experimenters on a short time basis 
from beam monitors. For this laser-driven photoemission source, the reversal of the laser 
polarization is sufficiently free of unwanted effects such that all systematic errors from the 
source were negligible. 



B. The Electron Beam 

1. Production and Transport 

The electrons produced by the polarized source were accelerated to energies between 9 
and 30 GeV in the linear accelerator. The electrons were then deflected through an angle of 
24.5° in the A-line beam transport and were directed onto the polarized target in the ESA. 

Because of its anomalous magnetic moment, the spin of the electron precesses by an 
angle larger than that of the bend angle of the beam, according to the formula: 
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where g is the gyromagnetic ratio, E is the energy (in GeV), m is the mass of the electron, 
and A0 is the angle between the electron spin and the momentum at the target. When A0 is 
an integral multiple of tt, the electron spin is longitudinal at the target. The experiment was 
run at energies of 9.71, 16.18 and 29.13 GeV, corresponding to Stt, Stt and 97r respectively. 
By varying the energy around the nominal value and measuring the longitudinal polarization 
in the M0ller polarimeter, we verified that the chosen energy produced the maximum 
polarization. 

2. Beam Monitoring 

The incident flux of electrons was measured independently in two identical precision 
toroidal charge monitors in the ESA. These were frequently calibrated with a known charge 
and agreed to better than 1%. The response of the toroids is independent of the polarity of 
the beam. 

The position of the beam at the target was monitored in two devices: a traveling-wave 
radio-frequency beam position monitor which was non-interfering and was placed just in 
front of the target, and a pair of secondary emission foil arrays with 1 mm spacing located 
10.8 m downstream from the target. The former provided a direct measurement of beam 
centroid position, and was used in an automatic feedback system to keep the beam on target; 
the latter allowed a measurement of both the position and the transverse dimensions of the 
beam by comparing the charge collected on the individual foils. 

3. Beam Rastering 

To minimize effects such as target depolarization from local beam heating and radiation 
damage, the beam was moved or "rastered" across the face of the target. The beam position 
was changed between pulses by means of a pair of air-core magnet coils located 67 m 
upstream of the target. The lack of iron in these magnets allowed the fields to be quickly 
changed under computer control. The beam at the target was rastered on a grid with a 
spacing of 1.2 mm inside a circle of radius 10.8 mm for a total of 253 points. Thus, each 
point in the target was illuminated only once every 2.1 seconds. The raster pattern skipped 
every other point and row, so that subsequent pulses did not overlap, and the entire raster 
pattern was completed in four passes. Because the dimensions of the beam (Gaussian a of 
2 mm horizontally and 1 mm vertically) were comparable to the raster spacing, the overall 
illumination of the target was quite uniform inside the circle of the raster. 

4- The Beam Chicane 

For the measurement of g2, the target was rotated by 90° so that the target nucleons were 
polarized transversely to the beam direction in the scattering plane. In this configuration, 
the electrons in the beam passed through / Bdl = 1.52 T-ni as they traversed the target. 
This was enough to deflect them through an angle of 0.90° at 29.1 GeV and to rotate the 
polarization vector through an angle of 60°. As a result, the beam after the target would 
have no longer been parallel to the nominal beam-line, and in fact would have been about 30 
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cm low at the exit of the ESA. More significantly, the deep-inelastic scattering would have 
taken place at a different average angle and longitudinal polarization than in the parallel 
case. 

To compensate for the effects of this magnetic field, we inserted four identical dipole 
magnets (the chicane) into the beam-line, three upstream and one downstream of the target. 
The first magnet deflected the beam down by 0.45°, and the second pair bent the beam back 
up by twice this amount. This caused the beam to arrive at the center of the target with 
both the momentum and polarization vectors horizontal. After exiting the target, the beam 
was tilted downward, and the fourth magnet returned the beam to the horizontal so that it 
left the ESA parallel to the nominal beam-line displaced vertically by only 3.5 cm at 29.1 
GeV. 



C. Beam Polarimetry 

A M0ller polarimeter was used to measure the beam polarization during the E143 
experiment. This is a practical and reliable approach based one + e^e + e scattering, a 
spin-dependent QED process with large a cross-section and analyzing power. The expected 
cross-section asymmetry can be calculated with high precision P3 and is not significantly 
modified by radiative processes 

For a beam with longitudinal polarization Pb and target with longitudinal polarization 
Pt, the beam polarization is measured by comparing the relative cross-section asymmetry 
for beam and target spins aligned parallel (H) and anti-parallel (ti): 

The relative cross-sections are determined by detecting either of the scattered electrons or 
both in coincidence. 



1. Layout 

A schematic of the polarimeter is shown in Fig. ^ The major components are a polarized 
electron target, an acceptance-defining collimator, a dipole magnet spectrometer, and two 
independent detector systems. One system detected the M0ller electrons in coincidence, 
whereas the other integrated the single electrons over the duration of the beam spill. 

The polarized electron target consisted of six magnetized ferromagnetic foils of different 
thickness which could be moved into the beam. The foils were magnetized to near saturation 
by Helmholtz coils providing nearly 100 gauss at the target center. The permendur (49% Fe, 
49% Co, 2% V) foils were 3 cm wide and varied in thickness from 20 fim to 150 ^m. 
The target electron polarization (typically 0.082) was determined to a relative accuracy of 



1.7% from foil magnetization measurements. 

The tungsten collimator which was 20 radiation lengths thick (see Fig. |^) had a central 
opening to allow the main beam to pass and wedge shaped apertures of constant azimuthal 
acceptance (0.2 radian top, 0.22 radian bottom) above and below the beam-line to select 
M0ller electrons scattered transverse to the bend plane of the downstream dipole magnet. 
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The vertical acceptance was 3.6-9 mrad in the lab (corresponding to 70-116°. in the 
center-of-mass) . The 2.1 T-m dipole field separated the scattered electrons according to 
momenta. Since M0ller scattering is elastic, the x and y position of the scattered electrons 
at the detector plane are correlated as shown in Fig. ^ 

The detector hut was situated 27 m downstream from the M0ller target. The single-arm 
detector package of three radiation lengths of lead and a single plane of position-sensitive 
silicon detectors was placed immediately in front of the coincidence detectors. The 
coincidence package consisted of two arrays of seven lead glass blocks (SF-6), each with 
a 10x10 cm^ entrance area and 25 cm of depth (~15 radiation lengths). 

2. The Measurements 

Measurements of the beam polarization were performed every one to two days. Each 
measurement period typically consisted of four runs using two target foils (thin and thick) 
and opposite target polarization directions. This made it possible to look for rate or helicity 
dependent effects in the data. For polarization measurements, the beam rastering was turned 
off and the beam focus was moved to the M0ller target. Otherwise, the beam conditions 
were identical to that of the main experiment. Data were obtained from almost 200 runs 
over a range of luminosities (more than a factor of 8) through different combinations of 
foil thickness and beam current. The coincidence polarimeter obtained a typical statistical 
precision of 0.010 (absolute) per run whereas the single arm polarimeter achieved 0.019 per 
run. Both detector systems took data at 29.1 and 16.2 GeV. The single arm collected data 
with the 9.7 GeV beam. 

3. The Coincidence Polarimeter 

The segmented lead glass arrays provided good energy and timing information and made 
it possible to accommodate the high instantaneous rates of several 10s of MHz characteristic 
of the low duty factor (10~^) at SLAG. The combination of Cerenkov light in the glass blocks, 
fast photomultiplier tubes, and a clipping circuit resulted in signal pulse widths as narrow 
as five nanoseconds. The signals were fed into an Ortec 935 constant fraction discriminator 
with the threshold set at 30-40% of typical M0ller signal amplitude. The discriminator 
output was then fed into a fast multiplexing circuit with a fanout to three or four TDG 
channels. The time of each event was recorded by a LeGroy 2277 multi-hit TDG which has 
a least significant bit time of one nanosecond. The multiplexor was required to decrease 
the dead-time and increase the maximum hit capability of the individual TDG channels. A 
laser pulser, triggered randomly <l/spill, was fed into each detector block simultaneously 
through fibers to provide both time calibration and detector dead-time information for the 
analysis. 

The data were recorded on tape on a spill-by-spill basis as a series of event times and 
corresponding TDG channels. The analysis identified coincidence events by the arrival of 
single event times within a predefined time window of ±4 ns, determined by the resolution of 
the TDG's. Frequently, the analysis encountered ambiguous coincidence pair combinations 
where a single event of one detector could be combined with events in two or more other 
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detectors. Such ambiguities arose as a result of cross-talk between adjacent detectors due 
to shower sharing or due to random coincidences, particularly at higher luminosities. In the 
case of ambiguities, the cluster of all possible coincidences was subjected to a decision making 
routine which selected the most probable combinations of events. In the case of cross-talk 
events the full weight of the single coincidence event was shared with the adjacent coincidence 
pairs. Background contributions consisted of random coincidences between M0ller or Mott 
type electron events. Their contribution, typically <1% was estimated from the product of 
the singles rates in each conjugate detector. 

The dead-time measurement was obtained using the laser pulser system which sent a 
known pulse to all 14 detectors simultaneously. The efficiency at which both detectors of a 
pair saw the pulser event yielded the live-time for that pair. It was also necessary to correct 
for the possibility of two M0ller events occurring in a given pair within the same coincidence 
window. In such a case the system is only capable of seeing one of the pairs, an inefficiency 
which would be unaccounted for in the dead-time correction. To correct for this effect an 
estimate of the number of M0ller coincidence events occurring during the pulser event was 
added to the known number of pulser events for each pair. 

A typical M0ller coincidence time difference spectrum is shown in Fig. |^a. Two views of 
a typical distribution of coincidence events in the two detector arms are given in Figs, ^d-c 
for a run at 29 GeV. True M0ller events were kinematically restricted to occur only in 11 (9 
for the 16 GeV data) of the possible (7x7) 49 pairs. Radiative effects did allow some true 
M0ller coincidences to occur on the low momentum side of the ridge seen in Fig. |^c, but 
these events were not considered in the total event yield due to poor signal-to-noise ratio 
and greater uncertainties in the analyzing powers. 

The beam polarization for each coincidence pair was determined from the asymmetry 
in the yield corrected for background, dead-time, charge asymmetry, the effective analyzing 
power of each pair, the target polarization and the target angle. The polarization for a run 
was taken from the weighted average over all the pairs. 

A Monte Carlo analysis |87] was used to determine the effective analyzing powers of 
each coincidence pair and to evaluate the sensitivity of the analyzing powers to possible 



systematic influences such as the atomic motion of the target electrons |p8[ and typical 



shifts in the beam position or focus. The analyzing powers were found to range from 0.776 
to 0.690 for the different pairs. The effect of the target electron motion was to increase the 
average analyzing power by <0.5%. Typical beam parameter shifts resulted in changes to 
the average analyzing power within ±0.6%. Since only one set of analyzing powers was used 
at each beam energy, the ±0.6% variation was included as a source of systematic error. The 
large acceptance of the detectors reduced the sensitivity of the analyzing powers to these 
systematic influences. 

The possibility of rate dependence was investigated in two studies. One study compared 
low and high luminosity runs taken during a run set where little variation in polarization 
was expected during the set. In this study the internal agreement between all measurements 
of a run set was very good, resulting in an average per degree of freedom (df) of 1.1 for 
all the run sets. Another study tested the effectiveness of the analysis routine in dealing 
with ambiguities in the data which were most prevalent at high luminosity. Data taken from 
subsequent spills of a low luminosity run were artificially superimposed to create a fictitious 
spill of high luminosity. After imposing the detector dead-time on the single events the 
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data were analyzed as a normal run, and the yields could be compared with the original 
luminosity analysis. In both studies it was possible to rule out a rate dependence at better 
than ±0.5%. 

Despite the excellent agreement of the polarization results within a run set, a large 
fluctuation in the average polarization values obtained from each coincidence pair was 
observed. Although the sources of these fluctuations likely cancel in the average, their 
origins are not clear. As a result a maximum error contribution of 1.3% was included. This 
contribution reduced the x^/ df of the pair-dependent polarization distribution to unity. This 
uncertainty was combined in quadrature with the uncertainties estimated for the analyzing 
powers and the limit on a possible rate dependence to obtain a total systematic uncertainty 
of 1.5% for the coincidence polarimeter measurement. 



4- The Single-Arm Polarimeter 

The single arm detectors had four silicon pad detectors above and below the beam height. 
A lead converter absorbed soft photon backgrounds and amplifled the M0ller signal. Each 
detector consisted of two 4 (x) by 6 (y) cm silicon devices approximately 300 /im thick. 
Each device was segmented into 7 pads (channels) 8.70 mm wide and 40 mm long. Only 
12 contiguous channels were instrumented in each detector. The detectors were tilted by 
— 10.5° (top) and +11.0° (bottom) to align the channels along the M0ller scattered electron 
stripe. Since each detector was formed from two silicon devices there was a 5.3 mm gap 
between channels seven and eight. 

The silicon detector channels were connected to 96 charge sensitive preampliflers which 
integrated over the entire 2300 ns beam pulse. The preamplifler outputs were brought into 
ADCs to measure the peak of the preamplifler signal and were recorded together with the 
sign of the beam polarization for each beam pulse. 

The M0ller analysis proceeded through two steps. The flrst-pass analysis calculated 
average pedestal subtracted pulse heights and errors for each channel from the pulse-by-pulse 
data. Separate averages were made for pulses tagged by right (R) and left (L) handed 
polarization bits. Correlations between channels were calculated and recorded. A very loose 
beam current requirement was made before including the pulse in the overall averages. A 
summary flle containing the ADC averages and errors as well as useful beam and polarimeter 
parameters was written for each run. A second-pass analysis read the summary flle, applied 
channel by channel gain corrections, and formed sum (R-fL) and difference (R-L) averages 
and errors for each channel. Typical (R+L) and (R-L) line-shapes are shown in Fig. ^ for 
data at 29.1 GeV. 

The background under the unpolarized (R+L) M0ller scatters was estimated by fltting 
the (R+L) line-shape to an arbitrary quadratic background plus the line-shape expected from 
unpolarized M0ller scattering. The technique for estimating the unpolarized line-shape used 



the observed (R-L) line-shape and angular smearing functions to generate a predicted 
(R+L) line-shape for M0ller scatters. The observed (R+L) distribution was then flt by this 
predicted line-shape and a quadratic background. Since the observed (R-L) line-shape is 
already broadened by multiple scattering in the target material, beam windows, air, and 



helium, only corrections to the line-shape which are different ||8^ for scatters from polarized 



and unpolarized target electrons are included in the smearing function. 
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An analyzing power for each detector was calculated from the target polarization and 
the expected M0ller asymmetry determined by Monte Carlo simulations of the scattering 
process and detector response. The effect of the target electron momentum distribution 
88| was to modify the expected asymmetries by 1.4%. The measured asymmetry for each 



detector was calculated from the ADC averages by: 

_ Y.i{R-L),-T.^{BKG)R.L 

Ameas. J^^^R + L), - j:^(B KG) n^,' ^ ^ 

where the sum is over the central five channels including the M0ller peak. The (R-L) 
background was estimated by averaging the channels far from the M0ller peak. The (R+L) 
background subtraction increased the measured asymmetry by 17-24%. The full covariance 
matrix calculated from the pulse-by-pulse data was used to determine the statistical error 
of Ameas- The beam polarization was calculated from the measured asymmetry divided by 
the analyzing power. 

To check for possible systematic biases in the single arm analysis, both the number of 
channels included in the sum over the M0ller peak and the shape of the background fit to 
the unpolarized (R+L) line shapes were varied. From the observed spread in calculated 
polarizations, the sensitivity of the single arm analysis to the choice of fit parameters 
was estimated to be 1.3%. The polarization determined individually by each detector 
agreed with the overall mean within statistical errors. The total systematic error of the 
single arm analysis includes contributions from: the detector analyzing power, known to 
~ 0.5%, possible nonlinearities in the preamplifier and ADC response which could change 
the computed polarization by < 1.0%, and the sensitivity to analysis parameters, 1.3% as 
discussed above. The total systematic error of the single arm analysis is estimated to be 
1.7%. 



5. Results 

The polarization values measured by the single and coincidence arm polarimeters were 
in good agreement, although the results from the coincidence system were on average 0.6% 
lower than the single arm. Both systems measured the same polarization dependence on 
the quantum efficiency of the polarized source, resulting in a linear decrease of source 
polarization with increasing quantum efficiency. The polarized source quantum efficiency 
time history is shown in Fig. ^. The variations in polarized source quantum efficiency were 
related to the frequent cesiation treatments which were applied to the source in order to 
maintain the source quantum efficiency at an acceptable level. 

To obtain the beam polarization for the main analysis, the average polarization value 
for each run set was computed separately for each polarimeter. The results from both 
polarimeters are shown in Fig. ^ plotted as a function of source quantum efficiency. The 
plotted errors are a combination of the computed statistical errors and an additional 0.8% 
systematic error to account for non-statistical fluctuations in the data. 

A linear fit to the single arm and coincidence data as a function of polarized source 
quantum efficiency (QE) yields: 

Pb = (0.866 - 0.34 x QE) ± 0.003 ± 0.022, (36) 
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where the first error term is statistical while the second and dominant term is systematic. 
The systematic error includes a contribution of ±0.8%, as discussed above, a ±1.6% 
contribution from the average of the single arm and coincidence M0ller systematic errors, 
and a ±1.7% contribution from the uncertainty assigned to the target foil polarization. The 
resultant systematic error is ±2.5%. 



D. The Polarized Target 



The polarized target required a high-power ^He evaporation refrigerator operating near 1 
K, and a 5 T superconducting split pair magnet.0 The target material, frozen ^^N ammonia, 
was polarized using dynamic nuclear polarization (DNP). A schematic diagram of the target 
is shown in Fig. ^ [^. The magnet is shown with its field direction along the beam 
momentum direction. The refrigerator is positioned vertically and along the axis of the 
magnet. It is connected to a large Roots blower pumping system. The target insert lies along 
the central axis of the refrigerator. This insert was slid up and down to position any one of 
four targets in the beam. The targets were (from the top position) ^^NDa, ^^NHs, an empty 
cell, and either carbon or aluminum. A fifth position having no target was also available. 
The target insert also carried coaxial cables for the NMR measurement, a waveguide to 
transmit microwaves to the target (s) for DNP, and various temperature sensors. A diagram 
of the target insert is shown in Fig. M. 



1. DNP and Ammonia 



The DNP process for polarizing protons, deuterons, or any nucleus possessing a magnetic 
moment, requires temperatures of ~1 K or less and large magnetic holding fields. For 
thermal equilibrium at 1 K and 5 T, the proton polarization is only about 0.5%. However, 
the polarization of the "free" electrons, associated with the paramagnetic radicals introduced 
into the target material, is greater than 99%. The electron polarization can be transferred 
to the proton through a hyperfine transition by irradiating the target with microwaves at 
appropriate frequencies. The two polarization directions for the proton are reached by 
irradiation at frequencies slightly above or below the electron Larmor frequency, ~ 140 
GHz at 5 T. Details of the DNP process can be found in the literature, e.g., Abragam and 
Goldman ||90| or Borghini ||91[]. In our case the magnetic field was held at 4.87 T to match 
the frequency range (136 - 137 GHz) of the microwave tubJ being used. 

Ammonia was chosen as the target material because of its relatively large dilution factor 
compared to most other polarized target materials, its high polarizability, and its resistance 
to radiation damage. Furthermore, ^^N ammonia (spin ^) was chosen over ^^N (spin 1) 
because in ^^N the spin is carried by an unpaired proton, in contrast to ^^N where the 
spin is carried by a proton-neutron pair. Using ^^NHa reduces the systematic errors on the 
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proton spin structure functions by eliminating unwanted contributions from the neutron 
asymmetry. In addition, the ^^N polarization is easier to measure. 

The ^^NHa and ^^NDs targets were both prepared in the same way: First, the ammonia 
gas was slowly frozen in a test tube; the resulting solid lump of ammonia ice was crushed 
while immersed in liquid nitrogen and sifted to select granules of approximately 2 mm size. 
Smaller pieces were recycled in the same apparatus. 

The paramagnetic radicals necessary for DNP were introduced by irradiation using 
various electron beams. Each sample was immersed in liquid argon and given a dose of 
about 3 — 5 X 10^^ electrons cm~^. Targets for E143 were irradiated at Bates (at an electron 
energy of 350 MeV), at the Naval Postgraduate School, Monterey (65 MeV), and at the 
High Energy Physics Laboratory at Stanford (30 MeV). Samples of ^^N ammonia were 
irradiated at Saskatoon (250 MeV) for the initial tests. All irradiated samples were packed 
into thin-walled (0.0127 cm) torlon cylinders with 0.0025 cm aluminum end-cap windows. 
Each cylinder was 3 cm long and 2.5 cm in diameter, and contained two NMR coils made 
from 70/30 Cu/Ni tubing of 0.5 mm outer diameter and 0.0178 cm wall thickness. A straight 
piece of tubing was used to measure the proton polarization in the NH3 cell and the residual 
proton polarization in the ND3 cell. A coil of three to four turns with a 1 cm diameter 
measured the deuteron polarization and ^^N polarization in the ND3 target, while a similar 
one measured the ^^N polarization in the NH3 cell. During the course of E143 only the 
proton and deuteron polarizations were measured; the ^^N and residual proton polarizations 
were checked after the experiment. 



2. Polarization Measurement and Performance 



The polarization was measured via NMR with a series-tuned Liverpool Q-meter |^ 



Each spin species in the targets was measured with its own separately tuned Q-meter. Only 
one Q-meter could measure at a given time, taking one polarization measurement per minute. 
The Q of the tuned circuit is changed by the presence of the appropriate polarized nuclei, 
and the integral of this response is proportional to the polarization. The response function 
was determined by subtracting the Q-curve measured when the magnetic field was moved off 
resonance from the Q-curve obtained when the magnetic field was moved on resonance. The 
integral is normalized by comparing to the signal area at thermal equilibrium (TE) where 
the polarization (Pte) can be calculated. For the proton, Pte = tanh , where /i is the 
magnetic moment of the proton and k is Boltzmann's constant. Therefore, Pte = 0.0034 
for 5 = 5 T and T = 1.5 K. 

The TE signal for the proton is relatively easy to observe and measure, but the deuteron 
TE signal is about 500 times smaller than this, and thus requires advanced techniques of noise 
and drift suppression and signal averaging for a credible measurement Measurements 
were made of the proton TE signal area to a precision of about 0.2%, but repeated 
measurements, over a period of many weeks, showed considerable fluctuations in the mean 
value. Including this scatter, the overall precision of measuring the TE polarization was ± 
2.5%. For the deuteron the precision of measuring the TE signal area was ± 3% and ± 4% 
overall. The fluctuation in signal area was attributed to small changes in the distribution 
of ammonia granules around the NMR coils. Typical TE signals for polarized protons and 
deuterons are shown in Fig. 0. 



23 



In our initial measurements of ^^N ammonia, the proton polarization performance was 
similar to that seen previously ||9^ (>90%), but the deuteron only reached 13%. The 
maximum deuteron polarization was expected to increase with in situ irradiation p5[ . 
Figure 0, which verifies this expectation, shows how the proton and deuteron polarizations 
performed as a function of beam dose for 5 x 10^^ electrons/sec rastered over the face of a 
target. 

The polarization decayed with beam dose as the ammonia became radiation-damaged. 
Once the polarization fell below a predetermined value, the other target was put into the 
beam until its polarization dropped to a specified level. Then both targets were annealed 
by warming them up to a temperature between 80 and 90 K. The sequence of polarization, 
irradiation and annealing affected ^^NHs and ^^NDs differently. For the proton, annealing 
brings the target back to its starting polarization, and there was no evidence of change over 
the period of irradiation. On the other hand, the deuteron polarization improved after each 
anneal, ultimately reaching a polarization of 42%. This value was obtained with frequency 
modulation of the microwave source. Previously, at CERN, in the SMC experiment ||96|| , 



frequency modulation had been found to improve the polarization of deuterated butanol by 
almost a factor of two. In deuterated ammonia the gain is more modest, with a factor of 
two improvement in the rate of polarization which leads to a gain in absolute polarization 
of 3 — 5%. The level of proton polarization in Fig. |1^ is lower than the expected maximum 
of more than 90% seen in the early measurements. This was because the NH3 target was 
situated below the ND3 target which absorbed some fraction of the microwave power. 

After the experiment the ^^N polarizations were measured as a function of both proton 
and deuteron polarizations. Residual proton polarizations were measured in the deuteron 
case. In addition, the protons in the torlon target cups became polarized once the electron 
beam created paramagnetic centers in that material. This led to a 3% correction of the 
proton polarization. The polarization values were also corrected for effects arising from 
inhomogeneities in target polarization due to local beam heating as discussed below. 

The average polarizations for the entire experiment were 0.70 with a relative precision 
of 2.5% for the proton and 0.25 with a relative precision of 4% for the deuteron. 



3. Beam Heating Corrections 



As the beam passes through the polarized target, the temperature of the ammonia 
granules increases, and the polarization drops. By rastering the beam over the face of 
the target, this depolarization effect is greatly reduced. The average polarization measured 
by the standard NMR technique is generally not the same as what the beam sees locally. 
One reason for this is that ammonia granules outside the raster radius do not experience 
the same depolarization from beam heating as the granules inside the raster radius. The 
measured polarization, however, reflects a combined polarization of all the target granules. 
Another reason is that the polarization during the beam spill may be lower than during the 
time between spills when no beam heats the target. This latter effect has been studied in 
detail ||^ and has been shown to be very small. Hence, it has been neglected in the present 
analysis. 

If z is the relative contribution of the rastered granules to the NMR signal, then 
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= zPt + (1 - z)P,, 



(37) 



where Pi is the initial polarization with no incident beam, P^ is the measured polarization 
with incident beam, and Pt is the true polarization of the rastered granules. We define a 
correction to the measured polarization Cheat as 

1 — Cheat = Pt/ Pm = ^75 ^ • (38) 

The parameter z depends on the geometry of the NMR coils, which is different for NH3 and 
ND3 targets, and on the direction of the target polarization (longitudinal or transverse). 
Values for Pi, Pm, Pt, z, and Cheat and the corresponding errors are given in Table |l|for a 
maximum beam intensity of 4 x 10^ electrons/pulse. The errors on z include uncertainties 
for the target granule settling effect and for the rastering radius due to the finite size of the 
beam spot. For the ND3 targets there is an additional uncertainty in the diameter of the 
4 turns of the NMR coil. The corresponding corrections at other beam intensities can be 
extracted using the knowledge that the measured target depolarization is proportional to 
beam intensity. 



E. Spectrometers 

Two large acceptance spectrometers, situated at 4.5° and 7.0°, were used to detect 
the electrons scattered from the polarized target. The momentum acceptance of each 
spectrometer arm ranged from 7 to 20 GeV/c. Each spectrometer contained two dipole 
magnets, bending in opposite directions in the vertical plane, two gas threshold Cerenkov 
detectors, two scintillation hodoscope packages, each consisting of several planes, and an 
array of lead glass total absorption shower counters which were 24 radiation lengths in depth. 
The 4.5° spectrometer also contained a quadrupole magnet which was needed to spread the 
scattered electrons over a larger detector area. A schematic of the spectrometers is shown 



in Fig. ^ 



The two-bend design was chosen to have maximum acceptance over a wide momentum 
range, and to shield the detectors from the considerable photon background produced by the 
electron beam interacting in the thick target. The Cerenkov detectors allowed discrimination 
against a large pion background. The hodoscopes were used to reconstruct the trajectory of 
each particle, which in turn could be used to determine the momentum and other kinematic 
variables. Finally the shower array provided the energy measurement as well as particle 
identification information. 



The spectrometers were almost identical to those used in the E142 experiment [10 
The magnets, however, were operated at somewhat higher fields to accommodate the larger 
momenta of scattered electrons due to the higher beam energy. Also, the hodoscopes were 
modified to handle a higher instantaneous rate. 

The two scintillator hodoscope arrays provided the track information of the incident 
particles for each spectrometer, and consisted of horizontal {y), vertical (x), and slanting 
(u) planes of fingers. The upstream hodoscope array contained four planes: u, x, and two y; 
the downstream hodoscope array contained an x, y, and u plane. The first y plane in each 
hodoscope array consisted of scintillator elements of 3.0 cm width, while the elements of the 
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second y plane in the upstream hodoscopes were 4.76 cm wide. The elements in the first x 
plane of the 4.5° spectrometer were 2.0 cm wide, and the remaining x plane elements were 
3.0 cm wide. The u planes contained elements that were 4.5 cm and 7.5 cm wide for the 
front and back hodoscopes respectively. Within each plane the hodoscope fingers overlapped 
by 1/3 of the width on both edges, resulting in a bin width of 1/3 of the element width. 
The moderately fine hodoscope segmentation (~210 scintillator elements per spectrometer) 
was chosen to tolerate the large photon and neutron backgrounds and to reconstruct with 
sufficient resolution the trajectory of the scattered particles. The signal from each finger was 
discriminated and fed into a multi-hit TDC which recorded all signals in a 100 ns window 
around each trigger. 

The separation of the two hodoscopes was 5.0 m in the 4.5° spectrometer arm and 5.1 m 
in the 7.0° arm. The scattering angle resolution at the target in the non-bend plane was 
0.3 mrad for both spectrometers, whereas for the bend plane, it was ±0.9 mrad for the 4.5° 
arm and ±0.3 mrad for the 7.0° arm. The ideal momentum resolution was dependent on the 
absolute value of momentum and varied from ±0.3% to ±3.2% for the 4.5° arm and from 
±0.6% to ±3.8% for the 7.0° arm. 

The upstream Cerenkov counters were 2.24 m long aluminum tanks filled with nitrogen 
gas at a pressure of 6.3 psi for a pion threshold of 9 GeV, and the downstream counters were 
4.3 m tanks containing nitrogen at 3.0 psi for a 13 GeV pion threshold. Pions below these 
threshold momenta did not emit Cerenkov light. The shorter tanks had inner radii of 60 
cm and effective radiator lengths of 2.0 m, while the larger tanks had inner radii of 80 cm 
to cover the large spectrometer acceptances, and had effective radiator lengths of 4.0 m. To 
minimize 5-ray production and multiple scattering effects, thin tank entrance/exit windows 
were made from 1 mm thick aluminum. 

Inside the tanks, spherical mirrors were positioned to reflect all of the emitted Cerenkov 
light back onto a single Hamamatsu R1584-01 five-inch photomultiplier tube coated with a 
p-terphenyl wavelength shifter and maintained at a base voltage of -2600 V. The mirrors had 
a radius of curvature of 1.63 m and 1.2 m for the large and small counters, respectively, and 
had reflectivity close to 90%. The large counters contained three mirrors vertically stacked 
and mounted on an adjustable frame for focusing purposes, and the small counters contained 
two mirrors mounted similarly. Signals from each photomultiplier tube were discriminated 
at four levels corresponding to 0.6, 1.5, 3, and 4 photoelectrons and fed into four channels 
of multi-hit TDC's as well as an ADC. 

An electromagnetic shower calorimeter was positioned at the downstream end of each 
spectrometer. Each detector consisted of two hundred 6.2 x 6.2 x 75.0 cm Schott type F2 
lead glass blocks stacked 10 wide and 20 high in a fly's eye conflguration. The glass had 
a radiation length of 3.17 cm and refractive index of 1.58. The incident electrons created 
showers via bremsstrahlung and /e" pair production in the lead glass. Electrons (and 
positrons) produced Cerenkov light in an amount proportional to the incident energy. The 
light was collected by phototubes attached to the back of the glass. To monitor the blocks, a 
high intensity Xe flash lamp system (Hamamatsu L2360) was installed in each calorimeter. 
The lamp delivered a luminous signal to each block via plastic optical flbers. The signal 
from one of the fibers, as well as that from a ^^^Am source, was read out by a monitoring 
photomultiplier tube and sent to ADCs to detect possible Xe lamp intensity fluctuations 
and to monitor ADC gain changes by looking at shifts in averaged signals. 
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F. Trigger 



The trigger consisted of a triple coincidence between discriminated signals from the 
two Cerenkov counters and the analog sum of the shower counter elements. The shower 
discriminator threshold was set to be greater than 99% efficient for the lowest energy 
electrons and the Cerenkov thresholds were set to be efficient for one photoelectron signals. 
Up to four triggers could be generated in each beam spill. Each shower and Cerenkov 
counter signal was fanned out to four separate ADCs, and each trigger gated a different set 
of these ADCs. The detector signals to the multi-hit TDC's were filtered by a sub-trigger 
to reduce noise hits. Additional triggers were used to record a small fraction of the pions 
and to measure detector efficiencies. 



G. Data Acquisition 



The data acquisition (daq) was distributed over a number of computers linked together by 
an ethernet-based network which implemented DECnet for communications. The distributed 
nature of the daq allowed us to build a system which could service interrupts at 120 Hz, 
read typically 3 KB of data for each interrupt, write data to tape at a sustained rate of 
nearly 300 KB/sec, control the electron beam position on target on a pulse-to- pulse basis, 
and analyze a substantial fraction of the event data online. 

A VAX4000.200 qbus computer, referred to as real time front end (RTFE), was 
interrupted at 120 Hz, read data from three CAMAC branches, built an event, and then 
sent it via network to the data logger computer, a VAX4000.60 workstation. The RTFE ran 
an application which was developed with DEC's VAXeln development toolkit. The VAXeln 
application was able to access hardware resources more efficiently than usually possible 
under the DEC VMS operating system. Also, task scheduhng was under programmer 
control. The data logger computer controlled two SCSI EXB8500 Exabyte tape drives. 
Event data received from the RTFE were packed into record-size buffers (approximately 32 
KB in size) and written to tape. The data logger computer also distributed a sample of 
the event data via network to two VAX4000.60 workstations, one for online data analysis 
in each spectrometer. The event data were analyzed and various histograms and tables 
were presented for viewing in Xll/Motif windows. Special purpose analyses could be 
performed by other VAX workstations which connected to the network. The electron beam 
was monitored and controlled by a micro VAX II computer which, like the RTFE, was 
loaded with a VAXeln application specially developed for this task. A VAX cluster boot 
node, VAX4000.300 computer, was used to control, monitor, and log information on the 
spectrometer magnets and their power supplies, detector high voltage power supplies, NIM 
and CAMAC crate voltages, scalers, target parameters, pedestals, etc. 
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IV. DATA ANALYSIS 



A. Event Selection 



Events which produced a trigger were further analyzed to identify electrons amid a 
background of mostly pions, and to determine the energy, momentum, and scattering angle. 
Raw asymmetries were then formed from the number of scattered electrons coming from 
each of the two states of incident electron polarization directions. These asymmetries are a 
function of x, Q^, and beam energy. 

Particle tracking was performed using the spatial and timing information provided by 
the hodoscopes and shower counter. Once a track was found, the particle's momentum and 
scattering angle were reconstructed. The efficiency of the hodoscope package was found to 
be 91% for the 4.5° spectrometer and 96% for the 7.0° spectrometer. The tracking efficiency 
was about 98% for the 4.5° and 99% for the 7° spectrometer. The hodoscope and tracking 
efficiencies were worse for the small angle spectrometer due to the higher count rate. 

The shower counter was used to measure the energy deposited by the incident particle 
and to provide electron identification. In order to use the shower counter for energy 
measurements, it was necessary to calibrate each block for differences in phototube, lead 
glass, and ADC channel responses to the electrons. This was achieved using clean electron 
events which were selected using knowledge from the other detectors. In an iterative process, 
a set of calibration constants for the glass blocks was determined by requiring that the total 
energy of the cluster be equal, on average, to the momentum of the event. 

Once calibrated, the shower counter was used to select electrons by comparing the energy 
of the particle as measured by the shower counter (£") to the momentum of the particle 
as measured by the hodoscope tracking system (P). Rejection of pions was achieved since 
typically electrons deposit all of their energy in the shower counter while pions do not. Thus, 
the electron events had an E'/P peak centered around unity, whereas the E'/P values for 
pions were in general much less than one. By making a cut around the electron peak of 
0.8 < E'/P < 1.25, we were able to reject the majority of pion contaminants left in our 
data sample. The E'/P requirement was approximately 96% efficient for electrons, and left 
a pion contamination of less than 1%. A sample plot of E'/P for this experiment is shown 
in Fig. 0. 

Electrons could also be identified over pions using the spatial profile of the shower formed 
by the incident particle. These profiles vary significantly depending on the type of incident 
particle. In particular, pion shower clusters are much smaller than electron clusters, and 
many of them are fully contained in one block. Electron clusters are typically contained in 
nine glass blocks. For electrons, the central block contained 50% to 90% of the energy and 
the eight neighbors contained the rest. We used a shower counter neural network algorithm 
| 99| which modeled a typical electron cluster profile to determine which events were electron 
events and which were background. The neural network was approximately 98% efficient for 
identifying electrons and left a pion contamination of about 0.5%. 

To further ensure a clean electron sample, spatial, timing, and pulse height cuts were 
made. A cut requiring a minimum Cerenkov ADC pulse height of 40 (more than two 
photoelectrons) was made yielding an efficiency of 95-99% for the four counters. Next, 
we required that the track used for the momentum measurement was within 40 mm, 
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horizontally and vertically, of the cluster from the shower counter that was used for the 
energy measurement. The track and the cluster were required to be within 10 ns of each 
other, and the track was required to point back to the target to within 13 mm to eliminate 
bad tracks. In addition, the few events with clusters on the outer edges of the shower counter 
were rejected due to the possibility of energy leakage from the sides of the counter. 



B. Asymmetries and Corrections 

Electrons passing the event selection cuts were binned in x such that the resolution in 
X was slightly finer than the binning. The electrons were also tagged according to their 
relative target and beam helicity states, iV^-'-^^^^ and which spectrometer they entered. The 
asymmetries A\\ and A±_ were formed: 

A,(o.^J^a(^^-C.)+A. (39) 

Here / is the dilution factor, and Pt are the beam and target polarizations, A^c is the 
radiative correction to the asymmetry, and Ci(2) are the corrections needed due to the 
presence of nitrogen in the targets, with C2 disappearing for the proton target. These 
corrections are discussed in more detail below. Here A^l(_r) is the number of left or 
right-handed helicity events corrected as 

^^(^^ = <m^^ (40) 

where c?l(_r) is the appropriate dead time correction and Ql(r) is the appropriate incident 
charge. 



1. Polarized Nitrogen and Residual Proton Corrections 

In measuring the proton and the deuteron asymmetries, it was necessary to correct for 
events which scattered from other polarizable nuclei in the target aside from the desired 
protons or deuterons ||100||. The targets were made of ^^NHs and ^^NDs, and both the ^^N 



and the ~2% contamination of ^^N were polarizable. In addition, the ^^NDs target contained 
~1.5% of unsubstituted or residual polarizable protons from ^^NHs. 

The polarization of ^^N and the residual protons was measured after the experiment. The 
unpaired proton in ^^N contributes to the measured proton asymmetry proportionally to the 
nitrogen polarization and with a negative sign because of the negative magnetic moment of 
^^N. For the target material ^^NHa, the following fit was used to express the ^^N polarization 
Pn in terms of the polarization of the protons Pp'. 

Pn = 0.136Pp-0.183Pp2 + 0.335 ^ 0.12. (41) 

The correction Cf to the proton asymmetry (Cf = 0) which is referred to in Eq. (|39|) is 
given by: 
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Cl = l- ll^^EMc(x) ^ 0.98. (42) 



Here g-EMc{'x) is the correction for the EMC effect ||101|| taken at atomic mass number 15. 
The first factor — | comes from Clebsch-Gordan coefficients involving the nitrogen wave 
function. The second factor | reflects the fact that ammonia has three hydrogen atoms for 
each nitrogen atom. The error on the second term in Cf was estimated to be about 20% 
relative which yields a systematic error of 0.004 on Cf . Here the contribution of ^^N to the 
asymmetry was neglected. 

For the target material ^^NDs, the corrections were more complicated because they 
account for both the residual protons and the unpaired proton in the ^^N. For each case, 
the correction involved the measured proton asymmetry. 

The ^^N polarization Pjq is given by 

= -0.40Prf, (43) 

where Pd is the polarization of the deuterons. The residual proton polarization is expressed 
as: 

P^^' = 0.191 + 0.683 Pd for Pd > 0.16 

= 1.875 Pd for Pd < 0.16 (44) 

The polarization of ^^N was obtained from the measured ^^N polarization by assuming that 
the polarization was equal and opposite in sign to that of ^^N. The corrections used in 
Eq. (p9D for the ^^NDs target are given by: 



(jd — ~ 1 02 

TJ 

Ci = TfrK^n - Dp){A'' - Arc) ~ -0.03(A^' - A,,). (45) 

Ap is the final proton (Born) asymmetry A\\ or A±, and by subtracting the appropriate 
proton radiative correction Arc we are left with the radiated asymmetry. Up and Ud are the 
radiative corrections to the unpolarized cross-sections. The remaining factors are defined as 

number of protons 015 

number of deuterons + number of protons ' 



Dn = Vn 



Pd 9 



pres p /\ 

D, = V^ir- + (2^ - 1)^^^^ , (46) 



number of ^"^N 
number of ^^N + number of ^^N 



The error on Cf was neglected since this value was very small and stable. The factor C2 
contains the proton asymmetry and was calculated for each x-bin using the measured proton 
asymmetry and its error. 
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2. Background Subtraction of Positrons and Pions 



The data collected in each of the spectrometers included background events coming 
from a small number of misidentified pions and from electrons produced in pair-symmetric 
processes (mostly 7r° — > 27, 7 — > e~ + e"'"). This background (mostly pair-symmetric) was 
responsible for up to 10% of the events in the lowest x-bin, h\it close to zero events for x > 0.3. 
To measure the background, data were taken with the spectrometer magnets' polarity flipped 
to measure tt"*" and e"*". The same cuts were applied to eliminate the majority of pions as in 
the electron runs. A positive particle asymmetry was formed and was corrected just as 
in the case of the electron asymmetry A_ for varying experimental conditions such as beam 
and target polarizations. This positive particle asymmetry was found to be consistent with 
zero. The background-corrected electron asymmetry was determined by 

A = A_ = A+ — — , (47) 

where N_, are the number of events per incoming charge for electron and positron 
runs. The misidentified pion background was subtracted along with that of the positron 
background since A+ also contained a measure of the misidentified pions and assuming 



3. False Asymmetries 



It is important to make sure that our experimental data are free from significant false 
asymmetries which could systematically shift the data. During the experiment, data were 
taken (either longitudinal or transverse) with the target B-field pointing in either one of the 
two possible directions. For each field direction, two different target polarization directions 
were used, parallel or antiparallel to the B-ficld. Wc then had four different configurations, 
and approximately the same amount of data were taken for each configuration, thus 
cancelling out the electroweak contributions to our measurement so that no correction to the 
data was necessary. The asymmetries in each configuration were compared by looking at the 

distributions of the asymmetry differences. For the proton, the distributions were all 
nicely centered at one, and the mean value of the asymmetry differences was approximately 
one standard deviation from zero. For the deuteron the results were slightly worse, yet still 
very reasonable. The distributions were centered around one with a few points greater 
than two, and the mean values of the asymmetry differences were within two standard 
deviations of zero. We conclude that there were no significant systematic effects on the 
asymmetry due to changes in target B-field or target enhancement field directions. Also, no 
statistically significant variation of the asymmetry was found for either NH3 or ND3 targets 
as a function of raster position. 



4- Dilution Factor 



In general, incident electrons will scatter both from polarized target nucleons and 
unpolarized nuclei that are part of the target assembly. These unpolarized materials include 



31 



liquid helium, ^^N, NMR pick-up coils, and vacuum windows. Scattering from unpolarized 
materials will dilute the measured asymmetry, and a correction must be applied. The 
dilution factor / is a function of x and Q^, and is defined as the ratio of the total event rate 
from polarizable nucleons, to the total event rate from all target materials. The measured 
asymmetry is then corrected for unpolarized events by dividing by the dilution factor. 

For a material of density p and thickness z, the event rate from Born processes was 
calculated as follows: 

r(x, Q2) = pz [ZF^ix, + NF^ix, Q^)] gEMc{x, Q'), (48) 

where Z and N are the number of protons and neutrons in the nucleus. and are 
unpolarized proton and neutron deep-inelastic structure functions. They were obtained 
from a parameterization of the NMC data [|102|] . The factor g corrects for the "EMC effect" 
which accounts for the difference in nucleon cross sections for free and bound nucleons. 
With this model for rates, the dilution factor can then be calculated as follows: 

where we are summing the rates from all unpolarized materials which contribute to the 
overall event rate. The factor Tc corrects the dilution factor for radiative effects and was 
typically less than a 5% correction. 

The target material was in the form of frozen granules which were tightly packed into a 
target cell. The volume fraction of the target cell which the target material occupied is known 
as the packing fraction, and was determined independently by three different methods. The 
first consisted of studying the difference in event rates for empty, carbon, and full target 
cells. The second method was a measurement of the attenuation of a mono-energetic X-ray 
beam as it passed through the target material. The attenuation of the incident beam is 
directly related to the thickness and attenuation coefficient of material it passes through, 
and was therefore sensitive to the packing fraction. Finally, the target material was weighed, 
and the packing fraction was determined using the known volume of the target cell. The 
measured packing fraction was different for each target used, and varied from 0.57 to 0.64. 

Over the kinematic range of interest, the dilution factor typically varied from 0.15 — 0.19 
for the NH3 target and from 0.23 — 0.25 for the ND3 target, with relative errors of 2% and 
1.5%, respectively. The error on the dilution factor comes from several sources. The packing 
fraction was known to 4% for both targets. The relative error from the cross-section ratio 
Cn/cTp was 1.0% | ]103|| . This was one of the dominant errors for the NH3 target and did not 



contribute to the ND3 target. The ratio of nuclear to deuterium cross sections, the EMC 
effect, is known to 1.5% relative and was another large source of error. This effect has a 1% 
overall normalization and another 1% uncorrelated error ||101|| . The small mass of the NMR 



coil (~0.1 gm) was known to 20% but did not contribute significantly to the overall error. 

5. Dead Time 

All the signals from various detectors went through discriminators before forming the 
various triggers. These discriminators have an output pulse width of 25 ns and a double 
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pulse resolution of 8 ns. They were operated in an updating mode such that a second 
signal entering the discriminator after 8 ns and before 25 ns produced a single output pulse 
with an extended width. The effective dead time was 32 ns instead of 25 ns, due to slight 
mis-timing between various signals and signal jitter, especially from the shower counters. 
Rates measured with each beam helicity were corrected separately. 

Using a Monte Carlo simulation, the probability matrix M{i,j) for observing i hits when 
there were really j hits was generated using a typical spill length of 2200 ns and a dead time 
of 32 ns as inputs. The observed trigger frequency distribution To(i) is related to the true 
distribution (without dead time) Tt{i) by 

T„(z) = ^M(z,j)xr,(j). (50) 

j 

Since there are practically no hits beyond 10, the sum is safely truncated at j = 16. The 
matrix M is inverted to solve for T((i), and the dead time correction factor d to the measured 
rates is defined as 

Here, the sum over Tg is split into two parts because only four triggers could be recorded 
per beam spill. 

The correction factor varies smoothly from 1 at very low rates to 1.07 at an average rate 
of 2 events/pulse. The systematic error on the corrections was calculated assuming upper 
and lower limits to the beam width of 2600 ns and 1800 ns. The dead time was found to be 
accurate to a few parts in 1000, and the error for the corrected asymmetry by applying these 
factors is found to be less than 2 x 10~^, which is completely negligible. No uncertainty in 
the dead time itself was considered because only the ratio of the beam spill length to the 
dead time is important. 



6. Radiative Corrections 

Our experimental goal was to measure a single photon exchange process (Born) at specific 
kinematics. In reality there are higher order contributing processes (internal), and the actual 
scattering kinematics can change due to energy losses in materials along the electrons' paths 
(external). The radiative corrections account for these unwanted effects. 

The radiative correction calculation is different for the unpolarized (cr") and polarized 
((jP) components of the helicity-dependent cross sections which are given by a^^ and a^^ for 
a longitudinally polarized target, and and a~^'^ for a transversely polarized target. The 
longitudinal and transverse asymmetries can be written as 



(52) 



which is equally valid for Born, internally radiated, or fully radiated cross sections and 
asymmetries. For the remainder of the radiative correction discussion, quantities which are 



33 



Born, internally radiated, or fully radiated are subscripted with 0, r, and R respectively. 
Also, for simplicity, references to a particular target polarization are dropped such that A 
could be either A\\ or A±, and could be either a^^ or cr^^. 

Calculation of the polarization-dependent internal corrections was done using code based 



on the work of Kuchto, Shumeiko, and Akusevich |p.04|| , who are also responsible for the 
development of their own code POLRAD. The two codes were carefully compared and found 
to be completely equivalent when the same input models were used. The calculation of the 
internally radiated cross sections can be decomposed into components: 

< = ao^(l+5.)+a^, + < + af„, 

< = ao"(l + 5.) + + < + <, (53) 

where 5^ includes corrections for the electron vertex and vacuum polarization contributions, 
as well as a term that is left after the infrared- divergent contributions are cancelled out. The 
vacuum polarization correction includes contributions from both leptons and light quarks. 
The terms aei, Cq, and ain are the radiative tails due to internal bremsstrahlung (occurring 
within the field of the scattering nucleon) for elastic, quasielastic, and inelastic scattering 
processes. 

The nuclear elastic tail for the deuteron was evaluated using fits to data | |105| | over a large 
range in Q^. The elastic tail for the proton (and quasielastic for the deuteron) were evaluated 
using various form factor models ||106|| which agree well with existing data over the kinematic 



region. Note that some of these models agree well with data for some of the four elastic form 
factors and not others, so different models were combined for the best representation of all 
four nucleon elastic form factors. Quasielastic cross sections were Fermi-smeared only for 
corrections to our resonance data since this smearing has a negligible impact on the radiative 
correction in the deep-inelastic region. Unpolarized cross sections were modeled using fits 
to structure function data in the deep-inelastic region ||102| , [78[| , and fits to cross sections in 



the resonance region | |107| , |108|| . The polarized component to the deep- inelastic cross sections 
was modeled using Q^-dependent fits to Ai as given in this paper. The polarized resonance 
region model was based on parameterizations of previous data and data presented here. For 
the transverse contributions, we used g2 = or (72 = which are both consistent 

with our data. 

The external corrections account for bremsstrahlung radiative effects which occur as the 
electrons pass through material in their path. Ionization effects were completely negligible 
at our kinematics. At any given interaction point within the target the radiative correction 
depends on the amount of material (in radiation lengths) the electron sees before (tb) and 
after (to) scattering. Because the radiation lengths before and after scattering did not 
vary significantly over the beam raster area, it was not necessary to integrate the external 
corrections over the raster area. Also, it was an excellent approximation to replace the 
target integration over the length of the target with the evaluation of the external radiative 
corrections at one point, namely the center of the target. At this point, U = 0.026, ta = 0.047 
for the 4.5° spectrometer, and ta = 0.040 for the 7° spectrometer. These radiation lengths, 
which are valid for both our proton and deuteron polarized target, are dominated by the 
target material, but also include contributions from various windows which are not part of 



the polarized target itself. The external corrections were thus calculated using ||109|| 
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aliEo,E'p,e) = I h{Eo,E,U)a^^{E,E\e)h{E\E'p,ta)[l-D{Eo,E,Z)]dEdE' 

al{Eo, E'p,e)= I h{Eo, E, h)a:;{E, E', e)h{E', E'^, QdEdE' (54) 

J Em 'in E m 



where Eq is the electron initial energy, E'p is the final scattered electron momentum, and 
Emin and E'^^^ are the minimum incident energy and maximum scattered energy as defined 
by elastic scattering. Ih{Ei, E2,t) is the probability ||109|| that a particle with initial energy 



El ends up with energy E2 after passing through a radiator of thickness t, and D[Eq, E, Z) 



is the electron depolarization correction [ p.lCI|| which corrects for the depolarization of the 



electron beam due to the bremsstrahlung emission of polarized photons. This correction 
depends weakly on the Z of the target material. 

An additive correction Arc to the data was formed by taking the difference between the 
fully radiated and Born model asymmetries 

Ac = Ao-An = ^-^. (55) 

Our fits to Ai and the radiative corrections were iterated until they converged. For 
the purposes of statistical error propagation on our measured asymmetries, a "radiative 
correction dilution factor" frc was evaluated. This dilution factor is simply a ratio of events 
coming from deep-inelastic processes to all events and multiplies the usual dilution factor in 
Eq. We only used frc for the error propagation and not for correcting the data directly. 
Systematic uncertainties were estimated by varying input models within reasonable limits 
and measuring how much the radiative correction changed. These uncertainties for the 
various models were then combined in quadrature for each x bin. Results for Arc are listed 
with data in Tables pKHVini, [XVI| - [XIX| , and |XXXV1| - [XXXVTT , 



C. Analysis of Resonance Region Data 

The resonance data [ |111|| were taken with a 9.7 GeV beam. The spectrometer angles 
of 4.5° and 7° corresponded to ~ 0.5 and 1.2 (GeV/c)^ in the resonance region 
[W"^ < 5 GeV^), respectively. We have extracted gi from the measured asymmetries A\\, 
and from the absolute cross-section differences given in Eq. ^ Each method has its own 
set of systematic errors. The difference method requires good knowledge of spectrometer 
acceptances, the number density of polarizable protons or deuterons in the target, and 
detector efficiencies. The asymmetry method requires knowledge of the dilution factor for 
the resonance region, which means an accurate model of the rapidly varying unpolarized 
cross sections is needed. We found that the two methods agreed to within a fraction of the 
statistical errors on each point (typically better than 3%). In our previous report we have 



used the difference method. The current reanalysis uses the asymmetry method, since we 
now believe that the systematic errors are slightly better in this case. Other improvements 



on Ref. [0 include better modeling of the resonance region for radiative and resolution 
effects. 

The resonance asymmetries were calculated as specified in Section IV. B for the 
deep-inelastic analysis. In the present case, we have determined the dilution factor / using a 
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Monte Carlo routine as described below. The term Arc also included a resolution correction 
in the resonance region. 

We developed a Monte Carlo code which simulated all relevant aspects of the experiment. 
It was used to predict total count rates and count rate differences from a set of tables 
of cross sections and asymmetries generated by the radiative corrections routine. The 
unpolarized cross sections came from parameterizations for the resonance region ||107|| and 
the deep-inelastic region [102|J75|] . The asymmetries contain both resonant and non-resonant 
contributions. The resonance contribution was calculated using the code AO ||112|| , which 



includes parameterizations of all of the existing resonance data; however, the helicity 

amplitudes Ai and As for 5*11 and D13 were tuned to agree with our data. The non-resonant 

2 2 

part came from a parameterization of all existing deep-inelastic data (Fit 111 of Ref. 0), 
which was extrapolated into the resonance region. Specifically, for W"^ < 2.5 GeV^, Ai 
was given by the tuned AO result alone; above W"^ = 3.0 GeV^, Ai was taken as the sum 
of the AO resonant contribution and the Fit III inelastic background; and in the region 
2.5 < W"^ < 3.0 GeV^ the two extremes were linearly interpolated. A2 was calculated using 
, which yields values close to zero. The model-dependence of this choice for A2 was 
determined by alternately considering g2 = and A2 = 0. 

The observed raw parallel asymmetry Ajj""' is proportional to the combination of photon 
asymmetries Ai + Therefore, we first extracted Ai + 7]A2 from the data, and then 

deduced gi from this using various assumptions about g2- The following steps were required 
to produce Ai + rjA2 and gy. 

1) The radiative corrections code was run with the options as specified above to create the 
Born cross sections, the Born asymmetries, and the predicted values of Ai, A2, gi and g2 
for both NH3 and ND3 targets at E143 kinematics. 

2) The radiative corrections code was run to create tables of cross sections and asymmetries 
over a wide range of kinematics, fully internally radiated, to use as input to the Monte Carlo 
generator. 

3) The Monte Carlo routine was run for both polarized protons and deuterons alone, and 
for full ND3 and NH3 targets. This simulation included external radiation, spectrometer 
acceptance, resolution, multiple scattering and Fermi motion, as well as the reconstructed 
kinematic variables and raw asymmetries. 

4) The raw data was corrected for efficiencies, polarization, polarized nitrogen and polarized 
protons in ND3 using the standard E143 procedure. Then, the data was corrected by the 
dilution factor (the ratio of Monte Carlo events from polarizable protons or deuterons 
to those from all target components), and the additive radiative correction term Arc 
(obtained from the difference between fully radiated Monte Carlo results and the model 
Born asymmetry) was applied in order to generate the fully corrected values of A\\. 

5) The ratios {Ai + riA2)/A\\ and gi/A\\ were used as calculated in the Born version of the 
radiative correction routine to find gi and Ai + riA2 for our data. 

This extraction method required that the Monte Carlo routine provide a detailed and 
realistic simulation of the data, including resolution effects which are very important in 
the resonance region. Therefore, we performed a series of tests to insure that the Monte 
Carlo simulation described the data well, and provided radiative and resolution corrections 
with sufficient precision compared to the statistical accuracy of our data. Without any 
normalization factors, the generated unpolarized counts versus W'^ agree with the data to 
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better than 2.2% (3.4%) in the 4.5° (7°) spectrometer. The shape as a function of W"^ is 
in even better agreement in the region of the resonances. This provides confidence that the 
acceptance and resolution of the spectrometer are properly modeled. In addition, we have 
found that the measured and simulated count-rate differences agree well with each other. 
These rates were integrated over the (quasi-)elastic region [W"^ < 1 GeV^) where model 
uncertainties are minimal because of reasonably good knowledge of elastic form factors and 
radiative corrections. The overall for the four degrees of freedom corresponding to p(4.5°), 
p(7°), d(4.5°) and d(7°) is 3.85. 



V. RESULTS AND COMPARISON TO THEORY 

Table |T| gives the results for y4j| and with the proton target for the beam energy of 
29.1 GeV and for the 4.5° and 7° spectrometers, respectively, along with the total radiative 
corrections to each point. Tables [rV[-|V| give the results for A^^ and radiative corrections for 
the beam energies of 16.2 and 9.7 GeV and for the two spectrometers. Tables [V^[VIII| give 



the corresponding results for the deuteron target. Since the 29 GeV data include both A\\ and 
Aj_, Eqs. ^are used determine the asymmetries Ai and A2, and Eqs. ^are used to determine 
the structure functions gi and g2 and the ratio of structure functions Qi/Fi. The NMC fit 
102|| was used for F2. The SLAG global analysis |]78[ was used for R. While the fit to R was 



made to data with a limited range and x > 0.1, it is consistent with recent measurements 
at lower x pi3| , |114 | and different [ |115| ]. Estimated errors on these unpolarized structure 



functions are given in Section |VB| . The neutron spin structure function can be extracted 
from the deuteron and proton results in a manner similar to that used for the unpolarized 
structure functions. For both gi and g2 we use the relation: 

^"(x, Q2) = 2g\x, g2)/(i _ 1.5^^) - /(x, (55) 
where ojd is the probability that the deuteron will be in a D-state. We use ojd = 0.05 ± 0.01 



116|| given by N-N potential calculations. No other nuclear contributions to ujd are included. 



The neutron asymmetries can then be calculated using Eqs. 



A. Ai and gi 

For beam energies of 16.2 and 9.7 GeV there are no Aj_ data available. We have assumed 
that g2{x,Q^) is given by either g^^ or g2=0, both of which are consistent with our g2 
data at 29 GeV. These different assumptions lead to very similar results. We have then 
determined Ai and gi using Eq. |. 

Tables |^ and ^ show the values of gi/Fi and Ai for deep-inelastic scattering (W^ > 



4 GeV^) for all three beam energies and both spectrometers using g2 = g^^ for the 16.2 



and 9.7 GeV data and A\\ and A± for the 29.1 GeV data. Figures |T3| and |T^ show gi/ Fi for 
proton and deuteron as functions of averaged into 8 x bins. Data from EMG 0, SMG 
0J||1|], SLAG E80 0, and SLAG E130 § are also included. The resuhs are consistent 
with gi/Fi and Ai being independent of for >\ (GeV/c)^. We fit all the deuteron 
and proton data, (including the SMG data at < 1) with the empirical parameterization 
ax"(l + hx + cx^)[l + C/(Q^)]. The coefficients of the fit are shown in Table and the fits 
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and |T^. We chose three forms for f{Q'^) shown in Column 3 of the 



are shown in Figs, 
table. They are 

• f{Q^) = 0: no dependence. 

• fiQ"^) = higher twist behavior. 

• /(g2) = ln(l/g2): pQCD behavior. 

The minimum of the fits is shown in Column 2 and is either 0.3 or 1.0 (GeV/c)^. 
For the proton data, the x^/df is less than unity for all the fits except fit I, indicating that 
there is dependence for < 1 (GeV/c)^. Fit II indicates that there is no need for any 
Q^-dependent term for > 1 (GeV/c)^, which is our cut-off for deep-inelastic scattering. 
For the deuteron data, the fits are not as good, but still have a confidence level of about 
10%. Fit V, which has Q^j„ = 0.3 (GeV/c)^, is used to evaluate and to iterate the 

radiative corrections described above. 



Also shown in Figs. ^ and |14| are the results from the E154 ||117|| leading order pQCD 
evolution fit to world data including preliminary results from this experiment. It is in good 
agreement with the data, including the data for < 1 (GeV/c)^ which was not used in the 
fit. However, for the proton it does have an exaggerated dependence at the highest x. 
Since gi/Fi and Ai are both consistent with being independent of for > I (GeV/c)^, 
we choose to combine our data at fixed x by averaging them over all measured values for 

> 1 (GeV/c)^. Tables pCII| and pCIIl| and Fig. |1^ show these averaged values as a function 



of X. The band at the bottom of Fig. represents the size of the systematic errors. Also 
shown are results from other experiments pj^,P,|Il] averaged in a similar way. The various 



experiments are in agreement with each other. Results for gi/Fi and Ai are similar at low 
X and diverge slightly at high x. For the proton, both gi/Fi and Ai are small and positive 
at low X and rise steeply toward unity as a; ^ 1. For the deuteron, both gi/Fi and Ai are 
close to zero at low x and increase slowly with increasing x. For the neutron, both gi/Fi 
and Ai are negative over most of the x region, showing almost no indication of becoming 
positive at high x as expected from earlier predictions p7 . 



Table pCIV| shows gi at the averaged measured value of obtained from the average 
value of gi/Fi. The quantity g^ was obtained using Eq. 



)6|. Figure |16| shows xgi as a 
function of log x. The area between the data and zero is the integral forming the sum rules 



/o ^(x)rfx = j^^xg{x)d\Yix. 



B. Systematic Errors 

The systematic errors were calculated for gi/Fi, Ai, and gi. Only the systematic error due 
to a" was considered since the systematic errors due to A± were negligible compared to the 
statistical errors. Some of the errors were multiplicative and independent of x while others 
were a;-dependent. The errors due to multiplicative factors (beam and target polarization) 
are shown in Table pCV| . The errors on gi and gi/Fi from these normalizations were obtained 
using a smoothed fit to gi/Fi. The breakdown of the major sources of error for a sample of 
our X bins is shown in Tables pCVl[ - [X.lX| for deuteron and proton targets for both gi and gi/Fi. 



The radiative correction error dominated at low x. The errors due to multiplicative factors 
were only significant when either gi or gi/Fi were large at middle and high x, respectively. 
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Multiplicative systematic errors: 



The error of the beam polarization Pb was estimated to be 0.024 (relative). See 
Section ffll C 5. 



The error of the target polarization Pt had a relative systematic error of 0.025 
for protons and 0.04 for deuterons. The error was assumed to be 100% correlated 
between runs, since the systematic error was obtained from the spread of the thermal 
equilibrium measurement results, each of which provided the calibration constants for 
large groups of runs. See Section |III D 2| . 



• The proton in nitrogen correction (see Eq. ^2D contributed with a 0.004 relative 
systematic error since the correction Ci was always around 0.02, while the relative 
error on Ci was estimated to be 0.2. 

Systematic errors dependent on x: 

• The error of the dilution factor / came from several sources. The component dependent 
on our experimental setup (the amount of ammonia in the target cell) was known to 



4% for both targets. The relative error from the cross-section ratio cx^/ was 1% ||103 
It was one of the dominant errors for the NH3 target and did not contribute for the 
ND3 target. The ratio of nuclear to deuterium cross sections, the EMC effect, is known 
to 1.5% relative and was another large error source. It has a 1% overall normalization 
and another 1% of uncorrelated error |101|. The small mass of the NMR coil was 



known to 20% but did not contribute significantly to the overall error. This leads to 
an average error on / of 2% for NII3 and 1.5% for ND3. 

The nitrogen correction was applied via two factors, Ci and C2 (see Eq. ^). The error 
on Ci was neglected since this value was very small and stable. The factor C2 (ND3 
only) contained the proton asymmetry and was calculated for each x bin using our 
measured proton asymmetry and its error. 

The systematic error on the radiative corrections was calculated for each x bin by 
varying several classes of input models. See Section [IV B 6| for details. It is shown in 



Tables PCVI[ - PCIX| for typical values of x at 29 Ge V beam energy for gi and gi / Fi 



The error due to the structure function R{x,Q^) contributed to gi/Fi and gi quite 
differently due to the relationship between Fi and F2. For gi the effect of R is negligible 
whereas for gi/Fi it is one of the significant errors. Its systematic error was taken from 
the SLAC global analysis ||78[ and ranged from 3% to about 7.5%. While this fit to 
R was made to data with a limited range and x > 0.1, it is consistent with recent 
measurements at lower x \TT^,TT^] and different | 115 ]. 



The error in the structure function F2 was obtained from the NMC fit [ 102 |. The error 
returned from the fit was taken as completely correlated point-to-point. 

Pion and charge-symmetric backgrounds were treated as statistical errors from the 
measurement with spectrometers set at opposite polarity. No systematic error was 
assigned to the model of charge symmetry. 
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When averaging gi/Fi over spectrometers and beam energies, the weight of each 
data point included statistical and point-to-point uncorrelated systematic errors. For the 
systematic error of the neutron structure function as well as of the difference Qi — Qi, the 
beam polarization error and the dilution factor errors due to the unpolarized cross sections 
were assumed to be 100% correlated, while the other errors were assumed to be uncorrelated. 

The systematic error on the integral was calculated using the separated correlated and 
uncorrelated systematic errors. The systematic errors of the low-x and high-x extrapolations 
were added together with the systematic error for the data region. The sum was then 
quadratically combined with the fit errors for the low- and high- a; extrapolations to yield 
the total systematic error on the integral. 



C. Integrals 

The Ellis- Jaffe (Eq. ^) and Bjorken (Eq. |20|) sum rules involve integrals over all values 
of x at a fixed Q^. The experimental results do not cov6r clII X cit cLiiy sin gle In the 
measured region of x we must either interpolate or extrapolate our results from the measured 

to some fixed Ql- In the regions of x above and below the measured region, we use 
model-dependent extrapolations. 



1. Measured Region 
Several methods have been used to determine gi at fixed Qo- 

• 1) Assume gi/Fi is independent of and determine gi from gi{x,Ql) = gi/Fi x 
Fi{x,Ql). 

• 2) Assume Ai is independent of and determine gi from [gi{x,Ql) = Ai{x) x 
F,{x,Ql)+^',g^'^{x,Q% 

• 3) Fit the data to a functional form which has semi-empirical dependencies on x and 

such as the fits described above. 

• 4) Do a pQCD fit to determine the quark and gluon distributions and then calculate 
the change in gi going from the measured to the desired kinematics. 

In this paper we will pursue the first two options with emphasis on the first. We note 
that the pQCD fits indicated in Fig. |13] and |l^ show little dependence (compared to 



the errors of the experiments) for gi/Fi at x < 0.5 in the relevant range. For x > 0.6, 
theoretical papers often use approximations in defining the relationship between F2 and Fi 
and sometimes use pQCD fits to F2 instead of empirical fits to the data. At x = 0.75 
typical pQCD fits [ ^ , |117|jTT8|| show gi/Ff differing by 30-50% between the measured 



~ 9.5 (GeV/c)^ and = 3 (GeV/c)^, but these predictions are questionable due to the 
assumptions used. 

Tables pCX[-pCXn| list gi as a function of x at fixed values of 2, 3, and 5 (GeV/c)^ 



for proton, deuteron, and neutron. These results were evaluated by method 1 {gi/Fi 
independent of Q^). Figure O shows the corresponding method 1 results for gi at 
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q2=3 (GeV/c)2. Results for /o°o^3 c/i(x)dx using methods 1 and 2 (Ai independent of Q^, 
g^^) at the same three values of are given in table pCXIII| . Method 1 yields slightly 



larger results in magnitude than method 2, but the difference is smaller than the total error 
for all targets at fixed > 3 (GeV/c)^. The components of the systematic error on the 
integral are shown in Table PCXIV| . The correlated systematic errors due to beam and target 
polarization and F2 dominate. The radiative correction errors tend to be anti-correlated 
between low and high x and thus partly cancel in the integral. 



2. Low X 



The evaluation of Jq '^^ gi{x,Ql)dx can be done by several methods. 



1) Using Regge trajectory- type behavior, gi = x°'g1 at fixed and low x. The 
difference gi-gi has isospin 1 and only one Regge trajectory contributes. The value 
of a is in the general range 0.5 > a > |p.l9|J120|J121|| . For the individual proton and 



deuteron targets there may be more than one pole contributing . There also may 
be even more complex behavior of the singlet term. 

2) Using SMC data from 0.003 < x < 0.03 and Regge extrapolations below x = 0.003. 



• 3) Using the form gi oc ln(l/a;) ||121| , |122(] . 

• 4) Using the parameterization II from Table ^ (the form is Regge inspired at low x). 

• 5) Using pQCD fits. 

The Regge method requires a choice of x range to determine the pole parameter and 
a choice of other possible Regge trajectories. In addition, if gi has Regge behavior at a 
given Q^, it will not have Regge behavior at other since gi evolves with differently at 
different values of x. Table pCXV| shows the results of various options, including using the 
Regge form at Q^=l and 3 (GeV/c)^. Constraining a = gives good fits at both Q'^=l and 
3 (GeV/c)^. However, requiring a = 0.5 gives a rather poor fit {x^/df ~ 2) for the proton. 
We take the average of the four fits in Table PCXV| with a = as the central value of the 
low X extrapolation. The error encompasses all the other models indicated. These averages 
are shown in Table |XXVI| . The values of the integral for proton, deuteron, neutron, and 
proton-neutron (p-n) may not add up exactly due to the non-linearity of the fits. 

"^■^ at low X. 



Recent results from SLAG [11 



indicate that g^ may be behaving as ~ x' 
If the proton behaves in a similar way, then the above extrapolations would be open to 
question. 



3. High x 

The extrapolation to high x was done by two methods: 1) assuming (71 oc (1 — x)^ 
123| | and fitting to the four highest x bins; 2) assuming Ai = 0.75(19 — 16F^/F2)/15 and 



A"^ = 0.75(2 - 3F^ / F^) / {5F^ / Fi) g^. For both gf and g^, both methods gave almost 
identical values of Jq ^ gi{x)dx. For the deuteron the value is 0.000 and for the proton 0.001. 
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We assign an error of ±0.001 in both cases. For the neutron the value of the integral is 
0.001 ± 0.001 where rounding errors account for inexact match with proton and deuteron 
results. The small value of the integral is mostly due to the small value of Fi and not the 



properties of the individual models. The average values are shown in Table PCXVI 



4- Total Integral 



Table PCXVI] shows Fi, the total integral from a; = to 1, in the last column for proton, 
deuteron, neutron, and the difference proton-neutron. The experimentally measured portion 
of the integral makes the largest contribution to 5Ti with the low x extrapolation error a 
close second. The correlation between the measurement errors at low x and the extrapolation 
errors is small compared to the model dependence of the extrapolation. 

The integrals from this experiment, E142 [l^, and SMC [[7tP,|l4| are compared in 



Table pCXVll] at values reported by the other experiments. Comparisons are made 



for the full range in a; of to 1, as well as for the common measured x range between 
experiments. For each comparison, we evaluated the integral using the same assumptions 
about dependence that the other experimenters used. (SMC results were calculated 
by us from their tables). In the experimental range there is good agreement between this 



experiment and the other results. Fi for each experiment, as shown in Table PCXVI1| , are also 
in excellent agreement with the caveat that different x ranges were measured and different 
extrapolations used in the unmeasured region. 

5. Ellis-Jaffe Sum Rule 

Ellis- Jaffe sum rule predictions for the integrals (Eq. ^1]) are shown at three values of in 
Table PCXVII1| . We use a,(M^) = 0.118 ±0.003 H, three active flavors, and 3F - D = 0.58 



with uncertainties of either 0.03 (small) or 0.12 (large) |]2D|. This larger error, as 



discussed in Section [II B 1| , is likely to be an overestimate. For the Ellis-Jaffe sum rule, the 
values for both the "invariant" and Q^-dependent singlet pQCD corrections (see Eq. are 
given. They differ by an amount which is larger than the theoretical error due to as- In the 
case of the deuteron, the experimental errors are comparable to the theoretical difference. 
The measured values of F^* and F^ are shown in the table along with the derived value of F". 
Using the small errors on 3F — D and the "invariant" singlet term the Ellis-Jaffe sum rules 
are violated by 0.023 ±0.007 (deuteron) and 0.032 ±0.012 (proton). A violation implies that 
there could be a significant SU(3) symmetry breaking effect or that there is strange and/or 
gluon spin contributing to the proton spin. If we consider the large 3-F — D errors (larger 
error due to possible symmetry breaking) combined with the Q^-dependent singlet term, the 
deviations from the sum rule reduce to 0.018 ±0.015 (deuteron) and 0.027 ±0.019 (proton). 

6. Bjorken Sum Rule 



The Bjorken sum rule integral (Eq. ^) is given in Table pCXVlll| for three different 
values of Q^. The theoretical value involves only non-singlet pQCD corrections and 
is thus independent of ambiguities associated with the singlet corrections (invariant or 
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Q^-dependent). Theoretical errors depend only on the uncertainty in The measured 
values are from this experiment only and the errors include all correlations. The experimental 
errors are considerably larger than the theoretical errors. Experiment and theory agree 
within one standard deviation. 

A more precise result can be obtained by combining all the experiments [[7|p|,|6|,p!I| which 
published a value of the integrals. SLAG experiment E154 [|11| on the neutron was not 
included because they did not publish an integral of their data alone. We consistently 
used the method with gi/Fi independent of to evolve the results to constant . In all 
experiments, the low-a: extrapolation errors were limited to gi being constant or approaching 
zero as X ^ 0. At Q'^=b {GeV/cf, the combined results are: rf^' = 0.170 ± 0.012, 
r? = 0.130 ± 0.006, and = -0.040 ± 0.008 with a very small xVdf. rf^' is one standard 
deviation from the theoretical value of 0.182 determined from Eq. |2^ with as{Mz) = 0.118. 
The addition of higher twist and other effects described below make the agreement even 
better. If we assume the Bjorken sum rule is true and solve Eq. |2^ for a<j we obtain 
«,(M,) = 0.1231°:°^°. 



D. Quark Polarization 



We used Eqs. p^-|15| to extract our measured value of Qq from the proton and deuteron 
first moments. Then using Eq. ^ (with AG = 0), we extracted the individual polarizations 
of the quarks. It is important to remember that these polarizations have meaning 
within the quark-parton model where ao = AS. In pQCD the interpretation becomes 
scheme-dependent and depends on whether AG{x) contributes to Fi. The results are shown 
in Table PCXIX| and Fig. |18|. Results for both the "invariant" and Q^-dependent pQCD 
singlet coefficients are shown in the table while the figure shows "invariant" results. The 
quantities oq, Au, and Ad are relatively insensitive to the values of F and D, but As is very 
dependent on them. The last two columns of Table X-XIX show the errors on As with two 
different estimates on the errors on 3F — D. As seen in Fig. |18| the results from the deuteron 
and proton targets are consistent with each other (there is only a small correlation between 
the errors). The differences between the "invariant" and Q^-dependent results are smaller 
than the present experimental errors. The negative polarization of the strange quark sea 
of about —0.08 is very significant only if the smaller estimates of F and D are used. Our 
averaged proton and deuteron results for Oq™ = 0.33 ± 0.06, while the world average yields 
Oq ~ 0.31 ± 0.04. The results for oq are significantly smaller than the naive parton model 
prediction of AS = 1, the relativistic parton model prediction of 0.75, the Ellis- Jaffe sum 



rule prediction of 0.58, and a Quenched Lattice Calculation of 0.60 ± 0.05 ||124|| . With such 
a low value of oq, angular momentum conservation (see Eq. |^) requires that the nucleon spin 
be dominated by a combination of gluon polarization and orbital angular momentum, or a 
large charm polarization not included in the formalism above (Eqs. p!^-pO|). There have been 
several approaches to understanding the low value of Qq. These are described in Section II. 
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E. A2 and g2 



Tables pCXXj - |XXXll| show A2 and (72 for each target and spectrometer (and the derived 
neutron) with beam energy of 29.1 GeV. Figures p!9| - |20| show A2 and xg2 for each target and 

m 



are 



spectrometer. The systematic errors are indicated by the bands. Also shown in Fig 
the results from SMC [p|,|l4| for the deuteron and proton and from E142 ]10|| and E154 [T^ for 
the neutron. All results are shown at their measured Q^. There is good agreement between 
the various experiments. As seen in Fig. |1^, the values of A2 are consistent with zero, while 
the values of A2 deviate significantly from zero for x > 0.1. For this experiment the average 
value for all x of A2 is 0.031 ± 0.007 for the proton, 0.003 ± 0.013 for the deuteron, and 
—0.03 ± 0.03 for the neutron. The measured A2 obeys the a/R bound within errors, and at 



almost all kinematics the absolute value of the measured values are significantly lower than 



P5| using gi evaluated from a 
curves for the 7° and the 4.5° 



the bound. The dashed curve is a calculation of from Eq. 
fit to world data discussed in Sec. |V A| of this paper. The g^^ 
kinematics are indistinguishable on the figure. The other theoretical curves are bag model 
predictions | pB|^ which include twist-2 and twist-3 contributions for = 5 (GeV/c)^. At 
high X the E143 results for ^ff indicate a negative trend consistent with the expectations for 
g^^ with a of 43 for 48 degrees of freedom. However, the results are also consistent with 
(^2=0 with a of 52. The deuteron and neutron results are less conclusive because of the 



larger errors and are also consistent with both g^^ and g2 
discussed below along with the moments of gi. 



0. The moments of g2 will be 



F. Higher Moments of gi and §2 



Using our results for both gi and g2, we have computed the third moment of the OPE 
sum rules (Eq. ^Bp, and solved for the twist-3 matrix element d2 and the twist-2 matrix 
element 02. For the measured region 0.03 < a; < 0.8, we evaluated gi, corrected the twist-2 
part of g2 to fixed Q^=5 (GeV/c)^ assuming gi/Fi is independent of Q^, and have averaged 
the two spectrometer results. Possible dependence of ^ was neglected. We neglect 
the contribution from the region < a; < 0.03 because of the suppression factor. For 
0.8 < X < 1, we assume that both gi and g2 behave as (1 — x)^ since at high x, g2 ~ —gi from 
Eq. ^ and Fi 0, and we fit the data for x > 0.56. The uncertainty in the extrapolated 
contribution is taken to be the same as the contribution itself. The results are shown in 



Table PCXXIII| . Our extracted values for d2 are consistent with zero, but the errors are large. 
For comparison, in Table pCXXlV| we quote theoretical predictions [^,^,^,0,0 for and 



dg- For ^2 the proton and neutron results were averaged and a deuteron D-state correction 
was applied. We note that the results for (if and (ig differ in sign from the theoretical QCD 
sum rule calculations [0,^^,^ . The bag model predictions ISBJS^ , however, are of the same 



sign as the data. Ali, Braun and Hiller ||125|| showed that g2 obeys an evolution equation in 
the limit that 00. However, this program of calculation has not been carried out yet. 

To test the Burkhardt-Cottingham sum rule, Eq. 28, we have evaluated the integrals 
Jo,03 g^ix)dx = -0.014 ± 0.028 and Jo,o3 9iix)dx = -0.034 ± 0.082 using the same high-x 
extrapolation as discussed above. These results are consistent with zero. To evaluate the 
integral for x < 0.03 is theoretically challenging. A double logarithmic approximation has 
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been used ||126|| to calculate both gi and g2 in the low-x region yielding g2 oc x '^•^^ at a fixed 
g2 of a few (GeV/c)2. Then J^'^^ g2{x)dx is negligible and the sum rule is confirmed. 

G. Higher Twist Effects 

We have compared our experimental integrals with theoretical predictions using pQCD 
for the finite corrections to various sum rules, which were originally derived at infinite Q^. 
At low it is possible that higher twist effects could also infiuence the evolution of gi. These 
terms generate a multiplicative term of the form {1 + C/[Q'^{1 — x)]} | |127| ]. When going to 
very high order in pQCD there is a confusion between resummation effects generating 
terms and the higher twist terms. There have been several calculations of the corrections 



to Fi using QCD sum rules and the bag model (see the reviews ||128|| ). These take the form 
of an additive correction to the sum rule of the form fi^/Q"^ where t=p, n, or d for proton, 
neutron, or deuteron. From QCD sum rules [^,0] the higher twist contribution to Fi is: 



= ^{a. + U2 + 4/,) + O(^) = ^ + 0{^) (57) 



where a2 (twist-2) and d2 (twist-3) have been calculated from data above and /2 is twist-4. 
The contribution to fi from 02 and d2 is 0.004 for the proton and 0.002 for the deuteron, 
which are quite negligible at our average = 3 (GeV/c)^ and small compared to the 
estimated contributions from /2. Table pCXXIV| shows calculated values of fi2 and /^g using 



bag models and QCD sum rules. The sum rule calculations average about —0.02 for the 
proton and —0.013 for the deuteron and thus would have an effect on the calculation of 
Fi at our average of 3 (GeV/c)^ comparable to our experimental error. The bag model 
calculations are similar in magnitude but opposite in sign to the sum rule calculations. A 
different type of calculation, using a diquark model ||130|| gives a higher twist contribution 



of a different form than Eq. ^ which numerically is 1% or less of both Ti'' and F^* for 
> 2 (GeV/c)^. Using data from this experiment, Ji and Melnitchouk ||8^ have extracted 



values for the twist-4 matrix element f2- Combining this with results for 02 and d2 they find 
for the proton = 0.04 ± 0.02, and for the neutron /ig = 0.03 ± 0.04. 

The QCD sum rule higher-twist correction and a Fade summation of the perturbative 
terms have been applied to the Bjorken sum rule (Eq. ^) by Ellis et al. | 129 |. They then 



use world data, including the preliminary results from this experiment, and find excellent 
agreement between experiment and theory. Working backwards, they determine the best 
value of Us is 0.117to;Q07 ± 0.002 where the first set of errors is experimental and the second 
theoretical, in excellent agreement with the world average of 0.118 ± 0.003. 

H. Pion Asymmetry 

The asymmetries for 7r+ and vr" for our primary energy of 29 GeV corresponding to target 



polarization parallel and anti-parallel to the beam direction (Ap are shown in Figs. 21-22. 
These data were measured using our polarized NH3 and ND3 targets, and were corrected for 
beam and target polarizations as well as dilution to obtain the asymmetries from polarized 
protons and deuterons. The asymmetries are small, but for the proton may be slightly 
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positive for both signs of pions at small momentum. Table PCXXV] gives the pion asymmetry 
at the beam energy of 29.1 GeV for proton and deuteron targets. Data for other beam 
energies and for the target spins oriented perpendicular to the beam direction have much 
larger statistical errors and are consistent with zero. 



I. Resonance Region 

Results for Ai + riA2 extracted via the asymmetry method for the resonance region 
{W^ < 5 GeV^) are shown in Table PCXXV1| and in Fig. |3[ Also shown in Fig. |2| is the 



data of Baum at al. [132|. The two data sets agree within the errors of both measurements. 
The asymmetry is negative and close to the expected value Ai = for the A resonance. 
In the region of the D13 and the Su resonances (W"^ ~ 2.34 GeV^) Ai + 77^2 is large and 
positive. Although 77 is large for our kinematics, a small value of A2 would imply that 
Ai+r]A2 is ^ Al. 



Table pCXXVlI lists the systematic errors on Ai + by category. The procedure for 



estimating these is as follows: a) radiative corrections (Arc)'- the maximum deviation in 
the radiative correction resulting from a 50% change in the input asymmetries; b) model 
dependence (model): the worst-case change due to using various cross-section models in 
the extraction of Ai + riA2 from A\\; c) central angle {6): uncertainty due to the location 
of the central angle of the spectrometers; d) energy calibration {E'): uncertainty due to 
the spectrometer energy calibration; e) spectrometer resolution (Resol): the maximum 
difference obtained by varying the width of the hodoscope fingers by 20% and re-running the 
Monte Carlo routine; f) polarization (Pf,Pf): combined uncertainty in the beam and target 
polarizations; g) dilution factor (/): uncertainty based on the variations in the calculated 
dilution factor with various cross-section models and the stated uncertainty in the target 
composition; h) R{x,Q^) (R): uncertainty arising from lack of knowledge of i? = aL/ax', i) 
no transverse data ((72): uncertainty in the extraction of gi from Ai + riA2 due to the lack 
of knowledge about A2. This was estimated using the maximum deviations in gi assuming 
A2 = ^ and g2 = 0. Even if A2 were as large as 0.3, the extracted values of gi would shift 
by less than 0.014, which is small compared to the statistical errors on each point. By far 
the largest error comes from radiative and resolution corrections. 

In addition to Ai + 77^2, Table pCXXVl| also shows the results for gi for the resonance 



region. Figure ^ shows gi for proton and deuteron (per nucleon) measured with the two 
spectrometers as a function of W"^. The data of Baum et al. [ |132|] are taken at similar 
kinematics and converted to gi for comparison by assuming ^2 = 0. Within errors, the two 
measurements agree well. Both data sets show a negative contribution in the region of the 
A(1232) resonance at GeV^, and a strongly positive contribution just above = 2 

GeV^ where the Su and D13 resonances are important. This peak is less pronounced for 
the deuteron. The solid lines show the Monte Carlo simulation. 

Figure ^ shows the integrals Tii^Q"^) for proton and neutron, evaluated at the average 



for the resonance region (M^ < < 4 GeV^). We summed our resonance results directly 
(where does not vary much) and then added a contribution from smaller x (larger W"^) 
at the same fixed by interpolating the 9.7 and 16 GeV data to the appropriate Q^. The 
neutron integrals were derived assuming a 5% D-state probability for the deuteron. The 
statistical errors assigned to the integral over the deep- inelastic region (F^^^) correspond to 



46 



the weighted average of the statistical errors on the corresponding 9.7 and 16.2 GeV data 
points used in the interpolation. Systematic errors on the total integral r*[°* were calculated 
using the systematic uncertainties for the measured gi in the resonance region added linearly 
to the systematic errors for the deep-inelastic region, which are highly correlated with each 
other. For the x < 0.03 extrapolation we simply took the overall parameterization of the 
data and integrated it from x = to a; = 0.03. Extrapolation errors for the region below 
the last measured datum at x = 0.03 were taken to be as large as the values themselves. 
Table pCXX Villi lists for each target the numerical values for the integrals in the resonance 
region alone (r^*^*^), in the deep-inelastic region (Ff^^), for the low-x extrapolation (F|^*), 
and for the combined total (F*[°*). 

Although several models for the evolution of Ti{Q^) exist |5||T3|,|T3|,|T35|,|T3§ , we 



show here only two representative ones, together with the evolution |]33| of the world's 
deep-inelastic data due to the changing coupling constant as- Although the GDH sum rule 
is strictly valid only at = where Fi(Q^) vanishes, it can be used to predict the slope 
of Fi(Q^) for small Q^. The solid line at low shows Fi = — k^Q^/SM^ in which k is 



the anomalous magnetic moment of either the proton or neutron. Burkert and loffe ||135 
consider the contributions from the resonances using the code AO, and the nonresonant 
contributions using a simple higher-twist-type form fitted to the deep-inelastic data. Their 
model is constrained to fit both the GDH and the deep-inelastic limits, and it describes 
the data quite well. Soffer and Teryaev [ |136| | assume that the integral over gi + g2 varies 
smoothly from high where ~ down to = 0. Using their simple prediction for this 
integral and subtracting the contribution from (72 using the Burkhardt-Cottingham sum rule 
||69| gives the dashed curves in Fig. |2^, which also agree quite well with our data. 

The present spin structure function data in the region of the nucleon resonances allow 
us to determine the integrals Fi(Q^) for the first time at below 2 (GeV/c)^. In contrast 
to the nearly flat behavior in the deep-inelastic region above = 2 (GeV/c)^, Fi varies 
rapidly below = 2 (GeV/c)^. Models that interpolate between the deep-inelastic and 
GDH limits describe the data quite well in this non-perturbative regime. 



VI. CONCLUSION 

In summary, we have presented final results from SLAG Experiment E143 on the spin 
structure functions gi and g2 for proton and deuteron targets covering a wide range of 
kinematics from the deep-inelastic to the resonance region. For deep-inelastic data the ratio 
gi/Fi is consistent with being independent of for > I (GeV/c)^, but also consistent 
with pQGD NLO fits which show a weak dependence. We have evaluated the first 
moments of gi, using a Regge form for the unmeasured low x region. The Ellis- Jaffe sum 
rules are a function of the SU(3) parameters F and D and the validity of the sum rules 
depend critically on the errors assigned to these parameters. We find = 0.33 ± 0.06, 
and in the parton model interpretation we find the average results: Au = 0.84 ± 0.02, 
Ad = -0.42 ± 0.02, and As = -0.09 ± 0.02 or ±0.05 depending on 6F/D. Gombined 
world data are consistent with the Bjorken sum rule at the one standard deviation level 
of 7%. Results for the twist-3 matrix element extracted from the higher moments of gi 
and g2 are consistent with calculations within the large errors. The resonance region data 
show the theoretically expected asymmetries at the A (1232) peak and larger than expected 
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asymmetries (at least at low Q'^) in the region of the 5*11 and D13 resonances. The first 
moment of gfiQ"^) decreases with decreasing at low toward the GDH sum rule limit 
as predicted by several models. The asymmetry of pions is close to zero. 

With the current round of experiments we now have good knowledge of the distribution 
of quark spins for x > 0.003. A complete understanding of the spin structure of the nucleon 
awaits experiments to measure directly the gluon spin distribution and to probe the quark 
spin distribution at lower x. 
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TABLE I. Quark helicity predictions from the nonrelativistic quark-parton model (NR QPM) 
where AG = and from the relativistic quark-parton model (R QPM) [^,22] with AG = and 
AG(Q2 = 1 (GeV/c)2) = 1.6 ±0.9^. 



NR QPM 



R QPM 



R QPM + gluons 



asAG/2TT 
Au - asAG/2Tr 
Ad - asAG/2Tr 
As - asAG/27r 
ao = AS - 3asAG/2TT 





1.33 
-0.33 


1.0 





1.0 

-0.25 


0.75 



0.13 ±0.08 
0.87 ±0.08 
-0.38 ±0.08 
-0.13 ±0.08 
0.36 ±0.24 



TABLE IL Beam heating correction results at beam intensity of 4 x 10^ electrons/pulse. 



Target NH3 NH3 ND3 ND3 

Polarization Long. Tran. Long. Tran. 

Pi (%) 75 ± 1.9 75 ± 1.9 30 ± 1.2 30 ± 1.2 

Pm {%) 68.3 ± 1.7 68.5 ± 1.7 24.9 ± 1.0 24.8 ± 1.0 

z 0.924 ±0.029 0.903 ± 0.033 0.912 ± 0.023 0.931 ± 0.021 

Pt {%) 67.7 ± 1.7 67.8 ± 1.7 24.4 ± 1.0 24.4 ± 1.0 

Cheat 0.0081 ± 0.0036 0.0103 ± 0.0040 0.0197 ± 0.0064 0.0157 ± 0.0053 
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TABLE III. Proton Aii and A± results with statistical errors for E=29.1 GeV at the measured 



in (GeV/c)^. Also shown are the radiative corrections A'j-c and A^jr^. which were applied to the 
data. 



X 


<Q2> 


^11 






A± 


rc 


e = 4.5° 


0.028 


1.17 


-0.026ib 0.054 




0.014 


0.031± 0.063 


0.004 


0.031 


1.27 


0.048± 0.026 




0.014 


O.OlOib 0.032 


0.004 


0.035 


1.40 


0.091± 0.019 




0.013 


0.012± 0.024 


0.004 


0.039 


1.52 


0.060± 0.016 




0.012 


0.007± 0.020 


0.004 


0.044 


1.65 


0.076± 0.015 




0.011 


0.004± 0.018 


0.004 


0.049 


1.78 


0.083± 0.014 




0.010 


0.008ib 0.017 


0.004 


0.056 


1.92 


0.082± 0.013 




0.009 


0.003± 0.016 


0.004 


0.063 


2.07 


0.082± 0.012 




0.008 


0.014± 0.015 


0.004 


0.071 


2.22 


0.086± 0.011 




0.007 


0.012± 0.014 


0.004 


0.079 


2.38 


0.102± 0.012 




0.006 


-0.009± 0.014 


0.004 


0.090 


2.53 


0.081± 0.012 




0.005 


-0.006± 0.014 


0.004 


0.101 


2.69 


0.114± 0.012 




0.004 


0.004± 0.014 


0.004 


0.113 


2.84 


0.108± 0.013 




0.004 


-0.012± 0.015 


0.004 


0.128 


3.00 


0.097± 0.013 




0.003 


O.OlOib 0.015 


0.004 


0.144 


3.15 


0.086± 0.013 




0.003 


-0.031± 0.015 


0.004 


0.162 


3.30 


0.113± 0.013 




0.003 


0.023± 0.016 


0.004 


0.182 


3.45 


O.llOib 0.014 




0.002 


-0.021± 0.016 


0.004 


0.205 


3.59 


0.097± 0.014 




0.002 


0.043± 0.017 


0.004 


0.230 


3.73 


0.118± 0.015 




0.002 


-0.005± 0.018 


0.004 


0.259 


3.85 


0.107± 0.015 




0.002 


0.017± 0.019 


0.004 


0.292 


3.98 


0.096± 0.016 




0.002 


-0.055± 0.020 


0.004 


0.329 


4.09 


O.llOi 0.018 




0.002 


-0.005± 0.022 


0.004 


0.370 


4.20 


0.080± 0.020 




0.002 


0.012± 0.024 


0.003 


0.416 


4.30 


0.140± 0.023 




0.002 


0.002± 0.028 


0.003 


0.468 


4.40 


0.140± 0.026 




0.002 


-0.048± 0.032 


0.003 


0.526 


4.47 


0.134± 0.031 




0.002 


-0.037± 0.038 


0.002 


0.592 


4.55 


0.066± 0.037 




0.003 


-0.029± 0.045 


0.002 


0.666 


4.63 


0.075± 0.045 




0.000 


-0.013± 0.055 


0.002 


0.749 


4.70 


0.128± 0.062 




-0.007 


-0.114zb 0.074 


0.004 


e = 7.0° 


0.071 


2.91 


0.261± 0.095 




0.018 


-0.122± 0.102 


0.006 


0.079 


3.17 


0.159± 0.043 




0.015 


-0.002± 0.049 


0.006 


0.090 


3.48 


0.115± 0.029 




0.012 


0.040± 0.034 


0.006 


0.101 


3.79 


0.143± 0.024 




0.010 


-0.016± 0.027 


0.006 


0.113 


4.11 


0.158± 0.022 




0.008 


O.Ollzb 0.025 


0.006 


0.128 


4.43 


0.164± 0.021 




0.006 


-0.026± 0.023 


0.006 


0.144 


4.78 


0.159± 0.020 




0.005 


0.007± 0.022 


0.006 


0.162 


5.13 


0.171± 0.019 




0.004 


0.018± 0.021 


0.006 


0.182 


5.49 


0.192± 0.019 




0.003 


0.052± 0.021 


0.006 


0.205 


5.86 


0.215± 0.019 




0.003 


0.012± 0.021 


0.006 
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0.230 


6.24 


0.150± 0.019 


0.002 


0.007± 0.022 


0.006 


0.259 


6.60 


0.250± 0.020 


0.002 


O.OOlib 0.022 


0.006 


0.292 


6.97 


0.197± 0.021 


0.002 


-0.023± 0.023 


0.005 


0.329 


7.34 


0.195± 0.022 


0.002 


-0.012± 0.025 


0.005 


0.370 


7.69 


0.190± 0.024 


0.002 


-0.018± 0.027 


0.005 


0.416 


8.04 


0.244± 0.026 


0.002 


-0.009± 0.030 


0.004 


0.468 


8.37 


0.223± 0.030 


0.002 


0.005± 0.034 


0.004 


0.526 


8.68 


0.233± 0.034 


0.002 


-0.015± 0.039 


0.003 


0.592 


8.99 


0.224± 0.041 


0.002 


-O.Olli 0.047 


0.003 


0.666 


9.26 


0.155± 0.051 


0.002 


-0.093± 0.063 


0.002 


0.749 


9.53 


0.223± 0.069 


0.005 


-0.097± 0.094 


0.001 
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TABLE IV. Proton results with statistical errors for E=16.2 GeV at the measured in 



(GeV/c)^. Also shown are the radiative corrections Arc which were applied to the data. 



X 


<Q2> 




A\\ 


4" 


e = 4.5° 


(J.Uzz 


0.47 




A AOO 1 A AO /I 

0.023± 0.024 


A AA r" 

0.005 


n no A 

L).L)z4 


0.51 




A riCO 1 A A1 A 

0.0d6± 0.019 


A AA r* 

0.005 


0.027 


0.55 




A ACO 1 A A1 O 

0.058± 0.018 


A AAC 

0.005 


0.031 


A r" A 

0.59 




A Ar" A 1 A A1 O 

0.054± 0.016 


A AAr" 

0.005 


0.035 


0.64 




A AO/^ 1 A A1 f 

0.036ib 0.015 


0.005 


0.039 


0.68 




0.041± 0.014 


A AA V 

0.005 


0.044 


0.73 




AATi 1 AAIO 

0.051zb 0.013 


0.005 


0.049 


A TO 

0.78 




A C^fl A 1 A A1 C\ 

0.074± 0.012 


A AAr" 

0.005 


0.056 


A OO 

0.83 




A Ar'O 1 A A1 C\ 

0.052± 0.012 


A AA A 

0.004 


0.0d3 


A O O 

0.88 




A A/^O 1 A A1 O 

0.0d8± 0.012 


A AA A 

0.004 


0.071 


A AO 

0.92 




AA/10 1 AA10 

0.043± 0.012 


A AA A 

0.004 


0.079 


A A'7 

0.97 




C\ C\ A C 1 AA10 

0.046± 0.012 


A AAO 

0.003 


0.090 


1 A1 

1.01 




A C\f*C\ 1 A A1 C\ 

0.0d2± 0.012 


A AAO 

0.003 


0.101 


1 A/^ 

1.06 




A C\T\ 1 A A1 C\ 

0.071± 0.012 


A AAO 

0.003 


0.113 


1.10 




f\ f\ A f\ 1 AA"tO 

0.049± 0.012 


0.003 


0.128 


1.14 




A A/^A 1 A A 1 O 

0.060ib 0.012 


0.002 


0.144 


1.18 




A C\r O 1 A A10 

0.056± 0.012 


A AAO 

0.002 


0.162 


1.22 




0.068± 0.013 


0.002 


0.182 


1.26 




C\ C\f -\ 1 AA10 

0.051± 0.013 


0.002 


0.205 


1.29 




A f\r*f\ 1 A A1 o 

0.060ib 0.013 


0.002 


0.230 


1.32 




0.047± 0.014 


0.002 


U.zoy 


i.OO 




U.U4izt U.U14 


U.UU2 


0.292 


1.38 




0.055± 0.015 


0.002 


0.329 


1.40 




0.043± 0.016 


0.003 


0.370 


1.43 




0.055± 0.017 


0.003 


0.416 


1.45 




0.079± 0.019 


0.000 


0.468 


1.46 




0.094± 0.021 


0.003 


0.527 


1.48 




0.082± 0.024 


0.010 


0.593 


1.49 




0.089± 0.028 


0.004 


0.668 


1.51 




0.052± 0.031 


-0.009 


0.752 


1.52 




0.258± 0.449 


-0.029 


e = 7.0° 


0.044 


0.98 




0.072± 0.390 


0.012 


0.049 


1.06 




0.055± 0.078 


0.012 


0.056 


1.16 




0.091± 0.043 


0.011 


0.063 


1.26 




0.034± 0.031 


0.010 


0.071 


1.37 




0.092± 0.026 


0.009 


0.079 


1.47 




0.082± 0.023 


0.008 


0.090 


1.58 




0.113± 0.022 


0.007 


0.101 


1.69 




O.lOlib 0.021 


0.006 
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0.113 


1.80 


O.lOSit 


0.019 


0.005 


0.128 


1.91 


O.llSib 


0.018 


0.005 


0.144 


2.03 


0.120± 


0.018 


0.004 


0.162 


2.14 


0.103± 


0.017 


0.003 


0.182 


2.26 


0.105± 


0.018 


0.003 


0.205 


2.36 


0.116± 


0.018 


0.003 


0.230 


2.47 


0.136± 


0.018 


0.002 


0.259 


2.57 


0.165± 


0.019 


0.002 


0.292 


2.67 


0.159± 


0.019 


0.002 


0.329 


2.76 


0.127± 


0.020 


0.002 


0.370 


2.85 


0.157± 


0.022 


0.002 


0.416 


2.94 


0.142± 


0.023 


0.003 


0.468 


3.02 


O.lSOib 


0.025 


0.003 


0.527 


3.08 


0.127± 


0.029 


0.004 


0.593 


3.15 


0.070± 


0.032 


-0.001 


0.668 


3.21 


0.136± 


0.038 


0.003 


0.752 


3.27 


-0.280± 


0.242 


0.009 
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TABLE V. Proton Ay results with statistical errors for E= 9.7 GeV at the measured in 



(GeV/c)^. Also shown are the radiative corrections Arc which were applied to the data. 



X 


<Q2> 




A\\ 


4" 


e = 4.5° 


n r\oo 

U.Uzo 


A O O 

0.28 




A A01 1 A r\ A C 

0.021± 0.046 


A AAO 

0.003 


n no 1 


A OA 

0.30 




A A01 1 A AOO 

0.021± 0.022 


A AAO 

0.002 


U.Uoo 


A O 1 

0.31 




C\ f\ AO 1 AA1A 

0.048± 0.019 


A AAO 

0.002 


L).Uo9 


A OO 

0.33 




A A1 1 1 A A1 ^7 

0.011± 0.017 


A AAO 

0.002 


A A /I 

(J.L)44 


A O C 

0.35 




0.046± 0.016 


A AA 1 

0.001 


A A /I A 

(J. 049 


A O/^ 

0.36 




AA/IO 1 AA1/^ 

0.043± 0.016 


A AA 1 

0.001 


0.056 


0.38 




AA/1A 1 AAIP* 

0.040zb 0.015 


0.001 


A A/^ O 

0.063 


A Af\ 

0.40 




A AO/^ 1 A A1 r 

0.026± 0.015 


A AA1 

0.001 


A ATI 

0.071 


A y< 1 

0.41 




A AOA 1 A A1 C 

0.030± 0.015 


A AA1 

0.001 


A AO A 

0.080 


A /I O 

0.43 




A AOO 1 A A1 CT 

0.032± 0.015 


A AA 1 

0.001 


A AAA 

0.090 


0.44 




AA10 1 AA1/f 

0.018± 0.014 


A AAA 

0.000 


A 1 A 1 

0.101 


A A^ 

0.45 




A AO A 1 A A 1 /I 

0.024± 0.014 


A AAA 

0.000 


A 1 1 O 

O.llo 


A AIT 

0.47 




AAy<1 1 AA10 

0.041± 0.013 


A AAA 

0.000 


A 1 OO 


A AO 

0.48 




A AAO 1 A A1 O 

0.002± 0.013 


A AAA 

0.000 


0.144 


0.49 




A AAO 1 A A "t O 

0.002ib 0.013 


0.000 


0.162 


A r A 

0.50 




AA10 1 AA10 

0.018± 0.013 


A AAA 

0.000 


A 1 OO 

0.182 


ATI 

0.51 




AAO/1 1 AA10 

0.024± 0.013 


A AAO 

—0.002 


A OAr" 

0.205 


0.52 




A A1 f 1 A A1 O 

0.015± 0.013 


A AAO 

—0.008 


A OO 1 

0.231 


A r" o 

0.53 




A Ar" A 1 A A1 A 

0.054± 0.014 


A AA A 

—0.004 


0.259 


0.53 




0.047± 0.014 


0.015 


U.zyZ 


U.o4 




U.Uoiih U.Uiu 


U.Ui4 


0.329 


0.55 




0.020± 0.018 


0.001 


0.370 


0.55 




-0.002± 0.019 


-0.007 


0.417 


0.56 




-0.004± 0.021 


-0.015 


0.469 


0.56 




-0.034± 0.021 


-0.020 


0.527 


0.57 




0.007± 0.029 


-0.034 


0.594 


0.57 




-0.013± 0.064 


-0.043 


0.669 


0.57 




O.OlOi 0.179 


-0.039 


0.753 


0.58 




0.021± 0.214 


-0.036 


0.847 


0.58 




0.031± 0.257 


-0.032 


e = 7.0° 


0.063 


0.60 




-0.013± 0.093 


0.003 


0.071 


0.64 




0.031± 0.038 


0.003 


0.080 


0.69 




0.076± 0.025 


0.003 


0.090 


0.74 




0.064± 0.019 


0.003 


0.101 


0.78 




0.057± 0.016 


0.003 


0.113 


0.82 




0.073± 0.015 


0.003 


0.128 


0.86 




0.048± 0.015 


0.003 


0.144 


0.90 




0.047± 0.014 


0.002 


0.162 


0.93 




0.065± 0.013 


0.002 
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0.182 


0.97 


0.087± 


0.013 


0.002 


0.205 


1.00 


0.068± 


0.013 


0.002 


0.231 


1.03 


0.084± 


0.013 


0.002 


0.259 


1.06 


0.070± 


0.013 


0.003 


0.292 


1.09 


0.055± 


0.013 


0.004 


0.329 


1.12 


0.088± 


0.014 


0.005 


0.370 


1.14 


0.088± 


0.014 


-0.004 


0.417 


1.16 


0.082± 


0.016 


0.001 


0.469 


1.18 


0.084± 


0.016 


0.018 


0.527 


1.21 


0.086± 


0.018 


0.010 


0.594 


1.22 


-O.OOlib 


0.018 


-0.004 


0.669 


1.23 


0.022± 


0.018 


-0.026 


0.753 


1.25 


-0.044± 


0.047 


-0.061 


0.847 


1.26 


-0.018± 


0.149 


-0.064 
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TABLE VI. Deuteron A\\ and A_i results with statistical errors for E=29.1 GeV at the measured 



in (GeV/c)^. Also shown are the radiative corrections A'j-c and A^jr^. which were applied to the 
data. 



X 


<Q2> 


^11 






A± 


rc 


e = 4.5° 


0.028 


1.17 


-0.042ib 0.075 




-0.004 


-0.138± 0.175 


0.002 


0.031 


1.27 


0.035± 0.030 




-0.004 


-0.114± 0.146 


0.002 


0.035 


1.40 


0.004± 0.021 




-0.004 


0.060± 0.083 


0.002 


0.039 


1.52 


0.043± 0.019 




-0.004 


0.038± 0.046 


0.002 


0.044 


1.65 


-O.Ollzb 0.017 




-0.004 


0.050± 0.038 


0.002 


0.049 


1.78 


-0.009± 0.016 




-0.004 


0.022± 0.034 


0.002 


0.056 


1.92 


0.012± 0.015 




-0.003 


-0.043± 0.031 


0.002 


0.063 


2.07 


O.OlOib 0.014 




-0.003 


-0.040± 0.029 


0.002 


0.071 


2.22 


0.014± 0.013 




-0.003 


0.009± 0.028 


0.002 


0.079 


2.38 


0.023± 0.013 




-0.003 


-0.009± 0.029 


0.002 


0.090 


2.53 


0.038± 0.014 




-0.003 


0.032± 0.029 


0.002 


0.101 


2.69 


0.028± 0.014 




-0.003 


-0.031± 0.030 


0.002 


0.113 


2.84 


0.037± 0.015 




-0.002 


0.012± 0.031 


0.002 


0.128 


3.00 


0.079± 0.015 




-0.002 


-O.OOlib 0.032 


0.002 


0.144 


3.15 


0.053± 0.016 




-0.002 


-0.023± 0.033 


0.002 


0.162 


3.30 


0.046± 0.016 




-0.001 


0.043± 0.035 


0.002 


0.182 


3.45 


0.054± 0.017 




-0.001 


-0.042± 0.036 


0.002 


0.205 


3.59 


0.049± 0.017 




0.000 


-0.013± 0.038 


0.002 


0.230 


3.73 


0.020± 0.019 




0.000 


0.008± 0.040 


0.002 


0.259 


3.85 


0.021± 0.020 




0.000 


0.095± 0.043 


0.002 


0.292 


3.98 


0.054± 0.021 




0.001 


-0.016± 0.047 


0.002 


0.329 


4.09 


0.078± 0.023 




0.001 


-0.025± 0.052 


0.002 


0.370 


4.20 


0.072± 0.026 




0.001 


-0.026± 0.059 


0.002 


0.416 


4.30 


0.063± 0.030 




0.001 


-0.050± 0.068 


0.001 


0.468 


4.40 


O.OlOzb 0.036 




0.001 


-0.071± 0.081 


0.001 


0.526 


4.47 


0.065± 0.043 




0.001 


0.077± 0.098 


0.001 


0.592 


4.55 


0.057± 0.052 




0.001 


-0.116± 0.120 


0.001 


0.666 


4.62 


0.023± 0.066 




-0.002 


0.160± 0.150 


0.002 


0.749 


4.70 


-0.190zb 0.091 




-0.006 


0.150± 0.203 


0.004 


e = 7.0° 


0.071 


2.91 


0.044± 0.108 




-0.004 


0.088± 0.237 


0.003 


0.079 


3.17 


-0.020± 0.049 




-0.004 


0.081± 0.094 


0.003 


0.090 


3.48 


0.044± 0.033 




-0.004 


-0.024± 0.062 


0.003 


0.101 


3.79 


0.044± 0.027 




-0.004 


0.075± 0.050 


0.003 


0.113 


4.11 


0.023± 0.025 




-0.004 


-0.036± 0.046 


0.003 


0.128 


4.44 


0.053± 0.023 




-0.003 


-0.014± 0.043 


0.003 


0.144 


4.78 


O.llOib 0.023 




-0.003 


-0.038± 0.041 


0.003 


0.162 


5.13 


0.051± 0.022 




-0.003 


-0.027± 0.041 


0.003 


0.182 


5.49 


0.133± 0.022 




-0.002 


0.035± 0.040 


0.003 


0.205 


5.86 


0.067± 0.022 




-0.002 


-0.085± 0.041 


0.003 
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0.230 


6.23 


0.088± 0.023 


-0.001 


-0.056± 0.043 


0.003 


0.259 


6.60 


O.OSlib 0.024 


0.000 


0.065± 0.045 


0.003 


0.292 


6.97 


0.102± 0.026 


0.000 


-0.028± 0.048 


0.003 


0.329 


7.33 


0.108± 0.028 


0.001 


-0.034± 0.052 


0.002 


0.370 


7.69 


0.143± 0.030 


0.001 


-0.024± 0.057 


0.002 


0.416 


8.03 


0.089± 0.034 


0.001 


-0.018± 0.064 


0.002 


0.468 


8.37 


0.125± 0.039 


0.000 


-0.003± 0.073 


0.002 


0.526 


8.67 


0.172± 0.046 


0.000 


-0.023± 0.086 


0.002 


0.592 


8.98 


0.094± 0.056 


-0.001 


0.235± 0.108 


0.002 


0.666 


9.26 


0.086± 0.070 


-0.001 


-0.125± 0.147 


0.002 


0.749 


9.52 


0.193± 0.096 


0.000 


-0.068± 0.212 


0.002 
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TABLE VII. Deuteron results with statistical errors for E=16.2 GeV at the measured in 



(GeV/c)^. Also shown are the radiative corrections Arc which were applied to the data. 



X 


<Q2> 




A\\ 


4" 


e = 4.5° 


(J.Uzz 


0.47 




— O.Oloi 0.064 


n nn/^ 

—0.006 


U.L)z4 


0.51 




r\ nor i n r\Ar\ 

— 0.035± 0.049 


n nn/^ 

—0.006 


0.027 


0.55 




— 0.012± 0.045 


n nnc 

—0.005 


0.031 


0.59 




0.012± 0.041 


n nnr 

—0.005 


0.035 


0.64 




n nnn i n noo 

— 0.009± 0.038 


—0.004 


0.039 


0.68 




n nr/^ i n no/f 

0.056± 0.034 


n nn a 

— 0.004 


0.044 


0.73 




n noi 1 n noo 

— 0.021± 0.028 


—0.004 


0.049 


0.78 




0.046ib 0.026 


n nno 

—0.003 


0.056 


0.83 




n nno i n nor 

0.003± 0.025 


n nno 

—0.003 


0.0d3 


0.87 




n noi 1 n nocr 

0.021± 0.025 


n nno 

— 0.003 


0.071 


0.92 




n no/^ 1 n nocr 

0.026± 0.025 


n nno 

—0.003 


0.079 


0.97 




n ni *^ 1 n nocr 

— 0.01d± 0.025 


n nno 

—0.002 


0.090 


1.01 




n noo 1 n nor 

0.038ib 0.025 


n nno 

—0.002 


0.101 


1.06 




n ni n i n nor 

0.019ib 0.025 


n nno 

—0.002 


0.113 


1.10 




n m r' i n nor 

0.016ib 0.025 


—0.001 


0.128 


1.14 




n n/^o 1 n nor 

0.063ib 0.025 


—0.001 


0.144 


1.18 




n nrn i n nor 

0.050± 0.025 


n nni 

—0.001 


0.162 


1.22 




n nor i n nor 

— 0.025± 0.025 


—0.001 


0.182 


1.25 




n n /( o 1 n no/^ 

0.043ib 0.026 


0.000 


0.205 


1.29 




0.042ib 0.027 


0.000 


U.zoU 


i.OZ 




U.UiDih U.U2/ 


U.UUU 


0.259 


1.35 




0.022± 0.029 


0.001 


0.292 


1.37 




0.062± 0.031 


0.001 


0.329 


1.40 




0.023± 0.033 


0.001 


0.370 


1.42 




-0.031± 0.036 


0.001 


0.416 


1.44 




0.013± 0.040 


-0.001 


0.468 


1.46 




0.089± 0.045 


0.001 


0.527 


1.48 




0.014± 0.053 


0.003 


0.593 


1.49 




0.071± 0.063 


-0.002 


0.668 


1.50 




0.037± 0.076 


-0.013 


e = 7.0° 


0.049 


1.06 




-0.075± 0.201 


-0.005 


0.056 


1.16 




0.062± 0.106 


-0.005 


0.063 


1.26 




-0.003± 0.075 


-0.004 


0.071 


1.36 




0.054± 0.060 


-0.004 


0.079 


1.47 




-0.073± 0.054 


-0.004 


0.090 


1.58 




0.121± 0.049 


-0.004 


0.101 


1.69 




0.054± 0.046 


-0.003 


0.113 


1.80 




0.027± 0.041 


-0.003 


0.128 


1.91 




0.018± 0.038 


-0.002 
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0.144 


2.03 


0.039± 


0.036 


-0.002 


0.162 


2.14 


0.060± 


0.036 


-0.002 


0.182 


2.25 


0.097± 


0.036 


-0.001 


0.205 


2.36 


0.045± 


0.036 


-0.001 


0.230 


2.47 


0.041± 


0.038 


0.000 


0.259 


2.57 


0.041± 


0.039 


0.000 


0.292 


2.67 


0.115± 


0.040 


0.001 


0.329 


2.76 


0.160± 


0.042 


0.001 


0.370 


2.85 


0.039± 


0.045 


0.001 


0.416 


2.93 


0.055± 


0.049 


0.001 


0.468 


3.01 


0.168± 


0.054 


0.001 


0.527 


3.08 


0.139± 


0.062 


0.001 


0.593 


3.15 


0.055± 


0.073 


-0.002 


0.668 


3.21 


-0.013± 


0.088 


-0.001 
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TABLE VIII. Deuteron results with statistical errors for E= 9.7 GeV at the measured in 



(GeV/c)^. Also shown are the radiative corrections Arc which were applied to the data. 



X 


<Q2> 




A\\ 


4" 


e = 4.5° 


n r\oo 

U.Uzo 


A O O 

0.28 




A ArO 1 A A/^A 

0.053± 0.069 


A A A CT 

—0.005 


n no 1 


A OA 

0.30 




A rxfl 1 A AO yl 

— 0.0b7± 0.034 


A A A CT 

—0.005 


U.Uoo 


A O 1 

0.31 




A A1 A 1 A AOO 

— 0.014± 0.028 


A AA A 

—0.004 


U.Uo9 


A OO 

0.33 




A A1 T 1 A AOC 

0.017± 0.025 


A AA A 

—0.004 


0.044 


A O C 

0.35 




A AA1 1 A AO /I 

0.001± 0.024 


A AA A 

— 0.004 


0.050 


A O/^ 

0.36 




A AOO 1 A AOO 

0.022± 0.023 


A AA A 

— 0.004 


0.056 


0.38 




A Al ^ 1 A AOO 

— 0.017± 0.022 


—0.003 


0.063 


A Af\ 

0.40 




A AAO 1 A A01 

— 0.002± 0.021 


A AAO 

—0.003 


C\ ATI 

0.071 


A y< 1 

0.41 




A AOO 1 A AOO 

0.032± 0.022 


A AAO 

—0.003 


f\ AO A 

0.080 


A /I O 

0.43 




A AA A 1 A AOO 

— 0.004± 0.022 


A AAO 

— 0.003 


A AAA 

0.090 


0.44 




A AAO 1 A A01 

0.008± 0.021 


A AAO 

— 0.002 


A 1 A 1 

0.101 


A A^ 

0.45 




A AOO 1 A AOA 

0.038± 0.020 


A AAO 

—0.002 


A 1 1 O 

O.llo 


A AIT 

0.47 




A AAO 1 A A1 A 

— 0.002± 0.019 


A AAO 

—0.002 


A 1 OO 


A AO 

0.48 




A A1 IT 1 A A1 A 

— 0.017± 0.019 


A AAO 

—0.002 


0.144 


0.49 




f\ r\r\ A 1 AA"tA 

— 0.004ib 0.019 


—0.002 


A 1 /^O 

O.lbz 


A r A 

0.50 




AAir 1 AA1A 

0.015± 0.019 


A AAO 

—0.002 


A 1 OO 

0.182 


ATI 

0.51 




A AOO 1 A A1 A 

0.028± 0.019 


A AAO 

—0.002 


A OAr" 

0.205 


0.52 




A A1 O 1 A AOA 

0.012± 0.020 


A AAf 

—0.005 


A OO 1 

0.231 


A r" o 

0.53 




A A1 1 A A01 

O.Oloi 0.021 


A AAO 

—0.003 


0.259 


0.53 




A AAA 1 A AOO 

0.009± 0.022 


0.005 


u.zyz 


U.o4 




U.UZoit U.UZ4 


U.UUD 


0.329 


0.55 




0.031± 0.026 


-0.004 


0.370 


0.55 




-0.032± 0.027 


-0.009 


0.417 


0.56 




0.029± 0.030 


-0.014 


0.469 


0.56 




-0.039± 0.032 


-0.019 


0.527 


0.57 




-0.005± 0.045 


-0.029 


0.594 


0.57 




-O.OOSzb 0.074 


-0.036 


0.669 


0.57 




-0.021± 0.110 


-0.039 


0.753 


0.58 




O.OOli 0.117 


-0.042 


0.847 


0.58 




O.OlOib 0.142 


-0.046 


e = 7.0° 


0.063 


0.60 




0.116± 0.133 


-0.005 


0.071 


0.64 




-0.036± 0.054 


-0.005 


0.080 


0.69 




0.013± 0.037 


-0.004 


0.090 


0.74 




-0.013± 0.028 


-0.003 


0.101 


0.78 




0.016± 0.024 


-0.003 


0.113 


0.82 




0.024± 0.022 


-0.003 


0.128 


0.86 




0.052± 0.021 


-0.002 


0.144 


0.90 




0.012± 0.020 


-0.002 


0.162 


0.93 




0.036± 0.020 


-0.001 
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0.182 


0.97 


0.027± 


0.019 


-0.001 


0.205 


1.00 


0.005± 


0.019 


0.000 


0.231 


1.03 


0.044± 


0.019 


0.000 


0.259 


1.06 


0.019± 


0.020 


0.001 


0.292 


1.09 


0.021± 


0.020 


0.002 


0.329 


1.12 


0.034± 


0.021 


0.002 


0.370 


1.14 


0.036± 


0.022 


-0.003 


0.417 


1.16 


0.068± 


0.023 


-0.001 


0.469 


1.19 


0.020± 


0.025 


0.007 


0.527 


1.21 


0.024± 


0.029 


0.001 


0.594 


1.22 


0.022± 


0.033 


-0.011 


0.669 


1.24 


0.028± 


0.039 


-0.029 


0.753 


1.25 


O.OOOib 


0.068 


-0.054 


0.847 


1.26 


-0.018± 


0.118 


-0.066 
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TABLE IX. Results for gi/Fi in the DIS region {W"^ > 4 GeV^). There is an additional 
normalization uncertainty due to beam and target polarization shown in Table XV. 



X < Q2 > E(GeV) 5f/Ff ibstat ibsyst ± stat ± syst ± stat ± syst 



0.024 





51 


16 


2 





092 


± 


027 


zb 





014 


-0 


048 


zb 





068 


zb 


009 


-0 


205 


zb 


153 


zb 





023 


0.027 





55 


16 


2 





086 


± 


026 


zb 





013 


-0 


018 


zb 





065 


zb 


009 


-0 


131 


zb 


146 


zb 





022 


0.027 


1 


17 


29 


1 


-0 


032 


± 


068 


zb 





010 


-0 


058 


zb 





092 


lb 


009 


-0 


099 


zb 


225 


zb 





020 


0.031 





59 


16 


2 





084 


± 


026 


zb 





012 





019 


zb 





065 


lb 


008 


-0 


048 


zb 


144 


zb 





020 


0.031 


1 


27 


29 


1 





064 


± 


034 


zb 





009 





041 


zb 





040 


zb 


008 





021 


zb 


095 


zb 





019 


0.035 





31 


9 


7 





096 


± 


037 


zb 





018 


-0 


027 


zb 





055 


zb 


008 


-0 


161 


zb 


125 


zb 





024 


0.035 





64 


16 


2 





059 


± 


024 


zb 





Oil 


-0 


015 


zb 





062 


zb 


007 


-0 


097 


zb 


139 


zb 





018 


0.035 


1 


40 


29 


1 





123 


± 


025 


zb 





008 





008 


zb 





029 


zb 


007 


-0 


115 


zb 


069 


zb 





017 


0.039 





33 


9 


7 





024 


± 


037 


zb 





018 





035 


zb 





053 


zb 


008 





050 


zb 


122 


zb 





024 


0.039 





68 


16 


2 





072 


± 


024 


zb 





Oil 





099 


zb 





060 


zb 


007 





140 


zb 


136 


zb 





017 


0.039 


1 


52 


29 


1 





083 


± 


023 


zb 





008 





062 


zb 





026 


zb 


007 





049 


zb 


063 


zb 





016 


0.044 





35 


9 


7 





106 


± 


037 


zb 





019 





003 


zb 





054 


zb 


008 


-0 


109 


zb 


126 


zb 





024 


0.044 





73 


16 


2 





096 


± 


023 


zb 





Oil 


-0 


038 


zb 





052 


zb 


006 


-0 


191 


zb 


120 


zb 





017 


0.044 





98 


16 


2 





097 


± 


520 


zb 





010 


-0 


943 


zb 


1 


296 


zb 


008 


-2 


222 


zb 2 


964 


zb 





019 


0.044 


1 


65 


29 


1 





110 


± 


021 


zb 





008 


-0 


013 


zb 





025 


zb 


006 


-0 


150 


zb 


060 


zb 





014 


0.049 





36 


9 


7 





108 


± 


039 


zb 





020 





055 


zb 





057 


zb 


008 





003 


zb 


134 


zb 





025 


0.049 





78 


16 


2 





147 


± 


024 


zb 





Oil 





092 


zb 





052 


zb 


006 





045 


zb 


121 


zb 





016 


0.049 


1 


06 


16 


2 





078 


± 


109 


zb 





010 


-0 


104 


zb 





280 


zb 


008 


-0 


321 


zb 


646 


zb 





018 


0.049 


1 


78 


29 


1 





125 


± 


020 


zb 





008 


-0 


013 


zb 





024 


zb 


005 


-0 


168 


zb 


058 


zb 





013 


0.056 





38 


9 


7 





108 


± 


041 
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TABLE X. Results for Ai in the DIS region {W"^ > 4 GeV^) 
normalization uncertainty due to beam and target polarization shown in Table PCV 
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29.1 
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9.7 
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1.22 


16.2 
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2.14 


16.2 
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3.30 


29.1 
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5.13 


29.1 


0.182 
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9.7 


0.182 


1.25 


16.2 
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2.25 


16.2 
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3.45 


29.1 
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1.00 


9.7 
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1.29 


16.2 
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2.36 


16.2 


0.205 


3.59 


29.1 


0.205 


5.86 


29.1 


0.230 


1.03 


9.7 



0.057 ±0.055 ±0.027 
0.140 ±0.043 ±0.014 
0.175 ±0.034 ±0.013 
0.194 ±0.037 ±0.010 
0.160 ± 0.023 ± 0.010 
0.152 ±0.040 ±0.010 
0.087 ± 0.057 ± 0.030 
0.132 ±0.040 ±0.014 
0.217 ±0.037 ±0.012 
0.183 ±0.037 ±0.009 
0.237 ±0.026 ±0.011 
0.199 ±0.033 ±0.010 
0.167 ±0.059 ±0.033 
0.185 ±0.040 ±0.015 
0.164 ±0.041 ±0.013 
0.206 ± 0.037 ± 0.009 
0.243 ± 0.028 ± 0.012 
0.224 ±0.031 ±0.010 
-0.010 ± 0.064 ± 0.036 
0.129 ± 0.041 ± 0.015 
0.217 ±0.045 ±0.013 
0.233 ± 0.037 ± 0.010 
0.228 ±0.031 ±0.013 
0.243 ± 0.031 ± 0.010 
0.134 ±0.042 ±0.016 
0.218 ± 0.049 ± 0.013 
0.258 ±0.038 ±0.011 
0.227 ± 0.033 ± 0.014 
0.242 ± 0.030 ± 0.010 
0.205 ± 0.043 ± 0.017 
0.290 ± 0.054 ± 0.014 
0.237 ± 0.040 ± 0.013 
0.297 ±0.037 ±0.015 
0.269 ±0.030 ±0.011 
0.298 ± 0.045 ± 0.018 
0.237 ± 0.060 ± 0.014 
0.258 ± 0.043 ± 0.014 
0.326 ±0.041 ±0.015 
0.309 ±0.031 ±0.011 
0.251 ± 0.048 ± 0.017 
0.305 ± 0.066 ± 0.016 
0.304 ± 0.046 ± 0.015 
0.278 ± 0.045 ± 0.015 
0.370 ± 0.033 ± 0.012 
0.335 ± 0.050 ± 0.018 



0.028 ± 0.078 ± 0.018 

-0.030 ± 0.062 ± 0.007 
0.109 ±0.071 ±0.006 
0.208 ± 0.085 ± 0.005 
0.070 ± 0.027 ± 0.004 
0.061 ± 0.045 ± 0.004 
0.154 ±0.081 ±0.023 
0.037 ± 0.058 ± 0.008 
0.057 ± 0.077 ± 0.006 
0.097 ± 0.083 ± 0.004 
0.064 ± 0.030 ± 0.004 
0.055 ± 0.038 ± 0.004 

-0.012 ±0.086 ±0.028 
0.061 ±0.057 ±0.009 
0.052 ± 0.083 ± 0.007 
0.051 ± 0.078 ± 0.004 
0.079 ± 0.033 ± 0.005 
0.035 ± 0.036 ± 0.004 

-0.087 ±0.092 ±0.031 
0.143 ± 0.059 ± 0.011 
0.229 ± 0.090 ± 0.008 
0.037 ± 0.077 ± 0.005 
0.187 ±0.037 ±0.006 
0.080 ± 0.035 ± 0.005 
0.033 ± 0.060 ± 0.012 
0.199 ± 0.099 ± 0.009 
0.085 ± 0.078 ± 0.006 
0.141 ±0.041 ±0.007 
0.171 ±0.035 ±0.005 
0.115 ±0.063 ±0.014 

-0.108 ±0.108 ±0.011 
0.137 ±0.082 ±0.007 
0.109 ±0.045 ±0.008 
0.084 ± 0.035 ± 0.005 
0.094 ± 0.067 ± 0.016 
0.200 ± 0.121 ± 0.012 
0.237 ± 0.087 ± 0.008 
0.175 ±0.051 ±0.009 
0.214 ±0.037 ±0.006 
0.017 ±0.071 ±0.016 
0.213 ±0.135 ±0.014 
0.117 ±0.094 ±0.010 
0.159 ±0.057 ±0.010 
0.132 ±0.039 ±0.008 
0.178 ± 0.075 ± 0.017 



0.000 ±0.193 ±0.046 
-0.235 ±0.156 ±0.020 
0.054 ±0.173 ±0.018 
0.268 ±0.207 ±0.013 
-0.022 ± 0.071 ± 0.013 
-0.033 ±0.117 ±0.014 
0.262 ±0.203 ±0.055 
-0.069 ±0.147 ±0.021 
-0.121 ±0.189 ±0.018 
0.016 ± 0.205 ± 0.012 
-0.129 ±0.079 ±0.014 
-0.106 ±0.099 ±0.013 
-0.226 ±0.215 ±0.064 
-0.075 ±0.147 ±0.023 
-0.071 ± 0.208 ± 0.019 
-0.126 ±0.196 ±0.012 
-0.103 ±0.088 ±0.016 
-0.187 ±0.094 ±0.013 
-0.196 ±0.233 ±0.072 
0.193 ±0.152 ±0.026 
0.294 ±0.227 ±0.021 
-0.197 ±0.194 ±0.014 
0.176 ±0.098 ±0.018 
-0.106 ±0.093 ±0.014 
-0.085 ± 0.158 ± 0.029 
0.218 ±0.251 ±0.023 
-0.114 ±0.198 ±0.016 
0.064 ±0.109 ±0.020 
0.118 ±0.094 ±0.014 
0.027 ±0.166 ±0.033 
-0.638 ± 0.279 ± 0.026 
0.039 ±0.211 ±0.019 
-0.112 ±0.124 ±0.021 
-0.138 ±0.097 ±0.015 
-0.148 ± 0.178 ± 0.036 
0.198 ±0.315 ±0.029 
0.261 ±0.227 ±0.022 
0.013 ±0.140 ±0.023 
0.132 ±0.102 ±0.017 
-0.285 ±0.190 ±0.037 
0.137 ±0.355 ±0.032 
-0.109 ± 0.250 ±0.024 
0.032 ±0.160 ±0.025 
-0.167 ±0.110 ±0.019 
0.005 ± 0.204 ± 0.039 
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TABLE XI. Coefficients for fits to all available data with > Q^m the form 

arc"(l + bx + cx^)[l + Cf{Q^)], along with the for the indicated number of degrees of freedom, 
calculated with the statistical errors only. Fits I to IV are to gi/Fi, while fit V is to Ai. 



fit to 






a 


a 


b 


c 


C 


x' 


df 




0.3 


none 


0.62 


0.641 


2.231 


-2.666 


0.000 


145 


117 


n. g'JFf 


1.0 


none 


0.64 


0.749 


1.466 


-1.982 


0.000 


112 


106 


III. g{/Ff 


0.3 


i/Q' 


0.62 


0.762 


1.434 


-1.917 


-0.160 


116 


116 


IV. g^Ff 


0.3 


-HQ') 


0.66 


0.728 


0.850 


-1.384 


-0.100 


120 


116 


V. A{ 


0.3 


i/Q' 


0.66 


0.898 


0.595 


-0.371 


-0.180 


118 


116 


I. gf/Ff 


0.3 


none 


1.52 


2.439 
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0.867 


0.000 


122 


111 


11. gf/Ff 


1.0 


none 
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0.829 


0.000 


115 


100 


III. gf/Ff 

IV. gf/Ff 


0.3 
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-HQ') 


1.44 


2.342 
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0.902 


-0.260 


119 
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0.3 


1.48 


2.030 


-1.812 


0.979 


-0.100 


120 


110 


V. Af 


0.3 


i/Q' 


1.46 


2.493 


-1.915 


1.376 


-0.260 


119 


110 
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TABLE XII. Results for averaged gi/Fi for > 1 (GeV/c)^. 



X 


<Q2 > 


q? /F? ± Stat ± syst 


of/Ff± Stat ± syst 




n / -t Stat ± syst 


0.031 


1.27 
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0.276 ± 0.019 ± 0.012 


0.090 ± 0.025 ± 0.006 


-0. 


,139 ± 0.068 ± 0.017 


0.182 


4.34 


0.296 ±0.020 ±0.013 


0.200 ± 0.027 ± 0.007 


0. 


,102 ±0.074 ±0.019 


0.205 


4.15 


0.319 ±0.020 ±0.014 


0.109 ±0.027 ±0.009 


-0. 


,165 ±0.075 ±0.023 


0.230 


4.37 


0.322 ± 0.021 ± 0.015 


0.130 ±0.029 ±0.011 


-0. 


,103 ±0.083 ±0.026 


0.259 


5.26 


0.434 ± 0.026 ± 0.015 


0.103 ± 0.036 ± 0.011 


-0. 


,369 ±0.104 ±0.026 


0.292 


5.53 


0.405 ± 0.029 ± 0.016 


0.227 ±0.041 ±0.012 


0. 


,029 ± 0.121 ± 0.029 


0.329 


6.01 


0.427 ± 0.035 ± 0.020 


0.284 ± 0.048 ± 0.014 


0. 


,130 ±0.147 ±0.034 


370 


6 99 






n 


1 QO _|_ n 1 70 _|_ A noo 


0.416 


6.56 


0.592 ±0.047 ±0.021 


0.232 ±0.067 ±0.017 


-0. 


,371 ±0.221 ±0.039 


0.468 


6.79 


0.608 ± 0.056 ± 0.024 


0.319 ±0.083 ±0.018 


-0. 


,126 ±0.282 ±0.042 


0.526 


7.72 


0.671 ± 0.083 ± 0.024 


0.455 ±0.113 ±0.019 


0. 


,141 ± 0.411 ± 0.045 


0.592 


7.97 


0.600 ±0.107 ±0.025 


0.322 ±0.148 ±0.019 


-0. 


,191 ±0.557 ±0.045 


0.666 


9.26 


0.478 ±0.165 ±0.026 


0.251 ±0.226 ±0.019 


-0. 


,192 ±0.876 ±0.046 


0.749 


9.52 


0.744 ± 0.237 ± 0.031 


0.646 ± 0.331 ± 0.033 


0. 


,569 ± 1.233 ± 0.073 
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TABLE XIII. Results for averaged Ai for > i (GeV/c)^ 



X 


<Q^> 


± Stat ± syst 


A'} ± Stat ± syst 


AV- ± Stat ± syst 


0.031 


1.27 


0.063 ± 0.034 ± 0.009 


0.050 ± 0.040 ± 0.008 


0.044 ± 0.095 ± 0.019 


0.035 


1.40 


0.122 ± 0.025 ± 0.008 


0.002 ± 0.028 ± 0.007 


-0.126 ±0.069 ±0.017 


0.039 


1.52 


0.083 ± 0.023 ± 0.008 


0.059 ± 0.026 ± 0.007 


0.041 ± 0.063 ± 0.016 


0.044 


1.65 


0.110 ± 0.021 ± 0.008 


-0.019 ± 0.025 ± 0.006 


-0.162 ±0.060 ±0.014 


0.049 


1.78 


0.123 ±0.020 ±0.008 


-0.016 ± 0.024 ± 0.005 


-0.175 ± 0.058 ± 0.013 


0.056 


1.91 


0.130 ±0.019 ±0.008 


0.023 ± 0.023 ± 0.005 


-0.091 ±0.057 ±0.012 


0.063 


2.04 


0.123 ±0.018 ±0.008 


0.019 ± 0.022 ± 0.004 


-0.105 ±0.056 ±0.012 


0.071 


2.19 


0.152 ±0.018 ±0.009 


0.028 ± 0.022 ± 0.004 


-0.110 ±0.056 ±0.012 


0.079 


2.41 


0.180 ± 0.018 ± 0.009 


0.024 ± 0.022 ± 0.004 


-0.165 ±0.057 ±0.012 


0.090 


2.55 


0.168 ±0.016 ±0.011 


0.080 ± 0.022 ± 0.004 


0.004 ± 0.055 ± 0.014 


0.101 


2.85 


0.214 ±0.016 ±0.011 


0.063 ± 0.022 ± 0.004 


-0.110 ±0.057 ±0.014 


0.113 


3.13 


0.217 ±0.017 ±0.011 


0.058 ± 0.022 ± 0.005 


-0.135 ±0.059 ±0.015 


0.128 


3.41 


0.232 ±0.017 ±0.011 


0.128 ±0.023 ±0.005 


0.025 ± 0.061 ± 0.016 


0.144 


3.71 


0.238 ±0.018 ±0.012 


0.154 ±0.024 ±0.006 


0.081 ± 0.065 ± 0.017 


0.162 


4.03 


0.272 ± 0.019 ± 0.012 


0.086 ± 0.026 ± 0.007 


-0.142 ± 0.069 ± 0.018 


0.182 


4.34 


0.294 ± 0.020 ± 0.013 


0.204 ± 0.027 ± 0.008 


0.115 ±0.075 ±0.020 


0.205 


4.15 


0.314 ±0.020 ±0.014 


0.123 ±0.027 ±0.010 


-0.122 ±0.076 ±0.024 


0.230 


4.37 


0.325 ± 0.022 ± 0.015 


0.142 ±0.030 ±0.011 


-0.075 ± 0.085 ± 0.027 


0.259 


5.26 


0.437 ± 0.026 ± 0.016 


0.071 ± 0.038 ± 0.011 


-0.457 ±0.108 ±0.027 


0.292 


5.53 


0.432 ± 0.030 ± 0.017 


0.245 ± 0.043 ± 0.013 


0.036 ±0.127 ±0.031 


0.329 


6.01 


0.441 ± 0.036 ± 0.021 


0.308 ±0.052 ±0.015 


0.180 ±0.157 ±0.036 


0.370 


6.29 


0.471 ± 0.041 ± 0.023 


0.319 ±0.061 ±0.016 


0.162 ±0.192 ±0.040 


0.416 


6.56 


0.616 ± 0.049 ± 0.022 


0.258 ± 0.074 ± 0.018 


-0.314 ± 0.242 ± 0.042 


0.468 


6.79 


0.654 ± 0.059 ± 0.026 


0.383 ± 0.093 ± 0.019 


0.019 ± 0.314 ± 0.046 


0.526 


7.72 


0.713 ± 0.090 ± 0.027 


0.456 ± 0.135 ± 0.022 


0.064 ± 0.481 ± 0.051 


0.592 


7.97 


0.649 ±0.118 ±0.030 


0.159 ±0.182 ±0.023 


-0.851 ± 0.677 ± 0.054 


0.666 


9.26 


0.612 ±0.182 ±0.033 


0.423 ± 0.284 ± 0.024 


0.112 ± 1.078 ±0.059 


0.749 


9.52 


0.914 ± 0.273 ± 0.041 


0.769 ± 0.443 ± 0.039 


0.613 ± 1.605 ± 0.088 
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TABLE XIV. Results for gi obtained from average gi/Fi for > 1 (GeV/c)^). Not included 
in the systematic errors listed are normalization uncertaintird shown in Table pCV|. 



< a; > 




X 


range 


<Q2 


> 






± stat ± syst 






gf ± stat ± syst 






g"^ lb stat lb syst 


0.031 


0, 


,029- 


■0, 


,033 


1, 


,27 


0, 


,248 


±0, 


,132 


±0, 


,034 


0, 


,150 


±0, 


.147 


±0. 


,030 


0. 


,077 


±0. 


,345 


±0, 


.074 


0.035 


0. 


,033- 


■0, 


.037 


1, 


.40 


0, 


.436 


±0, 


,089 


±0, 


,032 


0, 


,027 


±0, 


.097 


±0. 


,025 


-0. 


,378 


±0. 


,227 


±0, 


.063 


0.039 


0. 


,037- 


■0, 


.042 


1, 


.52 


0, 


.269 


±0, 


,073 


±0, 


,025 


0, 


,191 


±0, 


.080 


±0. 


,021 


0. 


,145 


±0. 


,188 


±0, 


.054 


0.044 


0. 


,042- 


■0, 


.047 


1, 


.65 


0, 


.322 


±0, 


.062 


±0, 


,022 


-0, 


.036 


±0, 


.068 


±0. 


,017 


-0. 


,401 


±0. 


,160 


±0, 


.043 


0.049 


0. 


,047- 


■0, 


,053 


1, 


.78 


0, 


.327 


±0, 


,053 


±0, 


,019 


-0, 


.034 


±0, 


.059 


±0. 


,013 


-0. 


,405 


±0. 


,139 


±0, 


.036 


0.056 


0. 


,053- 


•0, 


,059 


1, 


.91 


0, 


,309 


±0, 


,045 


±0, 


,017 


0, 


,039 


±0, 


.052 


±0. 


,010 


-0. 


,224 


±0. 


,121 


±0, 


.029 


0.063 


0, 


,059- 


■0, 


.067 


2, 


,04 


0, 


,266 


±0, 


,039 


±0, 


,015 


0, 


,027 


±0, 


.045 


±0. 


,008 


-0. 


,230 


±0. 


,106 


±0, 


.023 


0.071 


0, 


,067- 


■0, 


.075 


2, 


,19 


0, 


,294 


±0, 


,034 


±0, 


,014 


0, 


,053 


±0, 


.041 


±0. 


,006 


-0. 


,182 


±0. 


,095 


±0, 


.021 


0.079 


0, 


,075- 


■0, 


.084 


2, 


,41 


0, 


,310 


±0, 


,031 


±0, 


,014 


0, 


,039 


±0, 


.037 


±0. 


,005 


-0. 


,245 


±0. 


,086 


±0, 


.019 


0.090 


0. 


,084- 


■0, 


,095 


2, 


,55 


0, 


.260 


±0, 


,024 


±0, 


,012 


0, 


,123 


±0, 


.031 


±0. 


,005 


0. 


,015 


±0. 


,073 


±0, 


.019 


0.101 


0. 


,095- 


■0, 


,107 


2, 


,85 


0, 


,299 


±0, 


,022 


±0, 


,012 


0, 


,082 


±0, 


.029 


±0. 


,004 


-0. 


,130 


±0. 


,067 


±0, 


.017 


0.113 


0, 


,107- 


■0, 


,120 


3, 


,13 


0, 


,272 


±0, 


,021 


±0, 


,011 


0, 


,067 


±0, 


.026 


±0, 


,003 


-0, 


,142 


±0. 


,061 


±0, 


.015 


0.128 


0. 


,120- 


■0, 


,136 


3, 


,41 


0, 


,262 


±0, 


,019 


±0, 


.010 


0, 


,135 


±0, 


.024 


±0. 


,004 


0. 


,023 


±0. 


,057 


±0, 


.017 


0.144 


0. 


,136- 


■0, 


,153 


3, 


,71 


0, 


.239 


±0, 


,018 


±0, 


,010 


0, 


,138 


±0, 


.022 


±0. 


,004 


0. 


,056 


±0. 


,052 


±0, 


.017 


0.162 


0. 


,153- 


■0, 


,172 


4, 


.03 


0, 


.253 


±0, 


,017 


±0, 


,010 


0, 


,075 


±0, 


.021 


±0. 


,003 


-0. 


,100 


±0. 


,049 


±0, 


.014 


0.182 


0. 


,172- 


■0, 


,193 


4, 


.34 


0, 


.243 


±0, 


,016 


±0, 


.009 


0, 


,147 


±0, 


.020 


±0. 


,005 


0. 


,064 


±0. 


,046 


±0, 


.018 


0.205 


0. 


,193- 


■0, 


.218 


4, 


.15 


0, 


.231 


±0, 


,014 


±0, 


.009 


0, 


,070 


±0, 


.017 


±0. 


,005 


-0. 


,088 


±0. 


,040 


±0, 


.015 


0.230 


0. 


,218- 


■0, 


.245 


4, 


.37 


0, 


,206 


±0, 


,014 


±0, 


,008 


0, 


,072 


±0, 


.016 


±0, 


,005 


-0, 


,047 


±0, 


,038 


±0, 


.014 


0.259 


0. 


,245- 


■0, 


.276 


5, 


.26 


0, 


,242 


±0, 


,014 


±0, 


,008 


0, 


,049 


±0, 


.017 


±0, 


,004 


-0, 


,140 


±0, 


,040 


±0, 


.013 


0.292 


0, 


,276- 


■0, 


,310 


5, 


.53 


0, 


,192 


± 0, 


,014 


± 0, 


,007 


0, 


,089 


± 0, 


.016 


± 0, 


,004 


0, 


,009 


± 0, 


,037 


± 0, 


.013 


0.329 


0, 


,310- 


■0, 


,349 


6, 


.01 


0, 


,168 


±0, 


,014 


±0, 


,006 


0, 


,091 


±0, 


.015 


±0, 


,004 


0, 


,032 


±0, 


,036 


±0, 


.013 


0.370 


0, 


,349- 


■0, 


,393 


6, 


,29 


0, 


,144 


±0, 


,013 


±0, 


,005 


0, 


,076 


±0, 


.014 


±0, 


,004 


0, 


,026 


±0, 


,034 


±0, 


.011 


0.416 


0. 


,393- 


■0, 


.442 


6, 


,56 


0, 


,148 


±0, 


,012 


±0, 


.005 


0, 


,045 


±0, 


.013 


±0. 


,003 


-0. 


,051 


±0. 


,030 


±0, 


.009 


0.468 


0. 


,442- 


■0, 


.497 


6, 


.79 


0, 


,114 


±0, 


,011 


±0, 


.004 


0, 


.045 


±0, 


.012 


±0. 


,002 


-0. 


,012 


±0. 


,027 


±0, 


.007 


0.526 


0. 


,497- 


■0, 


,559 


7, 


.72 


0, 


.084 


±0, 


,010 


±0, 


.003 


0, 


,043 


±0, 


.011 


±0. 


,002 


0. 


,009 


±0. 


,025 


±0, 


.006 


0.592 


0. 


,559- 


■0, 


,629 


7, 


.97 


0, 


.048 


±0, 


,009 


±0, 


.002 


0, 


,019 


±0, 


.009 


±0. 


,001 


-0. 


,007 


±0. 


,021 


±0, 


.003 


0.666 


0, 


,629- 


■0, 


.708 


9, 


.26 


0, 


,021 


±0, 


,007 


±0, 


.001 


0, 


,008 


±0, 


.007 


±0. 


,001 


-0. 


,004 


±0. 


,017 


±0, 


.001 


0.749 


0, 


,708- 


■0, 


.791 


9, 


,52 


0, 


,015 


±0, 


,005 


±0, 


,001 


0, 


,009 


±0, 


.005 


±0, 


,001 


0, 


,005 


±0, 


,011 


±0, 


.001 



TABLE XV. Normalization systematic uncertainties 



Parameter Proton Deuteron 
Beam Polarization 2.4 2.4 
Target Polarization 2.5 4.0 
Dilution Factor 1.2 1.4 
Nitrogen Correction 0.4 

Total 3?7 49 
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TABLE XVI. Significant systematic errors for gf for E=29 GeV. The beam and target 



polarization errors are each completely correlated and the F2 error includes a correlated 
normalization error of ~ 1%. 



X 




f 


F2 


R 


Pb 


Pt 


A 


TOTAL 


e = 4.5° 



0.031 1.30 0.0010 0.0009 0.0016 

0.035 1.40 0.0010 0.0009 0.0014 

0.050 1.80 0.0011 0.0012 0.0011 

0.080 2.30 0.0013 0.0015 0.0007 

0.125 2.90 0.0014 0.0017 0.0005 

0.175 3.40 0.0014 0.0018 0.0003 

0.250 3.80 0.0013 0.0017 0.0001 

0.350 4.10 0.0009 0.0014 0.0001 

0.500 4.40 0.0005 0.0008 0.0000 



e = 7.0° 



0.080 


3.40 


0.0013 


0.0016 


0.0029 


0.0024 


0.0040 


0.0067 


0.0089 


0.125 


4.40 


0.0014 


0.0017 


0.0016 


0.0026 


0.0044 


0.0027 


0.0064 


0.175 


5.30 


0.0014 


0.0018 


0.0009 


0.0026 


0.0044 


0.0028 


0.0063 


0.250 


6.40 


0.0013 


0.0017 


0.0005 


0.0023 


0.0039 


0.0035 


0.0061 


0.350 


7.50 


0.0009 


0.0013 


0.0003 


0.0017 


0.0028 


0.0032 


0.0048 


0.500 


8.60 


0.0004 


0.0006 


0.0001 


0.0007 


0.0012 


0.0018 


0.0024 


0.700 


9.30 


0.0001 


0.0001 


0.0000 


0.0002 


0.0003 


0.0003 


0.0005 



0.0018 
0.0018 
0.0020 
0.0023 
0.0025 
0.0025 
0.0023 
0.0017 
0.0009 



0.0030 
0.0030 
0.0033 
0.0038 
0.0041 
0.0042 
0.0038 
0.0029 
0.0015 



0.0298 

0.0249 
0.0130 
0.0044 
0.0024 
0.0031 
0.0039 
0.0037 
0.0018 



0.0301 

0.0252 
0.0137 
0.0066 
0.0058 
0.0062 
0.0063 
0.0053 
0.0027 
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TABLE XVII. Significant systematic errors for gf/Ff for E=29 GeV. The beam and target 
polarization errors are each completely correlated. 



X 


Q' 


/ 


R 


Pb 


Pt 


Arc 


TOTAL 


= 4.5° 


0.031 


1.30 


0.0003 


0.0005 


0.0005 


0.0008 


0.0080 


0.0081 


0.035 


1.40 


0.0003 


0.0006 


0.0005 


0.0009 


0.0074 


0.0075 


0.050 


1.80 


0.0004 


0.0009 


0.0008 


0.0013 


0.0052 


0.0055 


0.080 


2.30 


0.0008 


0.0020 


0.0014 


0.0023 


0.0027 


0.0044 


0.125 


2.90 


0.0013 


0.0040 


0.0023 


0.0039 


0.0023 


0.0066 


0.175 


3.40 


0.0018 


0.0053 


0.0033 


0.0055 


0.0041 


0.0095 


0.250 


3.80 


0.0025 


0.0063 


0.0046 


0.0077 


0.0080 


0.0138 


0.350 


4.10 


0.0032 


0.0077 


0.0059 


0.0099 


0.0124 


0.0189 


0.500 


4.40 


0.0037 


0.0084 


0.0068 


0.0114 


0.0139 


0.0213 


e = 7.0° 


0.080 


3.40 


0.0008 


0.0017 


0.0014 


0.0023 


0.0039 


0.0051 


0.125 


4.40 


0.0013 


0.0030 


0.0024 


0.0039 


0.0024 


0.0061 


0.175 


5.30 


0.0019 


0.0041 


0.0034 


0.0056 


0.0035 


0.0087 


0.250 


6.40 


0.0026 


0.0060 


0.0047 


0.0078 


0.0071 


0.0132 


0.350 


7.50 


0.0032 


0.0082 


0.0060 


0.0099 


0.0115 


0.0185 


0.500 


8.60 


0.0037 


0.0099 


0.0069 


0.0114 


0.0166 


0.0238 


0.700 


9.30 


0.0037 


0.0101 


0.0069 


0.0115 


0.0132 


0.0217 
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TABLE XVIII. Significant systematic errors for for E=29 GeV. The beam and target 



polarization errors are each completely correlated and the F2 error includes a correlated 
normalization error of ~ 1%. 



X 


Q' 


/ 


F2 


R 


Pb 


Pt 


A 


TOTAL 


e = 4.5° 


0.031 


1.30 


0.0052 


0.0062 


0.0104 


0.0071 


0.0074 


0.0320 


0.0361 


0.035 


1.40 


0.0053 


0.0061 


0.0085 


0.0070 


0.0073 


0.0269 


0.0310 


0.050 


1.80 


0.0051 


0.0059 


0.0049 


0.0068 


0.0071 


0.0156 


0.0206 


0.080 


2.30 


0.0047 


0.0056 


0.0023 


0.0064 


0.0067 


0.0094 


0.0153 


0.125 


2.90 


0.0043 


0.0055 


0.0014 


0.0060 


0.0063 


0.0060 


0.0128 


0.175 


3.40 


0.0041 


0.0054 


0.0008 


0.0057 


0.0060 


0.0047 


0.0117 


0.250 


3.80 


0.0036 


0.0050 


0.0003 


0.0052 


0.0054 


0.0038 


0.0104 


0.350 


4.10 


0.0028 


0.0039 


0.0001 


0.0041 


0.0043 


0.0033 


0.0083 


0.500 


4.40 


0.0016 


0.0022 


0.0000 


0.0024 


0.0025 


0.0017 


0.0047 


e = 7.0° 


0.080 


3.40 


0.0046 


0.0059 


0.0095 


0.0066 


0.0069 


0.0130 


0.0202 


0.125 


4.40 


0.0046 


0.0057 


0.0042 


0.0064 


0.0067 


0.0070 


0.0143 


0.175 


5.30 


0.0043 


0.0056 


0.0022 


0.0060 


0.0063 


0.0043 


0.0122 


0.250 


6.40 


0.0037 


0.0050 


0.0012 


0.0053 


0.0055 


0.0027 


0.0102 


0.350 


7.50 


0.0027 


0.0037 


0.0006 


0.0040 


0.0041 


0.0020 


0.0076 


0.500 


8.60 


0.0014 


0.0018 


0.0002 


0.0020 


0.0021 


0.0012 


0.0039 


0.700 


9.30 


0.0003 


0.0005 


0.0000 


0.0005 


0.0005 


0.0004 


0.0010 
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TABLE XIX. Significant systematic errors for (^/Ff for E=29 GeV. The beam and target 
polarization errors are each completely correlated. 



X 




/ 


R 


Pb 


Pt 


A 


TOTAL 


6 = 4.5° 


0.031 


1.30 


0.0013 


0.0030 


0.0018 


0.0019 


0.0082 


0.0092 


0.035 


1.40 


0.0015 


0.0034 


0.0020 


0.0021 


0.0076 


0.0089 


0.050 


1.80 


0.0019 


0.0040 


0.0026 


0.0027 


0.0059 


0.0083 


0.080 


2.30 


0.0027 


0.0062 


0.0037 


0.0039 


0.0055 


0.0102 


0.125 


2.90 


0.0037 


0.0095 


0.0052 


0.0055 


0.0052 


0.0138 


0.175 


3.40 


0.0048 


0.0111 


0.0068 


0.0071 


0.0055 


0.0165 


0.250 


3.80 


0.0062 


0.0122 


0.0089 


0.0093 


0.0066 


0.0199 


0.350 


4.10 


0.0079 


0.0146 


0.0114 


0.0119 


0.0091 


0.0251 


0.500 


4.40 


0.0096 


0.0175 


0.0140 


0.0146 


0.0100 


0.0301 


e = 7.0° 


0.080 


3.40 


0.0025 


0.0051 


0.0036 


0.0038 


0.0072 


0.0106 


0.125 


4.40 


0.0038 


0.0073 


0.0053 


0.0055 


0.0057 


0.0126 


0.175 


5.30 


0.0048 


0.0086 


0.0069 


0.0071 


0.0049 


0.0148 


0.250 


6.40 


0.0063 


0.0115 


0.0090 


0.0094 


0.0045 


0.0190 


0.350 


7.50 


0.0079 


0.0156 


0.0115 


0.0120 


0.0057 


0.0248 


0.500 


8.60 


0.0096 


0.0206 


0.0142 


0.0148 


0.0081 


0.0317 


0.700 


9.30 


0.0100 


0.0229 


0.0149 


0.0156 


0.0129 


0.0355 
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TABLE XX. Deuteron results for gi at fixed of 2, 3, and 5 (GeV/c)^ evaluated assuming 
gi/Fi is independent of Q^. In addition to the systematic errors shown, there are normalization 
uncertainties shown in Table XV. 



X 


gfiQ' = 


2) lb stat lb syst 


5f (Q2 = 


3) lb stat lb syst 


gf{Q' = 


5) zb stat zb syst 


0.031 


0.166 


lb 0.162 lb 0.033 


0.180 


lb 0.176 lb 0.036 


0.198 


zb 0.194 zb 0.040 


0.035 


0.029 


lb 0.104 lb 0.027 


0.031 


lb 0.112 lb 0.029 


0.034 


zb 0.124 zb 0.032 


0.039 


0.202 


lb 0.084 lb 0.022 


0.218 


lb 0.091 lb 0.024 


0.239 


zb 0.100 zb 0.027 


0.044 


-0.038 


lb 0.071 lb 0.017 


-0.040 


lb 0.076 lb 0.019 


-0.044 


zb 0.083 zb 0.020 


0.049 


-0.035 


lb 0.060 lb 0.014 


-0.037 


lb 0.064 lb 0.014 


-0.040 


zb 0.070 zb 0.016 


0.056 


0.040 


lb 0.052 ±0.011 


0.042 


lb 0.056 lb 0.011 


0.046 


zb 0.061 zb 0.012 


0.063 


0.027 


lb 0.045 lb 0.008 


0.028 


lb 0.048 lb 0.008 


0.031 


zb 0.052 zb 0.009 


0.071 


0.053 


lb 0.040 lb 0.006 


0.056 


lb 0.042 lb 0.007 


0.060 


zb 0.046 zb 0.007 


0.079 


0.038 


lb 0.036 lb 0.005 


0.040 


zb 0.038 lb 0.005 


0.044 


zb 0.041 zb 0.006 


0.090 


0.120 


lb 0.030 lb 0.005 


0.125 


lb 0.032 zb 0.005 


0.135 


zb 0.034 zb 0.005 


0.101 


0.079 


lb 0.028 lb 0.004 


0.083 


zb 0.029 zb 0.004 


0.089 


zb 0.031 zb 0.004 


0.113 


0.064 


lb 0.025 lb 0.003 


0.067 


zb 0.026 zb 0.004 


0.072 


zb 0.028 zb 0.004 


0.128 


0.128 


lb 0.023 lb 0.004 


0.133 


zb 0.024 zb 0.004 


0.141 


zb 0.025 zb 0.004 


0.144 


0.130 


lb 0.021 lb 0.004 


0.135 


zb 0.022 zb 0.004 


0.142 


zb 0.023 zb 0.004 


0.162 


0.071 


lb 0.020 lb 0.003 


0.073 


zb 0.020 zb 0.004 


0.076 


zb 0.021 zb 0.004 


0.182 


0.141 


lb 0.019 lb 0.004 


0.143 


zb 0.019 zb 0.004 


0.148 


zb 0.020 zb 0.005 


0.205 


0.068 


lb 0.017 lb 0.005 


0.069 


zb 0.017 zb 0.005 


0.070 


zb 0.017 zb 0.005 


0.230 


0.071 


lb 0.016 lb 0.005 


0.071 


zb 0.016 zb 0.005 


0.072 


zb 0.016 zb 0.005 


0.259 


0.049 


lb 0.017 lb 0.004 


0.049 


zb 0.017 zb 0.004 


0.049 


zb 0.017 zb 0.004 


0.292 


0.094 


lb 0.017 zb 0.004 


0.092 


zb 0.016 zb 0.004 


0.090 


zb 0.016 zb 0.004 


0.329 


0.100 


lb 0.017 lb 0.004 


0.096 


zb 0.016 zb 0.004 


0.092 


zb 0.016 zb 0.004 


0.370 


0.089 


lb 0.017 lb 0.004 


0.083 


zb 0.016 zb 0.004 


0.078 


zb 0.015 zb 0.004 


0.416 


0.057 


lb 0.016 lb 0.004 


0.051 


zb 0.015 zb 0.003 


0.047 


zb 0.014 zb 0.003 


0.468 


0.064 


lb 0.017 lb 0.004 


0.055 


zb 0.014 zb 0.003 


0.048 


zb 0.012 zb 0.002 


0.526 


0.074 


lb 0.019 lb 0.004 


0.060 


zb 0.015 zb 0.003 


0.049 


zb 0.012 zb 0.002 


0.592 


0.043 


lb 0.020 lb 0.002 


0.031 


zb 0.014 zb 0.002 


0.023 


zb 0.011 zb 0.001 


0.666 


0.026 


lb 0.024 lb 0.002 


0.017 


zb 0.016 zb 0.001 


0.011 


zb 0.010 zb 0.001 


0.749 


0.048 


lb 0.025 lb 0.002 


0.028 


zb 0.014 zb 0.001 


0.016 


zb 0.008 zb 0.001 
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TABLE XXI. Proton results for gi at fixed of 2, 3, and 5 (GeV/c)^ evaluated assumme; 
gi/Fi is independent of Q^. In addition to the systematic errors shown, there are normalization 
uncertainties shown in Table XV. 



X 


g'AQ' = 


2) lb stat lb syst 


gP{Q^ = 3) ± Stat ± syst 


g'AQ' = 


5) ± stat ± syst 


0.031 


0.272 


±0.145 ±0.037 


0.294 ±0.157 ±0.040 


0.322 


±0.172 ±0.044 


0.035 


0.468 


± 0.096 ± 0.032 


0.504 ±0.103 ±0.035 


0.552 


±0.113 ±0.041 


0.039 


0.283 


± 0.077 ± 0.026 


0.304 ± 0.083 ± 0.028 


0.333 


±0.091 ±0.031 


0.044 


0.334 


± 0.065 ± 0.022 


0.357 ± 0.069 ± 0.024 


0.390 


± 0.075 ± 0.027 


0.049 


0.334 


± 0.054 ± 0.019 


0.357 ± 0.058 ± 0.020 


0.389 


± 0.063 ± 0.023 


0.056 


0.312 


± 0.046 ± 0.017 


0.332 ± 0.049 ± 0.018 


0.362 


± 0.053 ± 0.020 


0.063 


0.267 


±0.039 ±0.014 


0.283 ±0.041 ±0.016 


0.308 


± 0.045 ± 0.017 


0.071 


0.292 


± 0.034 ± 0.013 


0.309 ± 0.036 ± 0.014 


0.336 


±0.039 ±0.016 


0.079 


0.305 


±0.030 ±0.013 


0.322 ± 0.032 ± 0.014 


0.349 


±0.034 ±0.015 


0.090 


0.257 


±0.024 ±0.011 


0.270 ± 0.025 ± 0.012 


0.292 


±0.027 ±0.013 


0.101 


0.292 


±0.022 ±0.011 


0.306 ±0.023 ±0.012 


0.330 


±0.025 ±0.013 


0.113 


0.264 


±0.020 ±0.010 


0.276 ±0.021 ±0.011 


0.296 


±0.022 ±0.011 


0.128 


0.253 


±0.019 ±0.010 


0.263 ±0.019 ±0.010 


0.281 


±0.021 ±0.010 


0.144 


0.229 


±0.017 ±0.009 


0.238 ±0.018 ±0.009 


0.252 


±0.019 ±0.009 


0.162 


0.242 


±0.016 ±0.009 


0.250 ±0.017 ±0.009 


0.263 


±0.018 ±0.009 


0.182 


0.232 


±0.016 ±0.008 


0.238 ± 0.016 ± 0.008 


0.248 


±0.017 ±0.009 


0.205 


0.225 


± 0.014 ± 0.008 


0.229 ±0.014 ±0.008 


0.236 


±0.014 ±0.009 


0.230 


0.202 


±0.013 ±0.008 


0.204 ± 0.013 ± 0.008 


0.208 


± 0.014 ± 0.008 


0.259 


0.241 


± 0.014 ± 0.008 


0.241 ±0.014 ±0.008 


0.242 


±0.014 ±0.008 


0.292 


0.197 


± 0.014 ± 0.007 


0.195 ±0.014 ±0.007 


0.192 


±0.014 ±0.007 


0.329 


0.179 


±0.015 ±0.006 


0.174 ±0.014 ±0.006 


0.169 


±0.014 ±0.006 


0.370 


0.161 


±0.014 ±0.005 


0.153 ±0.014 ±0.005 


0.146 


±0.013 ±0.005 


0.416 


0.177 


± 0.014 ± 0.005 


0.164 ±0.013 ±0.005 


0.151 


±0.012 ±0.005 


0.468 


0.151 


± 0.014 ± 0.005 


0.134 ±0.012 ±0.004 


0.118 


±0.011 ±0.004 


0.526 


0.139 


± 0.017 ± 0.004 


0.115 ±0.014 ±0.003 


0.095 


±0.012 ±0.003 


0.592 


0.105 


±0.019 ±0.003 


0.078 ± 0.014 ± 0.002 


0.059 


±0.011 ±0.002 


0.666 


0.070 


± 0.024 ± 0.002 


0.046 ±0.016 ±0.001 


0.030 


±0.010 ±0.001 


0.749 


0.080 


± 0.025 ± 0.002 


0.046 ±0.014 ±0.001 


0.026 


± 0.008 ± 0.001 
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TABLE XXII. Neutron results for gi at fixed of 2, 3, and 5 (GeV/c)^ evaluated from and 
gf assuming gi/Fi is independent of Q^. In addition there is a normalization uncertainty common 
to all data of 2.4% due to beam polarization. 



X 


5^(g2^2)±stat±syst 


5^(q2 = 3) -t Stat ± syst 


g\'{Q^ = 5) ± Stat ± syst 


0.031 


0.085 ± 0.378 ± 0.081 


0.091 ± 0.408 ± 0.088 


0.100 ±0.448 ±0.097 


0.035 


-0.406 ± 0.244 ± 0.067 


-0.437 ± 0.262 ± 0.072 


-0.479 ± 0.287 ± 0.080 


0.039 


0.153 ±0.197 ±0.056 


0.164 ±0.212 ±0.060 


0.180 ±0.232 ±0.066 


0.044 


-0.415 ±0.166 ±0.044 


-0.444 ±0.177 ±0.048 


-0.485 ±0.194 ±0.053 


0.049 


-0.413 ±0.141 ±0.036 


-0.440 ±0.151 ±0.039 


-0.480 ±0.164 ±0.043 


0.056 


-0.225 ±0.122 ±0.029 


-0.240 ±0.130 ±0.030 


-0.261 ±0.142 ±0.033 


0.063 


-0.229 ± 0.106 ± 0.023 


-0.243 ± 0.112 ± 0.024 


-0.264 ±0.122 ±0.027 


0.071 


-0.180 ±0.094 ±0.019 


-0.190 ±0.099 ±0.021 


-0.206 ±0.107 ±0.022 


0.079 


-0.240 ± 0.084 ± 0.017 


-0.252 ±0.088 ±0.018 


-0.273 ± 0.095 ± 0.020 


0.090 


0.015 ±0.071 ±0.017 


0.016 ±0.074 ±0.018 


0.017 ±0.080 ±0.019 


0.101 


-0.125 ±0.064 ±0.015 


-0.131 ±0.067 ±0.016 


-0.141 ±0.072 ±0.017 


0.113 


-0.136 ± 0.058 ± 0.013 


-0.141 ±0.061 ±0.014 


-0.151 ±0.065 ±0.015 


0.128 


0.022 ±0.054 ±0.015 


0.023 ± 0.056 ± 0.016 


0.024 ± 0.059 ± 0.017 


0.144 


0.053 ± 0.050 ± 0.015 


0.055 ±0.051 ±0.016 


0.058 ± 0.054 ± 0.016 


0.162 


-0.095 ± 0.047 ± 0.012 


-0.098 ± 0.048 ± 0.012 


-0.102 ±0.050 ±0.013 


0.182 


0.061 ± 0.044 ± 0.016 


0.062 ± 0.045 ± 0.016 


0.064 ± 0.046 ± 0.017 


0.205 


-0.086 ± 0.039 ± 0.014 


-0.087 ±0.039 ±0.014 


-0.089 ± 0.040 ± 0.014 


0.230 


-0.047 ± 0.038 ± 0.014 


-0.047 ± 0.038 ± 0.014 


-0.047 ± 0.038 ± 0.014 


0.259 


-0.144 ±0.041 ±0.012 


-0.141 ±0.040 ±0.012 


-0.140 ±0.040 ±0.012 


0.292 


0.009 ± 0.040 ± 0.013 


0.009 ±0.039 ±0.013 


0.009 ± 0.038 ± 0.013 


0.329 


0.036 ± 0.041 ± 0.013 


0.034 ± 0.039 ± 0.013 


0.032 ±0.037 ±0.012 


0.370 


0.031 ± 0.040 ± 0.012 


0.029 ±0.037 ±0.011 


0.027 ±0.035 ±0.011 


0.416 


-0.067 ± 0.040 ± 0.010 


-0.060 ± 0.036 ± 0.009 


-0.054 ± 0.032 ± 0.009 


0.468 


-0.018 ± 0.040 ± 0.009 


-0.015 ± 0.034 ± 0.008 


-0.013 ± 0.029 ± 0.007 


0.526 


0.016 ±0.047 ±0.009 


0.013 ± 0.037 ± 0.008 


0.010 ±0.029 ±0.006 


0.592 


-0.018 ± 0.051 ± 0.007 


-0.012 ± 0.036 ± 0.005 


-0.009 ± 0.026 ± 0.004 


0.666 


-0.015 ±0.067 ±0.005 


-0.009 ± 0.041 ± 0.003 


-0.006 ± 0.026 ± 0.002 


0.749 


0.034 ± 0.073 ± 0.007 


0.018 ± 0.039 ± 0.004 


0.009 ± 0.021 ± 0.002 
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TABLE XXIII. gi{x)dx at different by different metliods. The first error is statistical and 



the second is systematic. There are additional normalization uncertainties shown in Table |X 





method 


g2 = 2 (GeV/c)2 


q2 = 3 (GeV/c)2 


q2 = 5 (GeV/c)2 


Deuteron 


91 /Fi 


0.050 ± 0.004 ± 0.003 


0.046 ± 0.003 ± 0.003 


0.043 ± 0.003 ± 0.002 


Deuteron 


Ai 


0.047 ± 0.005 ± 0.003 


0.044 ± 0.004 ± 0.003 


0.043 ± 0.003 ± 0.002 


Proton 


91 /Fi 


0.129 ±0.004 ±0.006 


0.121 ±0.003 ±0.006 


0.117 ±0.003 ±0.006 


Proton 


Ai 


0.120 ±0.004 ±0.006 


0.116 ±0.003 ±0.006 


0.116 ±0.003 ±0.006 


Neutron 


91 /Fi 


-0.022 ±0.011 ±0.006 


-0.023 ±0.008 ±0.006 


-0.025 ± 0.007 ± 0.006 


Neutron 


Ai 


-0.019 ±0.013 ±0.005 


-0.021 ±0.009 ±0.005 


-0.023 ± 0.007 ± 0.005 



TABLE XXIV. Systematic errors on the measured integral at = 3 (GeV/c)^. 



SOURCE Deuteron Proton Neutron p-n 

¥2 0.001 0.003 0.001 0.003 

R 0.000 0.001 0.000 0.001 

Arc 0.001 0.002 0.003 0.005 

f 0.001 0.002 0.001 0.003 

Pb 0.001 0.003 0.001 0.003 

Pt 0.002 0.003 0.005 0.007 

TOTAL 0.003 0.006 0.006 0.010 



TABLE XXV. Estimates of ^^'^^ gi{x,Q'^)dx at 
Columns 3-5 have the Regge form gi = I3x'^ fitted to gi I 
with a shown in the first row, in the range 0.03 < x < 
fit of the form ln{\/x). Column 7 has results of global 
integral of the SMC data with flat Regge extrapolation 



= 3 (GeV/c)^ using various hypotheses. 
X, Q^) at the shown in the second row, 
Xcut shown in column 2. Column 6 has a 



fit II of Table XI. The last column is the 
(a = 0) below x = 0.003. 





Xcut 


a = 


a = 


a = 0.5 


ln(l/x) 


global II 


SMC 






g2=3 


Q2=l 


g2=3 


g2=3 


g2=3 


g2=3 


Deuteron 


0.10 


0.002 


0.002 


0.001 


0.003 


0.001 


-0.005±0.003 


Proton 


0.10 


0.009 


0.010 


0.004 


0.015 


0.018 


0.014±0.003 


Neutron 


0.10 


-0.005 


-0.006 


-0.002 


-0.009 


-0.016 


-0.025±0.007 


p-n 


0.10 


0.014 


0.016 


0.006 


0.024 


0.034 


0.039±0.009 


Deuteron 


0.06 


0.001 


0.001 


0.000 


0.002 






Proton 


.06 


0.010 


0.011 


0.005 


0.015 






Neutron 


.06 


-0.008 


-0.009 


-0.004 


-0.012 






p-n 


.06 


0.018 


0.020 


0.010 


0.026 
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TABLE XXVI. Integral of gi in the measured region as well as extrapolations to high and low 
X as described in the text. Slight differences between the measured targets and derived targets, n 
and p-n, are due to correlations among systematic errors. The structure function gi was calculated 
at fixed assuming gi/Fi independent of Q^. 





Measured 


high X 


low X 


Total 


(GeV/c)2 


r.8 

J.03 91 


/.8 91 


Jo 91 


/o 51 



Deuteron 2 0.050 ± 0.004 ± 0.003 0.000 ± 0.001 0.001 ± 0.006 0.051 ± 0.004 ± 0.006 

Proton 2 0.129 ± 0.004 ± 0.006 0.001 ± 0.001 0.011 ± 0.007 0.140 ± 0.004 ± 0.010 

Neutron 2 -0.022 ± 0.011 ± 0.006 0.001 ± 0.001 -0.009 ± 0.016 -0.030 ± 0.011 ± 0.017 

p-n 2 0.149 ±0.012 ±0.011 0.001 ± 0.001 0.020 ± 0.019 0.169 ± 0.012 ± 0.022 

Deuteron 3 0.046 ± 0.003 ± 0.003 0.000 ± 0.001 0.001 ± 0.006 0.047 ± 0.003 ± 0.006 

Proton 3 0.121 ± 0.003 ± 0.006 0.001 ± 0.001 0.011 ± 0.007 0.133 ± 0.003 ± 0.009 

Neutron 3 -0.023 ± 0.008 ± 0.006 0.001 ± 0.001 -0.010 ± 0.015 -0.032 ± 0.008 ± 0.016 

p-n 3 0.143 ±0.009 ±0.010 0.001 ± 0.001 0.021 ± 0.018 0.164 ± 0.009 ± 0.021 

Deuteron 5 0.043 ± 0.003 ± 0.002 0.000 ± 0.001 0.001 ± 0.006 0.044 ± 0.003 ± 0.006 

Proton 5 0.117 ± 0.003 ± 0.006 0.001 ± 0.001 0.012 ± 0.008 0.129 ± 0.003 ± 0.010 

Neutron 5 -0.025 ± 0.007 ± 0.006 0.001 ± 0.001 -0.010 ± 0.015 -0.034 ± 0.007 ± 0.016 

p-n 5 0.141 ±0.008 ±0.010 0.001 ± 0.001 0.022 ± 0.017 0.164 ± 0.008 ± 0.020 



TABLE XXVII. Comparison of integrals from this experiment and E142 |Tg] and SMC 
Note that the SMC and E143 results used gi/Fi independent of and E142 used Ai independent 
of to evaluate gi at fixed from measurements at different Q^, and that the different 
experiments had different mean Q^. 

<— other experiments — > 
SMC 0.128 ±0.006 

0.140 ±0.011 
SMC 0.043 ± 0.007 

0.039 ± 0.008 
E142 -0.028 ± 0.008 
-0.031 ±0.011 







method 


X range 


this experiment 


Proton 


5 


91 /Fi 


0.03 < x < 0.7 


0.115 ±0.006 








< X < 1 


0.129 ±0.010 


Deuteron 


5 


91 /Fi 


0.03 < a; < 0.7 


0.041 ± 0.004 








< X < 1 


0.044 ± 0.007 


Neutron 


2 


Ai 


0.03 < a; < 0.6 


-0.021 ±0.009 








< X < 1 


-0.030 ± 0.020 
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TABLE XXVIII. Experimental value of Ti compared to the Ellis-Jaffe sum rule and Bjorken 
sum rule (p— n). For the theoretical input we take as{Mz) = 0.118 ±0.003 and 3F — D =0.58 with 
uncertainties of either 0.032 (small) or 0.120 (large). The Ellis-Jaffe sum is evaluated with both 
the invariant and Q^-dependent pQCD singlet corrections. The Bjorken sum rule depends only on 

the uoii-siu"i('t CDrrcxMioii. 





<Q2> 




ydheory 
^ 1 


-^theory 
^ 1 


error 


error 




(GeV/c)2 




invariant 


Q^-dependent 


(small) 


(large) 


Deuteron 


2 


0.051 ± 0.008 


0.070 


0.065 


±0.004 


±0.014 


Proton 


2 


0.140 ±0.010 


0.161 


0.156 


±0.005 


±0.016 


Neutron 


2 


-0.030 ± 0.020 


-0.010 


-0.015 


±0.005 


±0.016 


p-n 


2 


0.169 ± 0.025 


0.171 


0.171 


±0.006 


±0.006 


Deuteron 


3 


0.047 ±0.007 


0.071 


0.066 


±0.004 


±0.014 


Proton 


3 


0.133 ±0.010 


0.165 


0.160 


±0.005 


±0.016 


Neutron 


3 


-0.032 ±0.018 


-0.012 


-0.017 


±0.004 


±0.016 


p-n 


3 


0.164 ±0.023 


0.177 


0.177 


±0.004 


±0.004 


Deuteron 


5 


0.044 ± 0.007 


0.072 


0.068 


±0.004 


±0.015 


Proton 


5 


0.129 ± 0.010 


0.169 


0.164 


±0.005 


±0.016 


Neutron 


5 


-0.034 ±0.017 


-0.014 


-0.018 


±0.004 


±0.016 


p-n 


5 


0.164 ±0.021 


0.182 


0.182 


±0.003 


±0.003 



TABLE XXIX. The evaluated quark spins using both the "invariant" and Q^-dependent pQCD 
singlet coefficients, with an assumed error on F/D of 0.016. If the more conservative estimate of 
5{?>F — D) = 0.12 is used, the only change is to the error on As which is shown in the last column. 

method S = Oq Au Ad As Conservative 5As 

Deuterium invariant 0.35±0.07 0.84±0.02 -0.42±0.02 -0.08±0.03 ±0.05 

Proton invariant 0.29±0.09 0.83±0.03 -0.43±0.03 -0.10±0.03 ±0.06 

Deuterium ^3 (GeV/c)^ 0.37±0.08 0.85±0.03 -0.41±0.03 -0.07±0.03 ±0.06 

Proton Q2 ^3 (GeV/c)2 0.32±0.10 0.83±0.03 -0.43±0.03 -0.09±0.04 ±0.06 
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TABLE XXX. Results for A2, §2 and ^ for the proton measured in the 4.5° and 7.0° 
spectrometers at the indicated average values of x and and beam energy of 29.1 GeV. The 
highest X bin shown is in the resonance region defined by missing mass W'^ < 4 GeV^. 



X interval < x > < > 




p 
92 


5f 


(GeV/c)2 


drstat zbsyst 


ibstat ibsyst 


ibstat ibsyst 



0.029-0.047 0.038 1.49 0.016 ib 0.018 ib 0.006 0.489 ib 0.980 ib 0.332 0.223 ib 0.983 ib 0.332 

0.047- 0.075 0.060 2.01 0.025 i 0.014 i 0.005 0.397 i 0.374 i 0.138 0.223 i 0.375 i 0.138 

0.075 - 0.120 0.095 2.60 0.004 ± 0.015 ± 0.006 -0.236 ± 0.203 ± 0.074 -0.295 ± 0.204 ± 0.074 

0.120 - 0.193 0.152 3.21 0.021 ± 0.021 ± 0.008 -0.136 ± 0.125 ± 0.049 -0.127 ± 0.127 ± 0.049 

0.193-0.310 0.241 3.77 0.091 i 0.032 ib 0.011 -0.046 i 0.079 i 0.026 0.027 i 0.080 ib 0.026 

0.310-0.498 0.379 4.22 0.135 ± 0.060 ± 0.014 -0.050 ± 0.048 ± 0.009 0.051 ± 0.050 ± 0.009 

0.498 - 0.799 0.595 4.55 0.061 ib 0.154 ib 0.028 -0.037 ib 0.020 ib 0.004 -0.005 ib 0.022 ib 0.004 

0.075 - 0.120 0.101 3^76 0.025 ib 0.025 ib 0.007 0.060 ib 0.366 ib 0.116 0.000 ib 0.368 i 0.116 

0.120 -0.193 0.155 4.97 0.048 ib 0.019 ib 0.007 0.171 ib 0.141 ib 0.049 0.198 ib 0.142 ib 0.049 

0.193 -0.310 0.243 6.37 0.053 ib 0.022 ib 0.007 -0.068 ib 0.070 ib 0.021 0.023 ib 0.071 ib 0.021 

0.310 - 0.498 0.382 7.76 0.077 ± 0.035 ± 0.008 -0.039 ± 0.034 ± 0.007 0.043 ± 0.035 ± 0.007 

0.498-0.799 0.584 8.85 0.106 i 0.083 ± 0.016 -0.022 ib 0.011 i 0.002 0.001 ib 0.012 ± 0.002 



TABLE XXXI. Results for A2, 52 and 52 for the deuteron measured in the 4.5° and 7° 
spectrometers at the indicated average values of x and and beam energy of 29.1 GeV. The 
highest X bin shown is in the resonance region defined by missing mass W'^ < 4 GeV^. 



X interval 



<x> < > 
(GeV/c)2 



ibstat ibsyst 



ibstat ibsyst 
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ibstat ibsyst 



0.029 - 0.047 0.038 1.49 

0.047- 0.075 0.060 2.01 

0.075 -0.120 0.095 2.60 

0.120-0.193 0.152 3.21 

0.193-0.310 0.241 3.77 

0.310 - 0.498 0.378 4.22 

0.498 - 0.799 0.595 4.56 



0.070 lb 0.045 lb 0.010 3.426 ib 2.157 ib 0.575 3.275 ib 2.161 ib 0.575 

-0.025 lb 0.028 lb 0.006 -0.655 ib 0.707 ib 0.157 -0.799 ib 0.709 ib 0.157 

0.008 lb 0.032 lb 0.010 0.008 ib 0.390 ib 0.118 -0.048 ib 0.392 ib 0.118 

0.005 ± 0.045 ± 0.016 -0.118 ± 0.243 ± 0.080 -0.095 ± 0.245 ± 0.080 

0.078 lb 0.072 lb 0.020 0.127 i 0.154 ± 0.041 0.134 i 0.156 ib 0.041 
-0.079 lb 0.144 lb 0.017 -0.127 ib 0.094 ib 0.010 -0.095 ib 0.096 ib 0.010 

0.327 lb 0.390 lb 0.044 0.037 ib 0.039 ib 0.004 0.027 ib 0.041 ib 0.004 



0.075 -0.120 0.101 3.77 

0.120-0.193 0.154 4.97 

0.193-0.310 0.242 6.37 

0.310 -0.498 0.381 7.76 

0.498-0.799 0.584 8.86 



0.024 lb 0.046 lb 0.009 0.172 ib 0.621 ib 0.133 0.084 ib 0.624 ib 0.133 
-0.007 ± 0.036 ± 0.011 -0.109 ± 0.235 ± 0.076 -0.086 ± 0.237 ± 0.076 
-0.043 lb 0.043 ± 0.015 -0.133 ib 0.116 ib 0.041 -0.117 ± 0.118 ib 0.041 
0.000 lb 0.073 lb 0.014 -0.042 ib 0.056 ib 0.011 -0.006 ib 0.057 ib 0.011 
0.235 lb 0.183 lb 0.030 0.000 ib 0.018 i 0.003 0.011 i 0.019 i 0.003 
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TABLE XXXII. Results for A2, 52 and 52 for the neutron measured in the 4.5° and 7° 
spectrometers at the indicated average values of x and and beam energy of 29.1 GeV. The 
highest X bin shown is in the resonance region defined by missing mass W'^ < 4 GeV^. 



X interval 


< X > 


<Q2> 


An 
^2 














(GeV/c)2 


ibstat ibsyst 


±stat ±syst 


±stat ±syst 




0.029 - 


0.047 


0.038 


1.49 


0.143 ±0.105 ±0.023 


7.024 ± 4.777 ± 1.288 


6.963 ± 4.787 ± 1 


288 


0.047 - 


0.075 


0.060 


2.01 


-0.085 ± 0.065 ± 0.015 


-1.811 ± 1.574 ±0.367 


-1.948 ± 1.579 ±0 


367 


0.075 - 


0.120 


0.095 


2.60 


0.013 ±0.076 ±0.025 


0.254 ±0.868 ±0.266 


0.192 ±0.872 ±0 


266 


0.120 - 


0.193 


0.152 


3.21 


-0.013 ±0.114 ±0.040 


-0.118 ±0.541 ±0.180 


-0.079 ± 0.545 ± 


180 


0.193 - 


0.310 


0.241 


3.77 


0.067 ±0.192 ±0.056 


0.320 ± 0.342 ± 0.093 


0.263 ± 0.347 ± 


093 


0.310 - 


0.498 


0.378 


4.22 


-0.433 ± 0.419 ± 0.058 


-0.225 ±0.209 ±0.024 


-0.256 ± 0.214 ±0 


024 


0.498 - 


0.799 


0.595 


4.56 


0.926 ± 1.302 ±0.200 


0.116 ±0.086 ±0.010 


0.064 ± 0.091 ± 


010 


0.075 - 


0.120 


0.101 


3.77 


0.028 ±0.113 ±0.024 


0.315 ± 1.392 ±0.312 


0.184 ± 1.399 ±0 


312 


0.120 - 


0.193 


0.154 


4.97 


-0.078 ±0.091 ±0.029 


-0.407 ±0.528 ±0.172 


-0.384 ± 0.533 ± 


172 


0.193 - 


0.310 


0.242 


6.37 


-0.195 ±0.118 ±0.041 


-0.219 ±0.261 ±0.091 


-0.276 ± 0.265 ±0 


091 


0.310 - 


0.498 


0.381 


7.76 


-0.131 ±0.218 ±0.045 


-0.052 ±0.125 ±0.025 


-0.057 ± 0.128 ±0 


025 


0.498 - 


0.799 


0.584 


8.86 


0.544 ±0.627 ±0.120 


0.021 ± 0.041 ± 0.006 


0.023 ± 0.042 ± 


006 



TABLE XXXIII. Results for the moments vf^ and T^^^ evaluated at = 5 (GeV/c)^, and 
the extracted twist-3 matrix elements ^2 for proton (p), deuteron (d), and neutron(n). The errors 
include statistical (which dominate) and systematic contributions. 

02/2 = rf^ X lO'^ ^ X 10^ d2 X 10^ 

p 12.4 ± 1.0 -6.3 ± 1.8 5.8 ±5.0 

d 4.6 ±0.8 -1.4 ±3.0 5.1 ±9.2 

n -2.4 ±1.6 3.3 ±6.5 5.0 ± 21.0 



TABLE XXXIV. Theoretical predictions for the twist-3 matrix elements ^2 and for proton 
and d2 and /2 for deuteron. Also shown is fi, the higher twist correction to Ti described in the text. 



Bag models 
Ref. |g Ref. |6[ Ref. [^3| 



Ref. [p7| Ref. g| 



QCD sum rules 

Ref. m Ref. p9| 



Q2 (GeV/c)2 
£ X 103 



/I 



X 10^ 
X 10^ 
4 X 10^ 

/i^ X 10^ 



5 


5 


1 


1 


1 


1 




17.6 


6.0 


21 


10 


-6 ±3 


-3 ±3 








35 


28 


-37 ±6 


-50 ± 34 


-69 ±5 






27 


15 


-15 ±7 


-20 ± 13 


-27 ±2 


6.6 


2.9 


11 


5 


-17± 5 


-13± 5 








17 


14 


-25 ±4 


-34 ± 20 


-38 ±5 






13 


7 


-10 ±3 


-13 ±8 


-15 ±2 
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TABLE XXXV. Pion asymmetries versus momentum E' for proton and deuteron targets at 



E=29.1 GeV. 



E'(GeV) 




AP{Tr+) 




Ap+) 


e = 4.5° 


8.04 


0.028± 0.010 


0.016± 0.019 


O.OOli 0.011 


0.013± 0.021 


10.65 


0.019± 0.012 


0.043± 0.022 


-0.017 ± 0.013 


0.005± 0.025 


13.57 


0.019± 0.016 


0.019± 0.031 


-0.015 ± 0.018 


0.016± 0.036 


16.58 


O.OOlib 0.023 


0.065± 0.043 


-0.017± 0.025 


0.020± 0.052 


19.41 


0.002± 0.036 


0.041± 0.069 


0.108± 0.041 


0.031± 0.087 


21.88 


-0.044 ± 0.056 


0.028± 0.111 


-0.108± 0.068 


0.031± 0.146 


23.90 


0.018± 0.063 


0.273± 0.127 


-0.084± 0.079 


-0.172 ± 0.171 


25. U 


0.058± 0.082 


0.057± 0.168 


0.057± 0.111 


0.094± 0.242 


e = 7.0° 


5.61 


0.123± 0.069 


-0.010 ± 0.135 


0.022± 0.072 


-0.272 ± 0.195 


7.72 


O.OlSib 0.010 


0.027± 0.019 


-0.021 ± 0.010 


-0.038 ± 0.031 


10.29 


O.OOli 0.012 


0.081± 0.023 


-0.019 ± 0.013 


O.OOOi 0.037 


13.19 


0.016± 0.022 


0.057± 0.043 


0.063± 0.024 


0.034± 0.071 


16.19 


0.085± 0.036 


0.071± 0.075 


-0.021 ± 0.042 


0.142± 0.129 


19.06 


0.096± 0.052 


-0.035 ± 0.107 


0.070± 0.063 


0.231± 0.196 


21.59 


-0.034 ± 0.064 


0.106± 0.136 


0.023± 0.085 


0.060± 0.263 
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TABLE XXXVI. The measured virtual photon-nucleon asymmetry Ai + r]A2 and the spin 
structure function gi for the resonance region. The values of W"^ and are given at bin centers. 



The dilution factor / and applied correction term Arc, which for these data also includes a resolution 


correction, are from Eqs. 


I49I and [ 


^Fii T'Vip vi^liip of (ii 


in t,bp l?ist, mlnTTiTi is 


r?ilml?it,pr] froin A w nnrlpr 


the assumption that A2 


= 0. 












/ 


Arc 


Ai + 7]A2 
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GeV2 


(GeV/c)2 






± stat.ib syst. 


± stat.ib syst. 








Proton 


9 = 4.5° 




1.31 


0.55 


0.180 


-0.0444 


-0.086±0.126±0.225 


-0.011±0.017±0.030 


1.69 


0.54 


0.164 


-0.0284 


-0.453±0.125±0.127 


-0.112±0.031±0.032 


2.06 


0.53 


0.155 


0.0101 


0.461±0.106±0.104 


0.139±0.032±0.031 


2.44 


0.52 


0.153 


0.0204 


0.694±0.091±0.086 


0.344±0.045±0.035 


2.81 


0.50 


0.154 


-0.0051 


0.222±0.078±0.046 


0.142±0.050±0.027 


3.19 


0.49 


0.144 


-0.0049 


0.242±0.079±0.092 


0.142±0.046±0.050 


3.56 


0.48 


0.144 


-0.0047 


0.090±0.072±0.030 


0.061±0.049±0.018 


3.94 


0.47 


0.147 


-0.0032 


0.002±0.064±0.011 


0.002±0.054±0.010 


4.31 


0.46 


0.144 


-0.0013 


0.134±0.059±0.017 


0.129±0.057±0.008 


4.69 


0.45 


0.143 


-0.0013 


0.105±0.056±0.014 


0.110±0.058±0.009 








Deuteron 


9 = 4.5° 




1.31 


0.55 


0.247 


-0.0242 


-0.173±0.278±0.236 


-0.015±0.024±0.021 


1.69 


0.54 


0.235 


-0.0221 


-0.305±0.231±0.078 


-0.061±0.047±0.015 


2.06 


0.53 


0.232 


-0.0009 


0.290±0.157±0.095 


0.107±0.058±0.034 


2.44 


0.52 


0.227 


0.0115 


0.184±0.152±0.072 


0.078±0.064±0.030 


2.81 


0.50 


0.232 


0.0030 


0.021±0.128±0.055 


0.011±0.068±0.029 


3.19 


0.49 


0.230 


-0.0033 


0.147±0.113±0.021 


0.093±0.072±0.011 


3.56 


0.48 


0.230 


-0.0040 


0.023±0.102±0.011 


0.017±0.074±0.008 


3.94 


0.47 


0.230 


-0.0037 


-0.017±0.096±0.010 


-0.014±0.078±0.008 


4.31 


0.46 


0.233 


-0.0031 


-0.034±0.091±0.009 


-0.030±0.080±0.008 


4.69 


0.45 


0.229 


-0.0028 


0.055±0.086±0.010 


0.053±0.084±0.008 








Proton 


9 = 7.0° 




1.56 


1.26 


0.159 


-0.0524 


-0.143±0.128±0.113 


-0.015±0.013±0.012 


1.94 


1.23 


0.154 


-0.0012 


0.349±0.110±0.105 


0.038±0.012±0.012 


2.31 


1.20 


0.156 


0.0251 


0.795±0.087±0.088 


0.177±0.019±0.020 


2.69 


1.18 


0.158 


0.0036 


0.593±0.077±0.058 


0.163±0.021±0.014 


3.06 


1.15 


0.154 


0.0046 


0.507±0.069±0.078 


0.177±0.024±0.025 


3.44 


1.12 


0.150 


0.0029 


0.262±0.066±0.035 


0.101±0.025±0.010 


3.81 


1.10 


0.148 


0.0027 


0.299±0.063±0.039 


0.136±0.028±0.008 


4.19 


1.07 


0.148 


0.0038 


0.433±0.059±0.062 


0.221±0.030±0.012 


4.56 


1.04 


0.148 


0.0035 


0.324±0.056±0.054 


0.185±0.032±0.017 


4.94 


1.02 


0.147 


0.0033 


0.237±0.053±0.037 


0.150±0.034±0.010 








Deuteron 


9 = 7.0° 




1.56 


1.26 


0.232 


-0.0312 


0.043±0.243±0.103 


0.003±0.017±0.007 


1.94 


1.23 


0.232 


-0.0131 


0.163±0.179±0.128 


0.018±0.019±0.014 


2.31 


1.20 


0.230 


0.0108 


0.183±0.154±0.033 


0.030±0.025±0.005 


2.69 


1.18 


0.231 


0.0110 


0.237±0.134±0.072 


0.054±0.030±0.016 
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3.06 


1.15 


0.234 


0.0019 


0.389±0.114±0.046 


0.112±0.033±0.012 


3.44 


1.12 


0.228 


0.0000 


0.081±0.105±0.017 


0.028±0.037±0.006 


3.81 


1.10 


0.228 


-0.0005 


0.121±0.099±0.019 


0.049±0.040±0.006 


4.19 


1.07 


0.230 


-0.0002 


0.178±0.094±0.030 


0.080±0.042±0.009 


4.56 


1.04 


0.230 


-0.0006 


0.145±0.091±0.026 


0.073±0.046±0.008 


4.94 


1.02 


0.230 


-0.0006 


0.201±0.085±0.036 


0.113±0.048±0.012 
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TABLE XXXVII. Systematic errors (absolute) on Ai +r]A2 by category for the resonance region 
(See text for details). 





A 


model 


e 


E' 


Resol 




PbPt 


/ 


R 
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Proton 


e = 4.5° 










1.31 


0.0083 


0.1842 


0.0033 


0.1156 


0.0328 




0.0354 


0.0260 


0.0188 


0.0057 


1.69 


0.0006 


0.1245 


0.0051 


0.0163 


0.0052 




0.0019 


0.0014 


0.0169 


0.0046 


2.06 


0.0837 


0.0082 


0.0077 


0.0541 


0.0019 




0.0145 


0.0106 


0.0224 


0.0040 


2.44 


0.0548 


0.0614 


0.0070 


0.0114 


0.0098 




0.0211 


0.0155 


0.0520 


0.0088 


2.81 


0.0339 


0.0019 


0.0014 


0.0201 


0.0040 




0.0106 


0.0078 


0.0191 


0.0042 


3.19 


0.0544 


0.0832 


0.0008 


0.0118 


0.0027 




0.0106 


0.0078 


0.0349 


0.0040 


3.56 


0.0146 


0.0256 


0.0004 


0.0003 


0.0008 




0.0045 


0.0033 


0.0150 


0.0036 


3.94 


0.0102 


0.0015 


0.0000 


0.0022 


0.0020 




0.0009 


0.0006 


0.0003 


0.0032 


4.31 


0.0011 


0.0028 


0.0006 


0.0012 


0.0010 




0.0055 


0.0040 


0.0156 


0.0031 


4.69 


0.0062 


0.0017 


0.0004 


0.0008 


0.0001 




0.0043 


0.0032 


0.0115 


0.0029 










Deuteron 


9 = 


4.5° 








1.31 


0.1965 


0.1531 


0.0242 


0.1203 


0.0375 




0.0236 


0.0132 


0.0027 


0.0033 


1.69 


0.0638 


0.0613 


0.0129 


0.0331 


0.0048 




0.0015 


0.0008 


0.0259 


0.0049 


2.06 


0.0220 


0.0861 


0.0180 


0.0197 


0.0044 




0.0161 


0.0090 


0.0234 


0.0058 


2.44 


0.0107 


0.0702 


0.0059 


0.0073 


0.0010 




0.0046 


0.0026 


0.0126 


0.0040 


2.81 


0.0074 


0.0542 


0.0017 


0.0047 


0.0013 




0.0004 


0.0002 


0.0016 


0.0022 


3.19 


0.0110 


0.0047 


0.0067 


0.0046 


0.0016 




0.0099 


0.0056 


0.0108 


0.0034 


3.56 


0.0037 


0.0099 


0.0014 


0.0003 


0.0002 




0.0027 


0.0015 


0.0017 


0.0022 


3.94 


0.0029 


0.0094 


0.0003 


0.0001 


0.0002 




0.0003 


0.0002 


0.0016 


0.0019 


4.31 


0.0021 


0.0084 


0.0011 


0.0008 


0.0002 




0.0009 


0.0005 


0.0037 


0.0016 


4.69 


0.0024 


0.0056 


0.0027 


0.0007 


0.0001 




0039 


0.0022 


0.0060 


0.0016 










Proton 


e = 7.0° 










1.56 


0.0186 


0.1033 


0.0010 


0.0126 


0.0386 




0.0163 


0.0090 


0.0030 


0.0010 


1.94 


0.0165 


0.0186 


0.0086 


0.1006 


0.0010 




0.0163 


0.0089 


0.0137 


0.0019 


2.31 


0.0742 


0.0028 


0.0150 


0.0171 


0.0141 




0.0258 


0.0141 


0.0251 


0.0081 


2.69 


0.0103 


0.0287 


0.0088 


0.0196 


0.0012 




0.0237 


0.0130 


0.0370 


0.0064 


3.06 


0.0342 


0.0630 


0.0057 


0.0105 


0.0043 




0.0200 


0.0110 


0.0385 


0.0064 


3.44 


0.0106 


0.0196 


0.0029 


0.0039 


0.0024 




0.0100 


0.0055 


0.0262 


0.0033 


3.81 


0.0075 


0.0030 


0.0030 


0.0021 


0.0013 




0.0115 


0.0063 


0.0353 


0.0040 


4.19 


0.0054 


0.0056 


0.0040 


0.0004 


0.0004 




0.0164 


0.0090 


0.0591 


0.0061 


4.56 


0.0071 


0.0235 


0.0028 


0.0005 


0.0005 




0.0127 


0.0070 


0.0463 


0.0046 


4.94 


0.0009 


0.0099 


0.0018 


0.0004 


0.0003 




0.0089 


0.0049 


0.0341 


0.0039 










Deuteron 


9 = 


7.0 


p 








1.56 


0.0887 


0.0335 


0.0008 


0.0355 


0.0288 




0.0216 


0.0091 


0.0030 


0.0018 


1.94 


0.0144 


0.1146 


0.0172 


0.0523 


0.0004 




0.0158 


0.0066 


0.0016 


0.0019 


2.31 


0.0206 


0.0020 


0.0111 


0.0185 


0.0051 




0.0064 


0.0027 


0.0113 


0.0028 


2.69 


0.0072 


0.0686 


0.0095 


0.0040 


0.0004 




0.0098 


0.0041 


0.0163 


0.0043 


3.06 


0.0085 


0.0136 


0.0209 


0.0052 


0.0110 




0.0229 


0.0096 


0.0269 


0.0083 


3.44 


0.0045 


0.0129 


0.0051 


0.0033 


0.0032 




0.0049 


0.0020 


0.0056 


0.0031 


3.81 


0.0045 


0.0101 


0.0068 


0.0003 


0.0030 




0.0073 


0.0030 


0.0115 


0.0031 


4.19 


0.0046 


0.0095 


0.0101 


0.0062 


0.0004 




0.0113 


0.0047 


0.0232 


0.0045 
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4.56 0.0015 0.0075 0.0079 0.0032 0.0004 0.0091 0.0038 0.0208 0.0037 
4.94 0.0063 0.0059 0.0109 0.0082 0.0011 0.0130 0.0054 0.0289 0.0054 



TABLE XXXVIII. Integrals ri{Q^) of the structure functions gi for the proton (p), deuteron 
and neutron (n) at low Q^. Listed are the measured sums P^*^** for the resonance region {W'^ < 4 
GeV^) and F^^^ for the deep- inelastic region (W'^ > 4 GeV^ using data from 9.7 and 16 GeV beam 
energies), the low-x extrapolation Tf^^ for x < 0.03, and the combined total F^"*^. 



Q2 




pres 


rpis 


pext 


ptot 


(GeV/c)2 




ibstat.ibsyst. 


±stat.±syst. 




±stat.±syst. 


0.5 


Proton 


0.022 ± 0.007 ± 0.008 


0.017 ± 0.002 ± 0.003 


0.009 


0.047 ± 0.007 ± 0.015 


0.5 


Deuteron 


0.004 ± 0.010 ± 0.008 


0.004 ± 0.003 ± 0.002 


0.000 


0.008 ± 0.011 ± 0.009 


0.5 


Neutron 


-0.013 ±0.023 ±0.018 






-0.030 ± 0.024 ± 0.025 


0.5 


p-n 








0.077 ± 0.027 ± 0.036 


1.2 


Proton 


0.039 ± 0.003 ± 0.003 


0.051 ± 0.003 ± 0.003 


0.014 


0.104 ±0.005 ±0.016 


1.2 


Deuteron 


0.018 ± 0.005 ± 0.003 


0.019 ± 0.006 ± 0.002 


0.001 


0.037 ± 0.007 ± 0.006 


1.2 


Neutron 


-0.001 ± 0.010 ± 0.008 






-0.023 ± 0.016 ± 0.020 


1.2 


p-n 








0.127 ±0.018 ±0.034 
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FIG. 1. The layout of the lasers and the polarized electron gun at the accelerator injector is 
shown schematically. Two types of lasers are used, one for the SLC, which produces two 2 nsec 
pulses separated by 61 nsec, and one for the fixed target experiments, which produces one pulse 2 
/isec long. 
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FIG. 2. The layout of the M0ller polarimeter systems used in the EMS experiment (not to scale). 
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FIG. 3. (a) A typical M0llcr coincidence event time difference spectrum, (b)-(c) Two views of 
the distribution of coincidence events in the 7x7=49 possible combinations. View (c) shows the 
hidden back-side of view (b). True M0ller events were constrained to occur in the crest. 
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FIG. 4. Typical measured single-arm M0ller line-shapes in detector 7 at 29.1 GeV. The (R-l-L) 
data (points), (R-l-L) fit (histogram), and (R-l-L) background (shaded region) are plotted in (a) 
versus the relative y position of the channel. The (R-L) data (points) and (R-L) background 
(shaded region) are plotted in (b). 
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FIG. 5. Variation of the quantum efficiency (QE) of the polarized source over the course of the 
E143 experiment. 
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FIG. 6. Beam polarization versus QE for the the single (circles) and double-arm (squares) 
polarimeter systems. 
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FIG. 7. E143 target schematic. 
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FIG. 8. Target insert schematic. 
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FIG. 9. Typical TE signals measured from (a) polarized protons and (b) deuterons. The spikes 
in the deuteron signal are an artifact of the synthesized signal generator that was used. 
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FIG. 10. The polarization history over the course of a few days is shown as a function of received 
charge for (a) proton and (b) deuteron targets. 
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FIG. 11. A schematic of the E142/E143 spectrometer layout. 
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FIG. 12. Sample plot of events versus the ratio of energy to momentum as measured by the 
shower counter and tracking, respectively. The electron events are peaked at unity. The lower 
energy background pion events are removed when all cuts are applied to the data. 
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FIG. 13. g\/Ff as a function of for 8 different x bins. The data are from this experiment 
(soHd circle), SMC § (open circle), EMC @ (squares), SLAC E80 § (triangle), and SLAC E130 
1^] (diamond). The dashed and solid curves correspond to global fits II {gi/Fi Q^-independent) 
and III [gi/Fi Q^.dependent) in Table respectively. The E154 NLO pQCD fit |ll7) is shown 
as the dot-dashed curve. 
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FIG. 14. ^f/Ff* as a function of for 8 different x bins. The data are from this experiment 
(sohd circles) and SMC |M] (open circles). The dashed and solid curves correspond to global fits II 



{91/ Fi Q2-independent) and III {gi/Fi Q^.^ependent) in Table |Xl[ respectively. The E154 NLO 
pQCD fit |117| is shown as the dot-dashed curve. 
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FIG. 15. gi/Fi averaged over all beam energies and spectrometer angles as a function of x. The 
solid circles are 

at Q2 > 1 (GeV/c)2, while the diamonds are at < i (GeV/c)^ for x < 0.03. 
The data from SMC (squares) are shown for proton and deuteron, and the data from E154 
(open circles) are shown for neutron. The curves are NLO pQCD fits by Altarelli, Ball, Forte, and 
Ridolfi fS^ (solid) and Gliick, Reya, Stratmann, and Vogelsang |^ (dashed). 
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FIG. 16. xgi{x,Q^) evaluated at the average measured at each x. 
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FIG. 17. gi{x,Q'^) evaluated at = 3 (GeV/c)^ (assuming gi/Fi is independent of Q^) as a 
function of x. 
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FIG. 18. Measured singlet matrix elements a^'" from this experiment, E142 |T|], SMC and 
EMC/E80/E130 combined [^. These results were calculated from the published first moments 
of gi using up-to-date "invariant" singlet pQCD corrections and 6{3F — D) = 0.032. The dashed 
curve indicates the world average of 0.31 ± 0.04. 
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FIG. 19. The results for A2{x, Q ) as a function of x for this experiment (sohd squares for 4.5°, 



sohd circles for 7°). Also shown are SMC (open circles) |14|, E142 (open square) llOj], and E154 
(diamond) |jl5|. The solid and dashed curves correspond to the positivity constraint at 4.5° and 7° 
kinematics, respectively. Overlapping data have been shifted slightly in x to make errors clearly 
visible. The bands indicate systematic errors for the two E143 data sets. 
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FIG. 20. The results for xg2{x,Q'^) as a function of x for this experiment (squares for 4.5°, 
circles for 7°). Systematic errors are indicated by bands. Overlapping data have been shifted 
slightly in x to make errors clearly visible. The solid curve shows the twist-2 calculation for 

the kinematics of the 4.5° spectrometer. The same curve for 7° is nearly indistinguishable. The bag 
model calculations at = 5.0 (GeV/c)^ by Stratmann |^ (short dash) and Song and McCarthy 
(long dash) are indicated. The curves on the neutron plot are difficult to distinguish from zero. 
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FIG. 21. A\\ versus momentum measured from polarized protons for tt and 7r+, at £^ = 29 GeV 
in both spectrometers. 
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FIG. 22. 74|| versus momentum measured from polarized deuterons for vr and tt"*", at -E = 29 
GeV in both spectrometers. 
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FIG. 23. Extracted values for Ai + r]A2 (circles) in the resonance region for the proton at (a) 
4.5° and (b) 7.0°; and for the deuteron at (c) 4.5° and (d) 7.0°. Also shown are the Monte Carlo 
predictions (solid line) and the data of Baum et al. [|132| (diamonds). Error bars correspond to 
statistical errors only, whereas the bands below the data correspond to the systematic errors. 
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FIG. 24. Measurements of gi{x, Q^) in the resonance region as a function of W'^ for the proton 
at (a) 4.5° and (b) 7°; and for the deuteron at (c) 4.5° and (d) 7°. The present data (circles) are 
plotted together with the data of Baum et al. [|132| (diamonds) and our Monte Carlo simulation 
(solid line). The errors are indicated as in Fig. |2^. 
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FIG. 25. Integrals of gi{x,Q'^) at several fixed values of for (a) the neutron and (b) the 
proton. The present data (circles) are plotted together with data from E143 deep-inelastic results 
from this experiment (squares) and E142 (inverted triangle). The curves correspond to the 
evolution (3^ of the deep-inelastic results due to changing Us (solid) , the predictions of Burkert and 



loffe [|135| (short dash), the model of Soffer [136| (long dash), and the GDH approach to = 



(solid). The error bars are statistical, and the shaded bands correspond to the systematic errors. 
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